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ABSTRACT 


This  study  consists  of  two  parts.  The  first  part 
deals  with  a  laboratory  study  of  pattern  flooding.  A 
series  of  experiments  was  conducted  to  investigate  the 
effect  of  back  pressure  and  pattern  confinement  on  oil 
recovery  and  areal  sweep  efficiency. 

To  study  the  effect  of  back  pressure  a  normal 
five  spot  pattern  was  selected.  A  series  of  back  pressures 
ranging  from  3  cms .  to  70  cms .  of  water  column  was  studied. 
The  results  obtained  indicate  that  oil  recovery  and  areal 
sweep  efficiency  are  independent  of  back  pressure  within 
the  range  of  back  pressures  investigated. 

In  order  to  study  the  effect  of  pattern  confine¬ 
ment,  two  patterns:  a  direct  line  drive,  and  a  normal  five 
spot,  were  selected.  The  confinement  of  a  pattern  was 
achieved  by  surrounding  it  with  identical  patterns  around 
its  periphery  in  an  attempt  to  establish  isopotential 
boundaries  around  the  pattern  of  interest. 

The  direct  line  drive  pattern  was  surrounded  by 
two  rows  on  either  side  whereas  the  five  spot  pattern  was 
surrounded  by  eight  five  spots.  The  oil  recoveries  obtained 
in  these  two  cases  ranged  from  88.2%  for  a  direct  line 
drive  to  171.6%  in  case  of  the  five  spot. 

The  second  part  deals  with  an  examination  of  the 
Higgins-Leighton  method  for  predicting  the  performance  of 


■ 
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a  water  flood  in  an  oil  field.  Since  this  method  is  limi¬ 
ted  to  the  case  where  a  constant  pressure  drop  is  main¬ 
tained  between  the  injection  and  production  well,  it  was 
modified  and  extended  to  cover  the  following  cases: 

-  constant  injection  rate 

-  varying  injection  rate  and  pressure  drop. 

The  original  and  the  modified  Higgins-Leighton 

methods  were  compared  with  a  laboratory  study  of  a  five 
spot  pattern  reported  in  the  literature.  The  agreement 
between  the  experimental  and  the  calculated  results  was 


found  to  be  better  in  the  case  ofthe  modified  approach. 
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PART  ONE 

THE  EFFECT  OF  PATTERN  CONFINEMENT 

ON  OIL  RECOVERY  AND  AREAL  SWEEP 

EFFICIENCY 
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INTRODUCTION 

A  full  utilization  of  pilot  test  results  and 
reliable  prediction  of  large  scale  water1  floods  requires 
a  knowledge  of  the  relationship  between  the  recovery  per¬ 
formance  of  confined  and  unconfined  well  patterns.  It 
was  the  objective  of  this  investigation  to  examine  the 
effect  of  pattern  confinement  on  oil  recovery  and  areal 
sweep  efficiency  in  the  case  of  a  five  spot  pattern. 

To  date  only  limited  basic  information  is  avail¬ 
able  on  the  behavior  of  confined  pilot  water  floods.  Dalton, 
Rapoport  and  Carpenter  (1)  investigated  the  effect  of  the 
ratio  of  pressure  drawdown  at  the  producing  well  to  the 
pressure  build  up  at  the  injection  well  (it-  ratio),  on  oil 
recovery  in  confined  and  unconfined  pattern  floods.  Caudle 
and  Loncaric  (2)  on  the  other  hand  studied  the  effects  of 
changing  the  ratio  of  injection  rate  to  withdrawal  rate  on 
oil  recovery  for  confined  and  unconfined  five  spot  pilot 
floods.  However,  hardly  any  of  the  work  to  date  has  dealt 
explicitly  with  oil  recovery  and  areal  sweep  performance  of 
confined  pilot  floods. 

A  considerable  amount  of  experimental  work  in  this 
direction  has  been  conducted  at  the  University  of  Alberta. 
This  work  was  reviewed.  It  was  discovered  that  the  earlier 
workers  had  difficulties  with  regard  to  model  expansion, 
partial  well  penetration  and  cleaning  of  the  model.  Since  a 


2 


larger  number  of  wells  was  needed  to  achieve  confinement 

in  the  present  work  and  the  problems  encountered  in  the 
previous  work  were  to  be  avoided,  it  was  decided  to  con¬ 
struct  a  new  model. 

Finally,  it  was  decided  to  reinvestigate  the 
effect  of  back  pressure  on  oil  recovery  and  areal  sweep 
efficiency  in  order  to  resolve  the  controversy  expressed 
by  John  W.  Serra  (3)  and  S.K,  Bhatia  (4). 
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LITERATURE  SURVEY 


The  recovery  of  oil  from  a  pilot  flood  is  gov¬ 
erned  by  the  following  factors: 

1.  Amount  of  oil  originally  present  in  the  reservoir. 

2.  Rate  of  water  injection. 

3.  Viscosity  of  oil,  injection  water  and  connate 
water . 


4. 

5. 

6. 


Mobility  ratio  M  = 


K  y  - 

A  displacing  _  Aw  _  w  ° 


A  displaced 
The  original  free  gas  saturation. 


Ao 


”  K 
w  o 


Interfacial  tension  between  oil  and  water 


7.  Pore  size  distribution  which  together  with  inter¬ 
facial  tension  affects  the  capillary  pressure 
characteristics  of  the  rock. 

8.  Wettability  of  the  rock. 

9.  Gravity  effects. 

10.  Pattern  configuration. 

A  comprehensive  literature  review  covering  the 
effects  of  these  parameters  on  oil  recovery  may  be  found 
in  the  work  of  Neilson  (5),  Pritchard  (6),  Culham  (7), 

Serra  (8)  and  Bhatia  (9). 

Neilson (10)  investigated  the  effect  of  a  pre- 
established  free  gas  saturation  on  the  sweep ef ficiency 
of  an  isolated  inverted  five  spot  and  concluded  that  the 
ultimate  recovery  may  be  slightly  increased  by  establishing 
a  free  gas  saturation  prior  to  water  flooding.  He  also 
pointed  out  that  the  producing  water-oil  ratio  for  an 


4 


isolated  pattern  increases  very  slowly  after  breakthrough 
so  that  ultimate  areal  sweep  efficiencies  six  times  that 
at  breakthrough  are  possible. 

Culham  (11)  conducted  a  study  of  areal  sweep 
efficiency  as  a  function  of  water  injected,  for  a  wide 
range  of  mobility  ratios,  in  the  case  of  a  nine  spot  pat¬ 
tern  and  concluded  that  an  increase  in  mobility  ratio 
causes  the  areal  sweep  efficiency  and  oil  recovery  to  de¬ 
crease.  He  also  studied  the  effect  of  injection  rate  on 
the  performance  of  a  water-wet  two  dimensional  system  and 
was  of  the  opinion  that  oil  recovery  decreases  with  increa¬ 
sing  rate  until  stabilized  conditions  are  obtained. 

Pritchard's  (12)  work  was  concerned  with  obtaining 
more  information  on  oil  recovery  after  breakthrough  in  the 
case  of  an  isolated  nine  spot  pattern.  During  the  course 
of  his  investigation,  the  effect  of  injection  rate  on  oil 
recovery  was  also  studied.  The  results  obtained  from  the 
experimental  work  indicated  that  the  critical  rate  at  which 
oil  recovery  is  rate  independent  did  not  exist  for  the  sy¬ 
stem  under  study;  instead  an  optimum  rate  was  obtained  such 
that  the  oil  recovery  decreased  at  injection  rates  higher 
than  this  optimum  rate. 

Serra  (13)  investigated  the  effect  of  back  pressure 
on  the  production  performance  and  sweep-out  patterns  for 
the  case  of  a  single  normal  five  spot.  It  was  concluded 
from  the  experimental  results  that  the  area  contacted  by 
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the  displacing  fluid  is  sharply  reduced  as  the  back  pres¬ 
sure  on  the  system  is  increased.  He  also  pointed  out  that 
a  considerable  amount  of  production  could  come  from  the 
area  beyond  the  well  pattern;  hence,  a  single  isolated  pat¬ 
tern  may  not  give  accurate  estimates  of  the  expected  per¬ 
formance  of  a  fully  developed  flood. 

Bhatia  (14)  studied  the  effect  of  back  pressure 
on  oil  recovery  and  areal  sweep  efficiency  in  the  case  of 
the  isolated  normal  five  spot  and  four  inverted  five  spot 
patterns.  It  was  concluded,  on  the  basis  of  experimental 
results,  that  back  pressure  affects  neither  oil  recovery 
nor  areal  sweep  efficiency.  Bhatia  (15)  was  of  the  opinion 
that  the  effect  of  back  pressure  measured  by  Serra  was  due 
to  the  expansion  in  the  experimental  model. 
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EQUIPMENT  DESIGN 

After  reviewing  the  work  of  Serra  (16)  and 
Bhatia  (17) ,  it  was  decided  to  retain  the  major  features 
of  their  design.  That  is,  the  model  should  consist  of  a 
1/4  inch  thick  bead  pack  sandwiched  between  two  lucite 
sheets.  However,  if  the  patterns  were  to  be  confined 
a  greater  number  of  wells  would  be  required  and,  therefore, 
the  size  of  the  model  would  have  to  be  increased.  In  addi¬ 
tion;  the  wells  were  used  to  fasten  the  confining  lucite 
sheets  to  prevent  model  bulging. 


DESIGN  OF  OILFIELD  MODEL 

The  size  of  the  model  was  dictated  by  the 
availability  of  material  rather  than  engineering  design. 
Since  the  only  readily  available  thick  lucite  sheet  was 
4  feet  square  X  2  inches  thick,  it  was  decided  to  build 
the  model  of  that  material  rather  than  to  delay  the  pro¬ 
ject.  Consequently,  the  design  consisted  of  two,  two  inches 
thick  lucite  discs  having  a  diameter  of  four  feet. 

A  total  of  16,  1/8  inch,  holes  was  drilled  on 
a  circle  of  3  feet  11  1/4  inches  diameter  to  accodate  the 
wells  used  for  cleaning  the  model.  In  order  to  facilitate 
a  study  of  various  patterns,  a  total  of  185,  1/8  inch, 
holes  was  drilled  on  a  square  pattern  in  the  central  area 


7 


of  these  sheets.  These  holes  were  drilled  in  both  the 
sheets  simultaneously  to  avoid  misalignment  problems  while 
inserting  the  wells. 

These  two  lucite  sheets  were  separated  by  a 
1/4  inch  lucite  washer  having  an  outside  diameter  of  4  feet 
and  an  inside  diameter  of  3  feet  11  5/8  inches.  This 
assembly  was  then  fastened  together  around  the  periphery 
by  means  of  72,  1/4  inch,  brass  bolts.  A  seal  was  effected 
between  the  three  lucite  pieces  by  embedding  a  neoprene 
'o'  ring  in  each  of  the  thick  lucite  discs  opposite  the 
inner  washer. 

A  detailed  drawing  of  the  oilfield  model  design 
is  presented  in  Figure  I. 

DESIGN  OF  OIL  WELLS 

The  wells  were  designed  so  as  to  accomplish 
the  following  objectives: 

-  to  provide  a  facility  for  injecting  and  with¬ 
drawing  fluids  from  the  model,  and 

-  to  fasten  the  lucite  sheets  together  in  such 

a  manner  as  to  avoid  the  model  -  expansion  pro¬ 
blem. 

The  wells  consisted  of  two  parts: 

1.  Well  Head 


2.  Well  Body 


OIL  FIELD  MODEL  DESIGF- 


SECTION 
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1 .  Well  Head 

The  well  head  was  constructed  from  a  7/8  inch 
hexagonal  brass  stock.  This  piece  of  material  was 
drilled  on  one  end  to  accept  a  1/8  inch  N.P.T.  valve 
and  on  the  other  end  to  accept  a  1/8  inch  diameter  brass 
tube . 

2 .  Well  Body 

The  well  body  was  constructed  from  a  brass  tube 
1/8  inch  in  diameter  and  2  1/4  inches  long.  One  end 
of  this  tube  was  soldered  into  the  well  head  and  the 
other  extended  into  the  model  through  the  bead  pack. 

The  section  of  this  tube  exposed  within  the  pack  was 
drilled  with  4,  1/32  inch  diameter  holes,  disposed  on 
two  different  planes  and  situated  around  the  circum¬ 
ference  of  the  tube.  A  solid  1/8  inch  brass  rod  was 
soldered  onto  the  bottom  of  the  tube.  The  rod  extended 
through  the  lower  lucite  plate  and  its  protruding  end 
was  provided  with  #5-40  threads. 

The  entire  well  was  sealed  with  'o'  rings,  one 
at  the  top,  and  another  next  to  the  nut  and  washer  at 
the  bottom.  A  detailed  sectional  drawing  of  a  typical 
well  is  presented  in  Figure  II. 
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DESIGN  OF  INJECTION  SYSTEM 

In  order  to  permit  the  confinement  of  a  five  spot 
pattern,  it  was  felt  that  a  minimum  of  sixteen  injection 
lines  would  be  required.  The  injection  system  was,  there¬ 
fore,  designed  to  fulfill  this  condition.  This  was  achieved 
by  arranging  seventeen  stainless  steel  cylinders  according 
to  the  pattern  presented  in  Figure  III.  These  cylinders  were 

t 

bolted  to  a  bottom  plate. 

Of  these  seventeen  cylinder-piston  assemblies  the 
one  positioned  in  the  center  of  the  array  was  designed  to 
be  double  acting.  This  design  allowed  it  to  be  used  to  raise 
and  lower  the  top  plate, thereby  driving  the  pistons  in  the 
surrounding  cylinders  in  unison. 

Since  a  constant  injection  rate  was  desired,  the 
double  acting  master  cylinder  was  driven  by  means  of  a 
Ruska  proportioning  pump  (Cat.  No.  2249,  W  II,  Ser.  No. 

7681) .  This  pump  was  capable  of  providing  a  total  of  28 
constant  injection  rates  as  presented  in  Table  -  I. 

PRESSURE  TESTING  OF  THE  OILFIELD  MODEL 

Before  packing  the  model  with  glass  beads,  it  was 
necessary  to  pressure  test  it  in  order  to  ensure  that  no 
leakage  or  well  breakage  would  occur  during  the  experimen¬ 
tal  runs.  In  order  to  conduct  the  pressure  test,  the  wells 
were  mounted  on  the  model  and  held  in  position  by  tightening 
the  nuts  provided  at  their  base.  These  nuts  were  tightened 
very  carefully  to  avoid  the  development  of  an  uneven  stress 
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distribution.  The  model  was  then  completely  filled  with 
water  and  pressure  tested  using  the  Ruska  pump  mentioned 
earlier.  The  pressure  in  the  model  was  observed  on  a  mano¬ 
meter  connected  to  one  of  the  peripheral  wells. 

The  test  was  conducted  by  raising  the  test  pres¬ 
sure  step-wise.  If  and  when  various  well  assemblies  failed 
they  were  removed,  repaired,  and  reinstalled.  This  proce¬ 
dure  was  continued  until  a  test  pressure  of  9  p.s.i.  could 
be  maintained  without  leakage  or  mechanical  failure.  This 
pressure  of  9  p.s.i.  was  selected  on  the  basis  of  the  fol¬ 
lowing  considerations: 

1.  Pressure  drop  due  to  flow  in  the  old  model. 

2.  Pressure  due  to  a  column  of  70  cms .  of  water 
(maximum  back  pressure  imposed  on  the  model) . 

3.  A  safety  factor  of  2.0. 

PACKING  OF  THE  MODEL 

U.S.  30-40  industrial  glass  beads  were  used  to  pack 
the  model.  The  beads  were  first  placed  in  chromic  acid  for 
3  days ,  then  washed  with  an  aqueous  solution  of  sodium  hy¬ 
droxide  ,  and  finally  washed  repeatedly  with  distilled  water 

until  the  wash  water  was  neutral  to  litmus  paper.  These 
beads  were  then  dried  and  sieved  to  remove  clusters. 

The  model  was  thoroughly  cleaned  and  dried  by  pas¬ 
sing  air  through  it.  It  was  then  clamped  in  the  vertical 
position.  Two  vibrators  [Dynamite  models  222  and  223,  cap¬ 
acity  10  watts  each,  frequency  60  cycles/second]  were  con- 
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nected  to  the  model.  The  beads  were  weighed  and  poured 
into  the  model  in  batches.  Each  batch  was  allowed  to  set 
tie  for  3-4  hours  before  the  next  one  followed.  Just  be¬ 
fore  complete  fill  up,  the  model  was  tapped  gently  with  a 
lucite  hammer  to  further  promote  packing.  The  entire  pro 
cess  took  about  seven  days  and  resulted  in  a  homogeneous, 
uniform  pack. 
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EXPERIMENTAL  WORK 

GENERAL  PROCEDURE 

The  model  was  cleaned  and  brought  to  an  initial  wa¬ 
ter  saturation  condition,  using  the  peripheral  wells.  The 
water  used  for  injection  purposes  was  drawn  into  the  cylinder 
assembly,  from  an  overhead  water  reservoir  through  a  system 
of  quick  opening  toggle  valves.  The  water  injection  rate 
was  fixed  by  choosing  the  appropriate  gear  position  on  the 
Ruska  pump. 

The  water  injected  into  the  pattern  under  study  con¬ 
tained  a  fluorescent  dye.  The  flood  pattern  was  traced  by 
causing  the  dyed  water  to  fluoresce.  This  was  accomplished 
by  placing  an  ultra-violet  light  source  underneath  the  mo¬ 
del.  A  Pentacon-Six  Camera  was  used  to  take  photographs 
of  the  flood  front  from  time  to  time.  This  camera  used  a 
tri-X  120  black  and  white  film,  a  distance  of  5.5  ft.,  lens 
opening  between  2.8-4  and  exposure  time  of  5  seconds. 

The  injection  pressure  was  measured  using  mercury 
manometers  with  a  short  and  a  long  column  of  dyed  water. 

The  water  containing  the  dye  was  made  to  fluoresce  by  two 
ultra-violet  tubes  mounted  on  an  iron  stand  and  facing  the 
manometer  board.  These  manometers  were  photographed  using 
a  Rolleicord  Camera,  with  tri-X  120  black  and  white  film, 
aperture  between  2.8-4,  and  an  exposure  time  of  10  seconds. 

A  series  of  photographs  was  taken  during  each  experimental 
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run.  The  pressure  drops  were  then  estimated  by  measuring 
the  difference  between  colored  limbs,  using  a  travelling 
microscope.  A  schematic  diagram  of  the  equipment  is  pre¬ 
sented  in  Figure  IV  .  Figure  V  presents  a  photograph  of 
the  equipment,  showing  the  layout  of  various  units. 

PRELIMINARY  EXPERIMENTS  ON  THE  MODEL 

A  number  of  preliminary  experiments  was  conducted 
on  the  model  to  determine  the  porosity  and  absolute  permeability 
of  the  pack.  The  calculation  procedure  used  and  the  re¬ 
sults  obtained  are  presented  in  Appendix  A. 

POROSITY 

In  determining  porosity,  the  bulk  volume  of  the 
model  was  determined  by  carefully  measuring  the  amount  of 
water  necessary  to  fill  the  model,  prior  to  packing.  Know¬ 
ing  the  weight  of  the  beads  contained  in  the  pack  and  their 
specific  gravity,  it  was  possible  to  calculate  the  matrix 
volume  of  the  model.  The  bulk  volume  and  the  matrix  volume 
were  then  used  to  calculate  porosity.  Alternatively  pore 
volume  was  determined  by  measuring  the  volume  of  water  re¬ 
quired  to  completely  saturate  the  packed  model.  Both  me¬ 
thods  yielded  the  same  value  for  porosity,  namely  36.43%. 
ABSOLUTE  PERMEABILITY 

Since  the  model  was  completely  saturated  with  water, 
a  five  spot  pattern  was  selected  and  water  injected  into 
the  model  through  the  four  injectors.  The  pressure  drop 
due  to  flow  of  water  was  measured  and  the  absolute  permea- 
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bility  of  the  model  calculated  using  Muskat's  (18)  steady 
state,  radial  equation  for  a  five  spot  pattern.  The  absolute 
permeability  of  the  model  was  found  to  be  27.7  darcys. 

FLUID  PROPERTIES 

In  the  present  investigation  only  two  phases  viz  •  , 
oil  and  water  were  considered.  The  following  properties 
for  these  fluids  were  determined. 

1.  Density  -  using  a  specific  gravity  bottle. 

2.  Viscosity  -  using  a  Cannon  -  Fenske  Ostwald  Visco¬ 
meter  . 

3.  Interfacial  Tension  -  using  a  Dii-Nuoy  Tensiometer. 
The  values  obtained  at  a  temperature  of  78 °F  are 

presented  in  Table  II. 


TABLE  II 


FLUID  PROPERTIES  AT  78°F 


FLUID 


DENSITY 


Kerosene 
(Esso,  Regu¬ 
lar  grade) 

Distilled 

Water 


gms . /c . c . 
0 .7297 


0 .9900 


VISCOSITY 
cps  . 

1.311 


0.9680 


WATER- OIL 

INTERFACIAL  TENSION 
dynes ./cm. 


32.140 


Dyed  Water  0.9908 


0 .9700 


Injection  water  used  to  displace  the  oil  was  deae¬ 
rated  distilled  water  colored  with  a  fluorescent  dye  (Sod¬ 
ium  Fluorescein  ,  Uranine  B,  Fischer  Scientific).  This  dye 
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was  completely  soluble  in  the  water  phase  but  insoluble  in 
the  oil  phase.  A  concentrate  of  the  dye  was  made  with  5  gms . 
of  dye  per  500  c.cs.  of  distilled  water  and  250  c.cs.  of 
such  concentrate  was  mixed  in  8.25  litres  of  distilled  wa¬ 
ter  for  injection  purposes. 

SCALING 

The  model  was  scaled  so  that  the  performance  was 
not  rate  sensitive.  This  was  accomplished  by  using  the 
Rapoport,  Carpenter  and  Leas  (19)  scaling  coefficient  as 
the  criterion. 

The  critical  injection  rate  was  calculated  for 
the  system  under  study  and  the  injection  rates  used  in  var¬ 
ious  experimental  runs  were  kept  above  this  critical  value. 
The  calculations  for  the  critical  injection  rate  are  pre¬ 
sented  in  Appendix  B. 

PATTERN  CONFIGURATION 

The  following  pattern  configurations  were  selected 
for  conducting  the  experimental  runs. 

1.  Isolated  Normal  Five  Spot. 

2.  Confined  Normal  Five  Spot. 

3.  Confined  Direct  Line  Drive. 

The  arrangement  of  injectors  and  producers  in  each 
of  these  patterns  is  presented  in  Figure  VI. 

EXPERIMENTAL  PROCEDURE 

The  experimental  runs  were  carried  out  using  the 


following  steps: 


ARRANGEMENT  OF  WELLS  IN  VARIOUS  PATTERNS  STUDIED 
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1.  In  order  to  establish  connate  water  saturation,  the 
model  was  saturated  with  water  and  clamped  in  the 
vertical  position.  The  oil  was  injected  into  the 
model  through  the  two  uppermost  peripheral  wells, 
while  water  was  removed  from  one  of  the  lowermost 
wells.  The  water  produced  was  collected  in  a  ves¬ 
sel  kept  at  a  height  higher  than  the  injection  point 
in  the  model.  This  procedure  helped  in  taking  full 
advantage  of  the  density  difference  between  oil 

and  water,  thus  maintaining  the  displacing  front  as 
nearly  horizontal  as  possible.  The  injection  of 
the  oil  was  continued  until  no  more  water  was  pro¬ 
duced.  The  model  was  then  clamped  in  the  horizontal 
position . 

2.  The  flood  pattern  was  selected  and  water  was  removed 
from  the  injection  and  the  production  well  stems 

by  injecting  oil  into  the  model  through  one  of  the 
peripheral  wells.  The  production  lines  were  filled 
with  oil  before  starting  the  run. 

3.  Back  pressure  on  the  model  was  fixed  by  raising  or 
lowering  the  production  manifold  to  the  desired  le¬ 
vel.  The  elevation  of  the  manifold  above  the  center 
of  the  bead  pack  was  measured  by  means  of  a  catheto- 
meter . 

4.  The  injection  and  the  production  wells  were  opened 
and  sufficient  time  was  allowed  for  the  system 
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to  attain  pressure  equilibrium.  The  initial  pres¬ 
sure  drop  in  the  manometer  was  recorded.  The  gear 
position  on  the  Ruska  pump  was  adjusted  to  deliver 
the  desired  injection  rate. 

5.  The  pump  was  started  and  the  production  from  the 
model  was  collected  in  a  tube  bank  arranged  under¬ 
neath  the  production  manifold. 

6.  The  flood  front  was  traced  from  time  to  time,  by 
photographing  the  model.  The  volumes  of  fluids 
produced  at  each  instant  that  the  model  was  photo¬ 
graphed  were  noted.  The  experimental  run  was  con¬ 
tinued  until  boundary  effects  tended  to  distort 
pattern  growth. 

7.  After  the  run  was  over,  the  model  was  clamped  in 
the  vertical  position.  The  five  lowermost  periph¬ 
eral  wells  were  connected  to  the  injection  lines 
and  water  was  delivered  to  these  lines  from  the 
overhead  water  reservoirs.  The  oil  was  removed  from 
the  model  using  the  three  uppermost  peripheral  wells 
and  was  collected.  This  process  of  flushing  the 
model  with  water  was  continued  until  no  more  oil 
was  produced  and  the  model  was  completely  free  of 
dye . 

8.  Steps  1  through  7  were  repeated  for  the  next  run. 


.  . 
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The  photographs  obtained  during  the  run  were  plani- 
metered  to  find  the  area  of  the  porous  medium  contacted  by 
the  flooding  water  at  any  stage.  This  allowed  the  calcula¬ 
tion  of  areal  sweep  efficiency. 
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EXPERIMENTAL  RESULTS  AND  DISCUSSION 
MODEL  EXPANSION  STUDY 

A  preliminary  study  was  conducted  on  Bhatia's  model 
to  determine  the  validity  of  his  contention  (20)  that  the 
model  expanded  during  the  course  of  experimentation.  To  do 
this,  the  model  was  saturated  with  water  and  three  pointers 
were  attached  to  the  wells  along  a  diameter  on  the  top 
surface  of  the  model.  Similarly,  three  pointers  were  at¬ 
tached  to  the  bottom  surface  of  the  model,  directly  under¬ 
neath  the  pointers  on  the  top  surface.  The  injection  lines 
were  then  connected  to  the  model  in  the  same  five  spot 
pattern,  as  used  by  Bhatia  and  Serra.  Prior  to  starting 
the  run,  the  initial  elevations  of  the  pointers  were  mea¬ 
sured  using  a  cathetometer .  The  injection  rate  was  varied 
from  400  ccs/hr  to  960  ccs/hr  in  steps  and  at  each  step, 
the  elevation  of  the  pointers  was  recorded.  The  results 
obtained  are  tabulated  in  Table  III,  Appendix  C. 

The  net  expansion-contraction  curve  shown  in  Fig¬ 
ure  VII  indicates  that  the  net  expansion  in  the  model  in¬ 
creases  with  the  injection  rate  up  to  a  rate  of  640  ccs/hr, 
beyond  which  point  it  decreases.  It  is  interesting  to  note 
that  one  part  of  the  model  expands  while  the  other  contracts . 
This  can  perhaps  be  explained  on  the  basis  of  irregular 
stress  distribution  in  the  model.  A  similar  experiment 
was  conducted  on  the  new  model  and  no  measureable  expansion 
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or  contraction  could  be  detected.  Another  experiment  in 
this  direction  measured  the  .increase  in  pore  volume  at 
various  back  pressures.  The  results  of  this  experiment 
are  shown  in  Figure  VIII  and  tabulated  in  Table  IV,  Ap¬ 
pendix  C.  It  may  be  observed  that  3.05  c.cs  of  oil  was 
injected  into  the  model  at  a  back  pressure  of  70.0  cms . 
of  water.  This  represents  a  net  increase  of  0.0003  cms. 
in  the  thickness  of  the  model  and  is  beyond  the  accuracy 
of  the  cathetometer . 

An  alternative  approach  to  measure  expansion  in 
the  model  is  by  measuring  pressure  drop  at  various  injec¬ 
tion  rates,  with  back  pressure  as  a  parameter.  A  radial 
flow  experiment  was  set  up  for  this  purpose.  The  water 
was  injected  into  the  model  from  the  peripheral  wells  and 
the  production  obtained  through  the  central  well.  The 
results  of  this  experiment  are  tabulated  in  Table  V,  Ap¬ 
pendix  C,  and  plotted  in  Figure  IX.  The  linear  relation¬ 
ship  between  pressure  drop  and  flow  rate  indicates: 

—  that  laminar  flow  persists  in  the  porous  medium  at 
the  flow  rates  studied, 

—that  there  is  no  expansion  in  the  model. 

During  the  radial  flow  experiment,  colored  water  was 
injected  into  the  four  peripheral  wells  located  on  the  prin¬ 
cipal  axes  of  the  model.  Figure  X  presents  a  pictorial 
representation  of  the  pattern  development  from  start  to 
the  completely  developed  stage.  It  was  concluded  from 
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these  photographs  that  the  packing  in  the  model  was  uniform 
and  homogeneous . 

EFFECT  OF  BACK  PRESSURE 

As  mentioned  earlier,  the  oil  recovery  and  areal 
sweep  efficiency  of  a  pattern  flood  is  influenced  by  the 
following  factors. 

1)  Gravity  effects. 

2)  Capillary  effects. 

3)  Injection  rate. 

4)  Viscosity  of  oil,  connate  water  and  injection  water. 

5)  Wettability  of  the  porous  system. 

6)  Pattern  configuration. 

7)  Mobility  ratio. 

8)  Interfacial  tension  between  oil  and  water. 

9)  Amount  of  oil  originally  present  in  the  reservoir. 
Therefore,  in  order  to  investigate  the  effect  of  back 

pressure,  it  was  important  to  ensure  that  all  the  factors 
listed  above  remained  constant  in  the  experimental  runs. 
Since  the  thickness  of  the  porous  medium  used  to  carry 
out  this  investigation  was  1/4  inch,  it  could  be  safely 
assumed  that  gravitational  effects  were  negligible.  The 
effects  due  to  capillarity  and  rate  were  held  constant  by 
selecting  an  injection  rate  higher  than  the  critical,  cal¬ 
culated  from  the  scaling  coefficient  suggested  by  Rapo- 
port,  Carpenter,  and  Leas  (21) .  Parameters 
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4-8  were  held  constant  by  using  the  same  oil,  water  and 
pattern  configuration  in  the  various  experimental  runs. 

The  effect  of  back  pressure  was  investigated  by 
conducting  a  series  of  experimental  runs  (Runs  5-10). 

An  isolated  five  spot  pattern  was  selected  and  the  results 
obtained  are  tabulated  in  Tables  VI  -  XI,  Appendix  C. 

The  back  pressures  studied,  ranged  from  3.0  to  70.0  cms . 
of  water.  Due  to  physical  limitations,  higher  back  pres¬ 
sures  could  not  be  considered. 

In  order  to  establish  that  parameter  9  remained 
constant,  a  check  was  made  on  the  reproducibility  of 
experimental  results.  The  experimental  runs  8  and  10  were 
conducted  under  exactly  identical  conditions.  The  results 
obtained  are  plotted  in  Figure  XI.  It  may  be  observed 
that,  within  experimental  error  ,  the  results  of  Runs  8 
and  10  are  reproducible. 

The  results  of  experimental  runs  5-9  were  analysed 
by  plotting  oil  recovery  and  areal  sweep  efficiency  against 
cumulative  throughput  with  back  pressure  as  a  parameter. 
Figure  XII  presents  a  plot  of  cumulative  throughput  versus 
oil  recovery  for  five  different  back  pressures.  It  may 
be  noticed  that  the  scatter  in  the  data  points  in  Figure 
XII  becomes  pronounced  after  9.0  network  pore  volumes  of 
water  ha/e  been  injected.  This  is  perhaps  due  to  the  in¬ 
fluence  of  the  boundary  on  the  growth  of  the  pattern.  Thus 
it  may  be  concluded  that  within  experimental  error  ,  oil 
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recovery  is  independent  of  back  pressure. 

Figure  XIII  presents  a  variation  of  areal  sweep 
efficiency  as  a  function  of  cumulative  throughput  at  var¬ 
ious  back  pressures.  The  scatter  in  this  figure  is  not 
so  pronounced  which  is  due  to  the  fact  that  data  beyond 
10.0  network  pore  volumes  of  water  injected,  were  not 
considered.  At  a  cumulative  throughput  higher  than  10.0 
network  pore  volumes  there  was  a  significant  distortion 
in  the  pattern  with  the  result  that  it  was  difficult  to 
obtain  a  meaningful  value  for  the  area  swept.  It  may  thus 
be  concluded  that  within  experimental  error  ,  areal  sweep 
efficiency  is  not  influenced  by  the  pressure. 

Figure  XIV  shows  the  variation  of  the  instantan¬ 
eous  water-oil  ratio  as  a  function  of  cumulative  through¬ 
put.  The  average  curve  plotted  through  the  data  points 
shows  a  tendency  to  stabilize  at  an  instantaneous  water 
oil  ratio  of  3.5  and  a  cumulative  throughput  of  12.4  net¬ 
work  pore  volumes.  This  is  due  to  the  fact  that  there  is 
a  considerable  oversweep  in  the  case  of  an  isolated  pat¬ 
tern.  This  phenomenon  of  over sweep  causes  the  production 
of  additional  oil  from  freshly  contacted  areas,  thereby 
resulting  in  only  a  slight  increase  in  the  water-oil 
ratio.  However,  if  the  experimental  runs  were  continued 
for  a  longer  period  of  time,  the  oversweep  would  be  con¬ 
trolled  by  the  physical  dimensions  of  the  model  and  a 
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steady  increase  in  water-oil  ratio  could  be  noticed. 

An  example  of  the  photographs  taken  during  a  run 
showing  the  area  contacted  by  the  injected  water,  for  an 
isolated  five  spot  and  confined  five  spot  patterns  is 
shown  in  Figures  XV  -  XIX.  The  white  portion  in  these 
photographs  shows  the  area  contacted  by  the  water. 

EFFECT  OF  PATTERN  CONFINEMENT 

To  estimate  the  extent  of  oversweep  in  an  isolated 
pattern,  several  experimental  runs  (Runs  11,  12)  were  con¬ 
ducted  on  the  confined  five  spot  and  confined  direct  line 
drive  patterns.  The  data  obtained  are  reported  in  Tables 
XII  and  XIII,  Appendix  C. 

The  confinement  of  these  patterns  was  achieved  by 
surrounding  the  pattern  of  interest  with  similar  patterns 
so  that  equipotential  boundaries  were  established  on  all 
sides  of  the  confined  pattern.  Figures  XX  aid  XXI  repre¬ 
sent  the  effect  of  such  a  confinement  on  oil  recovery  and 
areal  sweep  efficiency  for  a  five  spot  pattern.  The  data 
required  to  plot  these  figures  appear  in  Tables  XIV 
and  XV  of  Appendix  C.  It  may  be  observed  in  Figure  XX 
that  the  oil  recovery  for  the  confined  five  spot  pattern 
stabilizes  at  a  cumulative  throughput  of  4 . 6network  pore 
volumes,  whereas  it  shows  an  increasing  trend  in  the  case 
of  an  isolated  five  spot  pattern.  This  may  be  explained 
on  the  basis  that  surrounding  a  pattern  with  similar 


PATTERN  GROWTH  IN  CASE  OF  ISOLATED  NORMAL  FIVE  SPOT  PATTERN 

RUN  NO.  -  7 


Injection  Rate:  345  ccs/hr/well 


T:  1 .38  mins 


T:2.09  mins 


T :  4 . 1 4  mins 


T:6.29  mins 


T:  refers  to  the  total  time  elapsed  after  start 


FIGURE  XV 


PATTERN  GROWTH  IN  CASE  OF  ISOLATED  NORMAL  FIVE  SPOT 

RUN  NO.  -  7 


Injection  Rate:  345  ccs/hr/well 


T:  18.21  mins  T:23.87  mins 


T:  refers  to  the  total  time  elapsed  after  start 


FIGURE  XVI 


PATTERN  GROWTH  IN  CASE  OF  CONFINED  FIVE  SPOT 

RUN  NO.  -  11 


Injection  Rate:  345  ccs/hr/well 


T:0.97mins  T:1.42mins 


T:  refers  to  the  total  time  elapsed  after  start 


T :0.60  mins 


FIGURE  XVII 


PATTERN  GROWTH  IN  CASE  OF  CONFINED  FIVE  SPOT 

RUN  NO.  -  n 


Injection  Rate:  345  ccs/hr/well 


T:  1 .95  mins 


T:2.47  mins 
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T :3.00  mins 


T:3.  73  mins 


T:  refers  to  the  total  time  elapsed  after  start 


FIGURE  XVIII 


PATTERN  GROWTH  IN  CASE  OF  CONFINED  FIVE  SPOT 

RUN  NO.  -  11 


Injection  Rate: 


T:7.00  mins 


345  ccs/hr/we 1 1 


T:4.70  mins 


T:5.68  mins  T:6.53  mins 


T:  refers  to  the  total  time  elapsed  after  start 


FIGURE  XIX 
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patterns  on  all  sides  results  in  the  establishment  of 
equipotential  boundaries  across  which  no  cross-flow  occurs. 
Thus,  the  area  contacted  by  the  injected  water  is  limited 
to  these  boundaries.  These  boundaries,  however,  do  not 
exist  in  the  case  of  an  isolated  pattern. 

A  total  of  1.55  hydrocarbon  pore  volumes  of  oil 
was  recovered  form  the  confined  five  spot  pattern  which 
indicates  that  in  spite  of  completely  surrounding  the  five 
spot  pattern  with  similar  patterns  there  was  an  oversweep. 

On  the  basis  of  the  total  quantity  of  oil  recovered  and 
irreducible  oil  saturation  at  the  end  of  the  experimental 
run,  the  oversweep  was  found  to  be  71.6%.  These  calculations 
are  reported  in  Appendix  D. 

Figure  XXI  presents  areal  sweep  efficiency  as  a 
function  of  cumulative  throughput.  The  areal  sweep  effic¬ 
iency,  in  this  figure,  was  calculated  on  the  basis  of  the 
inside  pattern  area  and  the  outside  pattern  area.  The 
inside  patter  area  (Area  1)  here  implies  the  actual  area 
of  the  pattern  under  study,  whereas  the  outside  pattern  area 
(Area  II)  refers  to  the  area  of  the  white  portion  visible 
in  the  photograph  (Figure  XIX) . 

Figure  XXII  shows  the  effect  of  pattern  confinement 
on  the  relationship  between  instantaneous  water-oil  ratio 
and  throughput.  It  may  be  noted  that  water-oil  ratio 
steadily  increases  as  more  and  more  water  is  injected. 
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It  is  expected  that  after  a  certain  cumulative  water  injec¬ 
tion,  the  water-oil  ratio  should  rise  sharply  and  no  more 
oil  should  be  produced  with  further  water  injection.  This 
behaviour  is  truly  exhibited  by  the  confined  five  spot  and 
confined  direct  line  drive  patterns.  However,  in  the  case 
of  the  isolated  five  spot  pattern,  the  increase  in  water  - 
oil  ratio  is  much  slower  as  compared  to  the  confined  case. 

As  explained  earlier,  this  is  believed  to  be  due  to  the 
oil  production  from  more  and  more  freshly  contacted  areas 
which  compensated  for  less  and  less  subordinate  production 
from  the  swept  area. 

EFFECT  OF  PATTERN  CONFIGURATION 

Figure  XXIII  presents  the  results  of  experimental 
runs  6,  11  and  12  in  a  concise  form  and  shows  the  effect 
of  pattern  configuration  on  oil  recovery.  The  obvious 
conclusion  that  oil  recovery  is  higher  for  an  isolated 
five  spot  pattern  as  compared  to  the  confined  five  spot 
and  confined  direct  line  drive  patterns,  is  misleading  be¬ 
cause  there  is  an  oversweep  in  the  case  of  the  isolated  five 
spot  and  confined  five  spot  patterns.  This  figure,  however, 
does  bring  out  the  fact  that  as  more  and  more  rows  of  simi¬ 
lar  patterns  are  used  to  surround  the  pattern  of  interest, 
there  is  a  tendency  towards  achieving  complete  confinement. 
In  the  case  of  the  confined  five  spot,  the  central  five  spot 
was  confined  with  one  row  of  five  spots  around  its  boundary 


FIGURE  XXIII.  EFFECT  OF  WELL  PATTERN  ON  OIL  RECOVERY 
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and  as  stated  earlier  there  was  still  some  amount  of  over¬ 
sweep.  However,  in  the  case  of  confined  direct  line  drive, 
two  rows  of  similar  patterns  were  used  and  the  improvement 
in  confinement  may  be  noticed  in  Figure  XXIII. 
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CONCLUSIONS 

The  results  of  the  present  investigation  led  to 
the  following  conclusions. 

1)  There  was  a  definite  expansion  in  the  model  used  by 
Serra  and  Bhatia. 

2)  The  expansion  problem  was  partially  solved  in  the 
new  model. 

3)  The  porous  medium  used  in  the  present  study 'was 
reasonably  homogeneous. 

4)  Reproducible  results  were  obtained. 

5)  Back  pressure  was  found  to  have  no  effect  on  oil 
recovery  and  areal  sweep  efficiency. 

6)  Surrounding  a  five  spot  pattern  completely  with 
similar  patterns  resulted  in  an  oversweep  to  the 
extent  of  71.6%. 

7)  Complete  confinement  of  a  pattern  requires  more  than 
one  row  of  similar  patterns  surrounding  the  pattern 
under  investigation. 
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RECOMMENDATIONS 

1)  The  well  design  should  be  altered  so  that  a  constant 
thickness  is  maintained  between  the  two  lucite  sheets 

2)  Use  of  sintered  metal  on  a  well  stem,  in  the  packed 
region,  would  provide  a  more  uniform  flow  of  fluids 
through  the  wells  thereby  avoiding  the  cusping  (Fig¬ 
ure  XIX)  due  to  directional  orientation  of  the  per¬ 
forations  . 

3)  Use  of  glass  beads  and  lucite  sheets  in  the  present 
model  resulted  in  non-uniform  wettability  characteris 
tics  of  the  medium  because  the  beads  were  water-wet 
and  the  interfaces  between  the  glass  beads  and  lucite 
sheets  were  oil-wet.  This  problem  can  be  avoided  by 
using  glass  sheets,  if  a  water-wet  medium  is  required 
or  lucite  beads  if  it  is  desired  to  have  an  oil-wet 
medium. 


4)  Use  of  consolidated  porous  media  should  be  attempted. 
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Average  Thickness  Correlation 
Calculation  of  Porosity 
Calculation  of  Absolute  Permeability 
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AVERAGE  THICKNESS  CORRELATION 

Prior  to  packing  the  model  with  glass  beads,  a 
variation  in  the  gap  between  the  confining  lucite  sheets 
was  noticed.  It  was  observed  that  the  gap  was  equal  to  the 
thickness  cf  the  lucite  spacer  ring  at  the  periphery  of  the 
model  and  decreased  towards  the  center.  This  variation  was 
attributed  to  the  following  reasons. 

1.  small  irregularities  in  the  inside  surfaces  of  the 
lucite  sheets 

2.  sagging  in  the  model  due  to  the  weight  of  the  lucite 
sheets . 

In  order  to  measure  this  variation,  a  travelling 
microscope  was  used  and  the  gap  between  the  lucite  sheets 
determined  by  focussing  the  microscope  on  the  bottom  sur¬ 
face  of  the  top  plate  and  then  on  the  top  surface  of  the 
bottom  plate.  This  measurement  was  carried  out  along  four 
diameters  of  the  model  and  the  results  averaged  to  get  a 
representative  figure  for  the  variation  in  thickness.  The 
data  obtained  were  as  follows. 


RADIUS  (From  Center) 


AVG.  THICKNESS 


Inches 


Cms  < 


22.6400 

18.8660 

15.0930 


0 .6158 
0.5565 
0.5108 


11.3200 


0.4830  con  t  d . 
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RADIUS  (From  Center) 
Inches 

7.5466 
3.7730 
0.0000 
3.7730 
7.5466 
11.3200 
15.0930 
18 . 8660 
22 . 6400 


AVG.  THICKNESS 
Cms . 

0.4486 

0 .4233 
0 .4222 
0.4210 
0.4494 
0.4600 
0.5104 
0 . 5425 
0.6133 


CALCULATION  OF  POROSITY 

1 .  Measurement  of  Matrix  Volume 

Total  weight  of  glass  beads  added  =  9122.39  gms . 

Specific  gravity  of  glass  beads  =  2.515 

.* .  Matrix  Volume  of  the  pack  =  9122.39 

2 . 515 

=  3630.42  c . cs . 


Bulk  Volume  of  the  model  =  5711.0  c.cs. 
(measured  by  filling  the 
model  with  water) 


Porosity,  % 


Bulk  Volume  -  Matrix  Volume 
Bulk  Volume 


x  100 


=  5711.0  -  3630.42 
5711.0 


x  100 


36 .43% 
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2 .  Measurement  of  Pore  Volume 

The  pore  volume  of  the  model  was  determined  by 

timing  the  process  of  saturating  the  model  completely  with 

water.  Knowing  the  injection  rate  and  the  total  time,  pore 

volume  could  be  easily  calculated. 

Time  required  to  completely  saturate  the  medium  = 

10.40  hrs . 

Injection  rate  =  200.0  c.cs/hr. 

.  .  Pore  Volume  =  200  x  10.40  c.cs. 

=  2080.00  c.cs. 

2080.0 

.  .  Porosity  =  57Yi~~0  x  100 

=  36.425% 

CALCULATION  OF  ABSOLUTE  PERMEABILITY 

The  porous  medium  was  completely  saturated  with 
water.  A  five  spot  pattern  was  then  used  to  measure  abso¬ 
lute  permeability.  The  calculations  were  performed  using 
Muskat ' s  steady  state  equation  for  a  five-spot  pattern. 

Muskat * s  Equation: 

K  =  Q  (li-L  (d/rwj  -  0-619) 

W  0.003541  AP  h 

where  , 

K  =  effective  permeability  to  water,  mds . 
w 

Q  =  rate  of  injection,  bbls/day 
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u  =  viscosity  of  water,  cps . 
w 

d  =  distance  between  injection  well  and  production  well, 
f  t . 

r  =  well-bore  radius,  ft. 
w  ' 

AP  =  Pressure  differential,  psi. 
h  =  formation  thickness,  ft. 

In  the  present  model , 

Q  =  532.0  c.cs/hr/well  =  2128.0  c.cs/hr 
=  0.3205  bbls/day 

AP  =  5.10  cms .  of  mercury 
=  0.987  psi . 

h  =  Average  thickness  for  11.3125  inch  side  five  spot 
=  0.4517  cms. 

=  0.01482  ft. 

y  =  0.970  cps. 
w  ^ 

r  =  0.0485  inch 
w 

=  0.00404  ft. 

d  =  8  inches 
=  0.6670  ft. 

d/r  =  0.6670/0.00404  =  165.20 

In  (d/r  )  =  In  (165.20)  =  5.230 
w 

*  K_  =  0.3205  x  0 .970  x  (5  ■  230  -  0 . 619) 

*  '  w  0.003541  x  0.01482  x  0 .987 


=  27,700  mds . 
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Calculation  of  critical  injection  rate. 
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BASIC  DATA  REQUIRED  TO  USE  THE  METHOD 

The  following  data  are  required  to  predict  the 
performance  of  any  pattern  using  the  modified  scheme: 

1.  Relative  permeability  versus  saturation  curves. 

2.  Pore  volumes  and  shape  factors  for  each  channel. 

3.  Initial  water  saturation  (connate  or  irreducible) 
and  irreducible  oil  saturation. 

4.  Maximum  water-oil  ratio  at  which  cut-off  is 
desired . 

5.  Pressure  drop  between  the  injection  and  the 
production  well  for  a  constant  pressure  drop  case 

6.  Injection  rate  for  a  constant  injection  rate  case 

7.  Time  versus  injection  rate  history  where  both 
injection  rate  and  pressure  drop  vary. 

8.  Absolute  permeability  of  the  medium. 

9.  Oil  and  water  viscosities. 
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3.5  x  10  3  x  32.14  x  (27700  x  0 . 3 6 4 3 )  0  *  5bbls/day/f t 

0.970 


3.5  x  10  3  x  32.14  x  100.45 

U.970 


bbls/day/f t 


=  11.649  bbls/day/ft. 


It  may  be  observed  that,  the  critical  injection 
rate,  as  expressed  here,  depends  upon  the  thickness  of 
the  porous  medium.  Thus  the  higher  the  value  used  for 
thickness,  the  higher  is  the  resulting  injection  rate.  A 
maximum  average  thickness  of  0.5000  cms .  was,  therefore , 
selected  for  calculation  purposes. 


Maximum  average  thickness  =  0.5000  cms. 

=  0.0164  ft. 

thus  q  =  11.649  x  0.0164 

=  0.1910  bbls/day 
=  1268.17  c.cs/hr. 

This  implies  that  if  the  injection  rate  is  kept 
above  1268.17  c.cs/hr.,  the  resulting  floods  will  be 
stabilized  and  will  no  longer  be  rate  -  sensitive.  The 
minimum  injection  rate  used  in  this  study  was  1380.0  c.cs/ 
hr. 


APPENDIX  -  C 


Measurement  of  Expansion  in  the  old  model 
Measurement  of  Expansion  in  the  new  model 
Results  of  Radial  flow  experiment. 
Experimental  data  on  isolated  five  spot 
system  at  five  different  back  pressures 
Calculation  of  Oil  Recovery  and  Areal 
Sweep  Efficiency  for  isolated  five  spot 
pattern . 

Calculation  of  Oil  Recovery  and  Areal 
Sweep  Efficiency  for  confined  five  spot 
pattern . 

Calculation  of  Oil  Recovery  for  confined 
direct  line  drive  pattern. 
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NOMENCLATURE  USED  IN  THE  TABLES 

N 

=  Initial  oil  in  place  in  the  model,  c.cs. 

L 

=  Length  of  one  side  of  the  well  -  pattern,  cms 

d 

=  Distance  between  the  injection  and  the  pro¬ 
duction  wells,  cms. 

r 

w 

=  Well  bore  radius,  cms. 

h 

avg 

=  Average  thickness  of  the  packing  in  the  well 
pattern 

q 

=  Injection  rate,  c.cs/hr. 

0 

=  Porosity,  %. 

s°r 

=  Residual  oil  saturation 

So  . 

1 

=  Initial  oil  saturation 

Sw 

c 

=  Connate  water  saturation 

N. W.P .V. 

=  Net-work  pore  volume 

H.C.P.V. 

=  Hydrocarbon  pore  volume 

W .  0  .  R . 

=  Water  -  oil  ratio 

TABLE  III 

MEASUREMENT  OF  EXPANSION  IN  THE  OLD  MODEL 
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TABLE  IV 

MEASUREMENT  OF  EXPANSION  IN  THE  NEW  MODEL 


Obs.  No.  PRESSURE  VOL.  OF  OIL  INJECTED 
cms .  of  Water  c . cs . 

1  18.50  0.856 

2  60.70  2.551 

3  77.90  3.427 

4  99.70  4.270 

5  174.93  7.726 

6  201.33  9.535 


NET  INCREASE  IN  THICKNESS  OF  THE  MODEL 

Referring  to  figure  VIII,  a  pressure  of  70  cms. 
of  water  causes  an  increase  of  3.05  c.cs.  in  the  bulk 
volume  of  the  model. 

Internal  diameter  of  the  model  =  45.281  inches 
.  .  ,  3.05 

.  .  Net  increase  m  thickness  =  - — - - p 

tt/4  x  (45.281  x  2.54)^ 


=  0.0003  cms. 
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TABLE  V 

RESULTS  -  RADIAL  FLOW  EXPERIMENT 


BACK  PRESSURE 

Obs .  No. 

FLOW  RATE 

AVG.  PRESSURE  DROP 

Cms .  of  Water 

c.cs/hr . 

cms.  of  Hg. 

1 

161.550 

0.370 

9.635 

2 

341.100 

0.770 

3 

436.620 

0.975 

4 

520.250 

1.180 

5 

606.880 

1.450 

1 

162.884 

0.350 

2 

327.350 

0.700 

21.240 

3 

428 . 373 

0.980 

4 

520.570 

1.230 

5 

604 . 670 

1.450 

1 

170.310 

0.350 

2 

334.910 

0.716 

33.740 

3 

433.120 

0.967 

4 

517.850 

1.175 

5 

607.820 

1.450 

1 

168 .400 

0.350 

2 

338 .930 

0 .730 

54.545 

3 

434 . 360 

0.950 

4 

518 . 670 

1.175 

5 

605 . 500 

1.450 

C-5 


RUN  NO .  5 

ISOLATED  NORMAL  FIVE  SPOT 


N 

=  1774.60  c.cs. 

So  = 
r 

0.0484 

L 

=  14.40  cms . 

So  .  = 

l 

0.8530 

d/r 

w 

=  82.60 

Sw  = 
c 

0.1470 

avg . 

=  0.4250  cms . 

q  = 

345  c.cs/hr/well 

0 

=  36.43% 

N.W.P.V.  = 

32 . 3116  c.cs. 

Avg.  Pattern  area  =  0.0086  units  H.C.P.V.  =  32.3116  x 

0.8530 

=  27.562  c.cs. 

Back  Pressure  =  3.335  cms .  of  water 


PRODUCTION  HISTORY 


Obs .  No . 

Water 

Prod . 

Oil 

Prod . 

Total 

Prod . 

Inst . 

Cum . 

Inst . 

Cum. 

Inst . 

Cum. 

1 

0.0 

0.0 

11.10 

11.10 

11.10 

11.10 

2 

0 . 0 

0.0 

6.50 

17.60 

6 . 50 

17.60 

3 

0 . 0 

0.0 

10.40 

28.00 

10.40 

28 . 00 

4 

0 . 0 

0.0 

10 . 30 

38 . 30 

10.30 

38.30 

5 

0.0 

0.0 

9.60 

47 .90 

9 . 60 

47.90 

6 

25.0 

25 . 0 

30.00 

77 . 90 

55.00 

102.90 

7 

55.0 

80.0 

24 . 50 

102.40 

79.50 

182.40 

8 

59.0 

139 . 0 

17.50 

119 .90 

76 . 50 

258.90 

9 

82.0 

221.0 

27 . 00 

146.90 

109.00 

367.90 

10 

81.0 

302.0 

23.00 

169.90 

104.00 

471.90 

11 

114.0 

416 . 0 

33.00 

202.90 

147.00 

618.90 

12 

139 . 0 

555.0 

38 .00 

240.90 

177.00 

795.90 

* 

28.0 

583.0 

7.00 

247 . 90 

35.00 

830.90 

*  Production  after  injection  ceased. 
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RUN  NO .  6 


ISOLATED  NORMAL  FIVE  SPOT 


N 

=  1781.40  c . cs . 

So 

r 

= 

0.0503 

L 

=  14.40  cms . 

So  . 

_ 

0.8562 

i 

d/r 

w 

=  82.60 

SWC 

= 

0.1483 

h 

avg . 

=  0.4250  cms . 

q 

= 

345  c.cs/hr/well 

0 

=  36.43% 

N.W.P.V. 

= 

32 . 3116  c . cs . 

H.C.P.V. 

= 

32.3116  x  0.8562 

= 

27.6652  c .cs . 

Back  : 

Pressure  =  13.335  cms. 

of  water 

Obs  . 

PRODUCTION  HISTORY 

No.  Water  Prod.  Oil  Prod. 

Total 

Prod . 

Inst 

.  Cum. 

Inst . 

Cum. 

Inst . 

Cum. 

1 

0 . 00 

0 . 00 

9.90 

9.90 

9.90 

9.90 

2 

0 . 00 

0.00 

7.00 

16.90 

7 . 00 

16.90 

3 

0.00 

0.00 

11.90 

28.80 

11.90 

28.80 

4 

0.00 

0.00 

14 . 00 

42.80 

14 . 00 

42.80 

5 

0 . 20 

0.20 

10 .00 

52.80 

10.20 

53.00 

6 

26.00 

26.20 

25.50 

78 . 30 

51.50 

104.50 

7 

52.00 

78.20 

24.50 

102.80 

76 . 50 

181.00 

8 

56 . 00 

134.20 

20 . 00 

122.80 

76.00 

257.00 

9 

79 . 00 

213.20 

24 . 00 

146.80 

103.00 

360.00 

10 

90.00 

303 . 20 

26 . 00 

172.80 

116 . 00 

476.00 

11 

108 . 00 

411.20 

29.00 

201.80 

137.00 

613.00 

12 

138.00 

549.20 

37.00 

238 . 80 

175.00 

788.00 

* 

37 . 00 

586 . 20 

8 . 00 

246.80 

45.00 

833.00 

*  Production  after  injection  ceased. 
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RUN  NO .  7 

ISOLATED  NORMAL  FIVE  SPOT 


N 

=  1831.70  c.cs. 

S°r 

= 

0.0716 

L 

=  14.40  cms . 

So. 

1 

= 

0.8804 

d/r 

w 

=  82.60 

Sw 

c 

= 

0.1196 

h 

avg . 

=  0.4250  cms . 

q 

= 

345  c.cs/hr/well 

0 

=  36.43% 

N.W.P.V. 

= 

32.3116  c.cs. 

Avg .  : 

Pattern  area  =  0.0075 

H.C.P. V. 

32.3116  x 

0.8804 

as 

28.450  c.cs. 

Back  Pressure  =  24.255  cms .  of  water 


PRODUCTION  HISTORY 


Obs 

.  No.  Water 

'  Prod. 

Oil  Prod. 

Total 

Prod . 

Inst . 

Cum. 

Inst . 

Cum. 

Inst . 

Cum. 

1 

0.0 

0.0 

1 . 50 

1.50 

1.50 

1.50 

2 

0.0 

0.0 

15.00 

16 . 50 

15.00 

16 . 50 

3 

0 . 0 

0.0 

15.20 

31.70 

15.20 

31.70 

4 

0 . 50 

0.50 

15.90 

47 .60 

16 . 40 

48.10 

5 

24 . 60 

25.10 

22.50 

70 . 10 

47.10 

95.20 

6 

32.00 

57.10 

17.50 

87.60 

49.50 

144 . 70 

7 

52 . 50 

109 . 60 

22.50 

110.10 

75.00 

219 . 70 

8 

75.50 

185.10 

22 . 50 

132.60 

98 . 00 

317 .70 

9 

79 . 00 

264 . 10 

22.20 

154.80 

101.20 

418 . 90 

10 

102.00 

366.10 

28 . 00 

182.80 

130.00 

548.90 

11 

163.00 

529.10 

42 . 00 

224 .80 

205.00 

753.90 

* 

52 .00 

581.10 

12.00 

236 . 80 

64 . 00 

817.90 

* 

Production 

after 

injection 
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RUN  NO .  8 

ISOLATED  NORMAL  FIVE  SPOT 


N  = 
L  = 


d/r  = 
/  w 


h 


avg 


1812.90  c . cs 


14.40  cms . 


82.60 


0.4250  cms 


0  =  36.43% 


So  = 


So .  = 


Sw  _ 


c 

q 

N.W.P.V. 


Avg.  Pattern  area  =  0.0077  uts .  H.C.P.V.  = 


0.0560 

0.8713 


=  0.1287 

345  c.c s/hr/we 11 
32 . 3116  c . cs . 


32.3116  x 
0.8713 


=  28.153  c . cs . 


Back  Pressure  =  45.455  cms.  of  water 


PRODUCTION  HISTORY 


Obs  . 

No.  Water 

Prod . 

Oil 

Prod . 

Total 

Prod . 

Inst . 

Cum. 

Ins  t . 

Cum. 

Inst . 

Cum. 

1 

0.00 

0.00 

9 . 60 

9 .60 

9 .60 

9.60 

2 

0.00 

0.00 

11.50 

21.10 

11.50 

21.10 

3 

0.00 

0.00 

12.40 

33.50 

12.40 

33.50 

4 

0.00 

0.00 

10.50 

44 . 00 

10.50 

44.00 

5 

23.00 

23 . 00 

30.00 

74 . 00 

53.00 

97.00 

6 

51.00 

74.00 

27.00 

101.00 

78.00 

175.00 

7 

56 . 50 

130 . 50 

23.50 

124.50 

80.00 

255.00 

8 

78.00 

208 . 50 

30.00 

154.50 

108.00 

363.00 

9 

82 . 00 

290 . 50 

27.00 

181.50 

109.00 

472.00 

10 

109 . 00 

399 . 50 

34.00 

215.50 

143.00 

615.00 

11 

142.00 

541.50 

39 . 00 

254.50 

181 . 00 

796.00 

* 

38.00 

579 . 50 

13.00 

267 . 50 

51.00 

847.00 

* 


Production  after  flection  ceased 
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RUN  NO .  9 

ISOLATED  NORMAL  FIVE  SPOT 


N 

=  1747.70  c.cs. 

I- 

So 

r 

= 

0.0392 

L 

=  14-40  cms . 

So  . 

l 

= 

0 . 8400 

d/r 

w 

=  82.60 

Sw 

c 

= 

0.1600 

h 

avg . 

=  0.4250  cms. 

q 

= 

345  c.cs/hr/well 

0 

=  36.43% 

N.W.P.V. 

- 

32.3116  c.cs. 

Avg . 

Pattern  area  =  0.0085 

units  H.C.P.V. 

= 

32.3116  x 

0 . 8400 

= 

27.142  c.cs. 

Back 

Pressure  =  69.275  cms 

.  of  water 

Obs  . 

PRODUCTION  HISTORY 

No.  Water  Prod.  Oil  Prod. 

Total  Prod. 

Inst . 

Cum. 

Inst . 

Cum. 

Inst . 

Cum. 

1 

0.00 

0.00 

6 . 50 

6.50 

6.50 

6.50 

2 

0 .00 

0.00 

11 .80 

18 . 30 

11.80 

18 . 30 

3 

0.00 

0 . 00 

13.50 

31.80 

13.50 

31.80 

4 

0 . 00 

0.00 

13 . 30 

45.10 

13 . 30 

45.10 

5 

21.00 

21.00 

28.00 

73.10 

49.00 

94.10 

6 

51.00 

72 .00 

27.50 

100.60 

78.50 

172 .60 

7 

76.00 

148 . 00 

30.50 

131.10 

106.50 

279.10 

8 

78 . 00 

226.00 

29 . 00 

160.10 

107 . 00 

386.10 

9 

104.00 

330.00 

32.00 

192.10 

136.00 

522.10 

10 

111.00 

441.00 

35.00 

227.10 

146.00 

668.10 

11 

136.00 

577 . 00 

39 .00 

266 . 10 

175.00 

843.10 

* 

27 . 00 

604 .00 

8 . 00 

274.10 

35.00 

878.10 

*  Production  after  injection  ceased 
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RUN  NO.  10 

ISOLATED  NORMAL  FIVE  SPOT 


N  = 

1835.80  c.cs. 

So 

r 

= 

0.0678 

L  = 

14.40  cms. 

So . 

_ 

0.8828 

l 

d/r  = 
w 

82.60 

Sw 

c 

= 

0.1172 

avg .  ~ 

0.4250  cms. 

q 

= 

345  c.cs/hr/well 

0  = 

36.43% 

N.W.P.V. 

= 

32 . 3116  c.cs. 

H.C.P.V. 

= 

32.3116  x 

0.8828 

=  28.50  c . cs . 

Back  Pressure  =  48.100  cms .  of  water 


PRODUCTION  HISTORY 


Obs . 

No.  Water  Prod. 

Oil 

Prod . 

Total 

Prod . 

Inst . 

Cum. 

Inst . 

Cum . 

Inst . 

Cum. 

1 

0.0 

0 . 0 

0.60 

0 . 60 

0.60 

0.60 

2 

0.0 

0.0 

15.80 

16.40 

15.80 

16.40 

3 

0.0 

0 . 0 

15.70 

32.10 

15.70 

32.10 

4 

0.0 

0.0 

13.50 

45.60 

13.50 

45.60 

5 

22.0 

22.0 

30.00 

75.60 

52.00 

97 .60 

6 

52 . 0 

74.0 

26.00 

101.60 

78.00 

175.60 

7 

57.0 

131.0 

24.00 

125.60 

81.00 

256.60 

8 

80.0 

211.0 

28 . 00 

153.60 

108.00 

364.60 

9 

80.0 

291.0 

23.00 

176.60 

103.00 

467.60 

10 

103.0 

394 . 0 

32.00 

208.60 

135.00 

602.60 

11 

150.0 

544.0 

45.00 

253.60 

195.00 

797 . 60 

* 

23.0 

567.0 

8 .00 

261.60 

31.00 

828.60 

* 


Production  after  injection  ceased. 


• 

CALCULATION  OF  OIL  RECOVERY 
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RUN  NO.  11 

CONFINED  NORMAL  FIVE  SPOT 


N 

=  1848.40  c.cs. 

So 

r 

= 

0.0704 

L 

=  14.40  cms. 

So . 

l 

= 

0.8884 

d/r 

w 

=  82.60 

Swc 

= 

0.1116 

h 

avg . 

=  0.4250  cms. 

q 

= 

345  c.cs/hr/well 

0 

=  36.43% 

N.W.P.V. 

= 

32 . 3116  c.cs. 

Avg .  ! 

Pattern  area  -  I  =  0.02115 

II  =  0.0838 

H ,  C  .  P  .  V . 

units 

= 

32.3116  x 

0.8884 

=  28.706  c.cs. 

Back  Pressure  =  13.005  cms .  of  water 


PRODUCTION  HISTORY 


Obs  . 

No.  Water  Prod. 

Oil  Prod. 

Total  Prod. 

Inst . 

Cum. 

Inst . 

Cum. 

Inst . 

Cum. 

1 

0.00 

0.00 

5.90 

5.90 

5.90 

5.90 

2 

0.00 

0.00 

7.80 

13.70 

7.80 

13.70 

3 

0.00 

0.00 

8.60 

22.30 

8.60 

22 . 30 

4 

3.20 

3.20 

7.10 

29 . 40 

10 . 30 

32.60 

5 

7 .80 

11.00 

4.40 

33.80 

12.20 

44.80 

6 

9.00 

20.00 

3.10 

36.90 

12.10 

56.90 

7 

10 . 50 

30.50 

1.50 

38.40 

12.00 

68.90 

8 

14.70 

45.20 

2.20 

40.60 

16.90 

85.80 

9 

19.80 

65.00 

2.60 

43.20 

22.40 

108.20 

10 

21.50 

86.50 

1.00 

44.20 

22.50 

130.70 

11 

18 .60 

105.10 

0.80 

45.00 

19.40 

150.10 

12 

11.10 

165.10 

0.00 

45.00 

11.00 

161.10 

TABLE  XII 
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RUN  NO.  12 

CONFINED  DIRECT  LINE  DRIVE 


N  = 

1692.0 

c  .  cs  . 

So 

r 

= 

0.0596 

L  = 

14  .40 

cms . 

So . 

_ 

0.8133 

l 

d/r  = 
w 

82 . 60 

Sw 

c 

= 

0.1867 

h  = 

avg . 

0.4250 

cms . 

q 

= 

345  c.cs/hr/well 

0  = 

36 .43% 

N. W.P.V. 

= 

54.672  c . cs . 

H.C.P.V. 

= 

54.672  x 

0.8133 

= 

44.480  c . cs . 

Back  Pressure  = 

13.005  cms. 

of  water 

PRODUCTION  HISTORY 

Obs  . 

No.  Water 

Prod . 

Oil 

Prod . 

Total  Prod. 

Inst . 

Cum. 

Inst . 

Cum. 

Inst . 

Cum. 

1 

0.000 

0 . 000 

6.000 

6.000 

6.00 

6.00 

2 

0.000 

0.000 

8.000 

14 . 000 

8.00 

14.00 

3 

4.932 

4.932 

5.568 

19.568 

10.50 

24.50 

4 

16 . 266 

21.198 

6.934 

26.502 

23.20 

47.70 

5 

20 .667 

41.865 

3.433 

29.935 

24.10 

71.80 

6 

23.131 

64.996 

2.869 

32.804 

26.00 

97 .80 

7 

32.792 

97.788 

2.408 

35.212 

35.20 

133.00 

8 

23.626 

121.414 

1.274 

36.486 

24.90 

157.90 

9 

24.099 

145.513 

1.001 

37.487 

25.10 

183.00 

10 

19 . 359 

164.872 

0.541 

38.028 

19.90 

202.90 

11 

32.878 

197.750 

0.622 

38.650 

33.50 

236.40 

12 

42.450 

240.200 

0.550 

39.200 

43.00 

279 . 40 
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APPENDIX 


D 


Calculation  of  displacement  efficiency  for: 

a.  Isolated  five  spot 

b.  Confined  five  spot 

Calculation  of  oversweep  in  confined  five  spot. 
Calculation  of  oil  recovery  from  areal  sweep 
efficiency  versus  cumulative  throughput  curve. 
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CALCULATION  OF  DISPLACEMENT  EFFICIENCY  FOR  ISOLATED  NORMAL 
FIVE  SPOT 


Displacement  Efficiency, 

E  _  Oil  Recovered  at  any  stage _ 

d  Area  contacted  by  flooding  water  at 
this  stage  x  Average  thickness  x 
Porosity  x  Initial  Oil  Saturation 


In  the  case  of  Run  No.  9, 

Porosity  =  0.3643 

Initial  Oil  saturation  -  0.8400 

Referring  to  Table  X  it  is  observed  that  at  a 
total  injection  of  5.3420  network  pore  volumes  of  water. 
Recovery  of  oil  =  3.7060  H.C.P.V. 

=  3.7060  x  27.142 
=  100.60  c . cs . 

Areal  sweep  efficiency  =  6.8 

.  2 
.  .Area  contacted  by  flooding  water  =  6.8  x  (14.4) 

=  1410.10  sq.  cms . 

Average  thickness  for  this  area,  from  Appendix  A 

=  0.4533  cms. 


100.60 _ _ _ 

1410.10  x  0.45333  x  0.3643  x  0.8400 


x  100 


51.50% 
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CALCULATION  OF  DISPLACEMENT  EFFICIENCY  FOR  CONFINED  FIVE 
SPOT 

In  the  case  of  confined  five  spot  the  areal 
sweep  corresponding  to  the  oil  recovered  is  undefined  be¬ 
cause  of  the  manner  in  which  pattern  confinement  was  achieved. 
As  such  the  following  definition  of  displacement  efficiency 
has  been  used. 


Displacement  Efficiency, 

„  Total  oil  recovered 

E^  =  - 

Oil  in  place 

, ,  Sw  So  x 
(1  -  c  -  r) 


(i  -  Swi> 


In  the  case  of  Run  No.  11 

Sor  =  0.0704 
Sw .  =  0.1116 

l 

p  =  (1  ~  0.1116  -  0.0704)  100 

.  .  nd 

(1  -  0.1116) 


0 .8180 
0 .8884 

-  92.1% 


100 


It  may  be  noted  that  this  value  of  displacement 
efficiency  is  applicable  only  when  complete  flooding  of  the 
pattern  has  occurred. 
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CALCULATION  OF  OVERSWEEP  IN  CONFINED  FIVE  SPOT 

In  the  case  of  Run  No.  11, 

Porosity  =  0.3643 

Residual  oil  saturation  =  0.0704 

Connate  water  saturation  =  0.1116 

H.C.P.V.  =  28.706 

Avg s  thickness  =  0.4250  cms . 

Referring  to  Table  XII  it  is  observed  that  oil  re¬ 
covery  stabilizes  at  1.5650  H.C.P.V. 

Thus  total  oil  recovered  from  the  pattern 
=  1.5650  x  28.706 
=  45.0  ccs . 

Oil  recovered  can  also  be  determined  from  the  following 
equation : 

Oil  recovered  =  Area  swept  x  thickness  x  porosity  x 

(1  -  Residual  oil  saturation  - 
connate  water  saturation) 

45.0 

Area  swept  =  - - 

0.4250  x  0.3643  x  (1  -  0.0704  -  0.1116) 

=  355.31  sq.  cms. 

Pattern  Area  =  14.4  x  14.4 

=  207.0  sq.  cms. 

•  .  _ _ .  355.31  ,A 

.  .  Areal  sweep  efficiency  =  -  x  100 

207.0 
=  171.6% 

Thus  oversweep  =  171.6  -  100 

The  extra  area  swept  on  the  basis  of  this  areal  sweep 
efficiency  is  indicated  by  the  shaded  portion  on  Figure  XXIV. 


CONFINED  FIVE  SPOT 
REPRESENTATION  OF  EXTRA  AREA  SWEPT 


FIGURE  -  XXIV 
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CALCULATION  OF  OIL  RECOVERY  FROM'  AREAL  SWEEP  EFFICIENCY 
VERSUS  CUMULATIVE  THROUGHPUT  CURVE 

In  the  case  of  Run  No.  7  on  Isolated  Normal  Five  Spot, 
Porosity  =  0.3643 

Irreducible  oil  satruation  =  0.0716 
H.C.P.V.  =  28.45  c.cs. 

Referring  to  Table  VIII,  it  is  observed  that 
at  a  total  injection  of  2.95  network  pore  volumes  of  water, 
Oil  recovered  =  2.4650  H.C.P.V. 

Areal  sweep  efficiency  =  4.70 
Thus,  area  contacted  by  flooding 

water  =  4.70  x  (14. 4)2 
=  974.59  sq.  cms . 

Average  thickness  for  this  area,  from  Appendix  A 

=  0.4400  cms. 

Amount  of  oil  displaced  when  974.59  sq.  cms.  of  area  is 
contacted  by  flooding  water 

=  Area  contacted  x  avg .  thickness 
x  porosity  x  (1  -  residual  oil 
saturation)  x  displacement  ef¬ 
ficiency  . 

.  .  Oil  displaced  =  974.59  x  0.440  x  0.3643  x  (1  -  0.0716) 

x  0.5150 


74.6  c.cs. 
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Comparing  it  to  the  oil  recovered: 

Oil  recovered  =  2.4650  H.C.P.V. 

=  2.4650  x  28.45 
=  70.2  c ♦ cs . 

Thus,  it  may  be  observed  that,  within  the  experi¬ 
mental  errors,  the  agreement  between  oil  recovered  and 
oil  displaced  is  reasonably  close. 
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INTRODUCTION 

The  objective  of  this  investigation  was  to  check 
the  compatibility  of  laboratory  results  against  the  results 
obtained  from  theoretical  prediction,  using  the  proper¬ 
ties  of  the  model.  Various  waterflood  prediction  methods 
which  have  been  proposed  in  the  literature  from  time  to 
time  were  reviewed.  It  was  decided  to  use  the  method 
suggested  by  Higgins  and  Leighton  (1)  due  to  its  simpli¬ 
city  . 

An  examination  of  the  Higgins-Leighton  (2) . 
method  revealed  that  it  is  applicable  only  to  the  case 
where  a  constant  pressure  drop  is  maintained  between  the 
injection  and  the  production  well.  Since  the  experimen¬ 
tal  work  was  conducted  at  constant  injection  rate  condi¬ 
tions,  a  few  modifications  had  to  be  made  in  the  original 
scheme  before  it  could  be  applied  to  the  present  work. 
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LITERATURE  SURVEY 

Waterflood  prediction  methods  give  estimates 
of  the  following: 

-  ultimate  production  to  be  attained 

-  time  rate  of  oil  production  both  before  and 
after  breakthrough 

-  time  rate  of  water  production  after  breakthrough 

Depending  on  the  purpose  of  the  waterflood  pre¬ 
diction,  this  estimate  of  performance  may  range  from  a 
thorough  analysis  of  all  factors  known  to  affect  the  flood 
to  a  simple  rule  of  thumb  for  preliminary  evaluation . The 
first  step  in  the  development  of  a  prediction  method  is 
to  overcome  the  problems  caused  by  significant  reservoir 
heterogeneities.  It  was  recognized  very  early  in  the 
literature  that  the  success  or  failure  of  a  flood  depended 
upon  the  permeability  stratifications  existing  within  the 
reservoir.  Later  on  it  was  discovered  that  the  injection 
well  and  its  injectivity  characteristics  were  fundamental 
to  the  study  of  such  a  project.  The  classical  methods 
which  developed  from  these  considerations  are  those  pro¬ 
posed  by  Dykstra-Parsons  (3), Stiles  (4),  Yuster  -  Suder  - 
Calhoun  (5)  and  much  later  by  Pratts,  Mathews,  Jewett, 
Baker  (6 ) . 

While  these  methods  were  primarily  concerned 
with  permeability  and  injectivity  variations,  there 
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developed  at  the  same  time  an  appreciation  of  the  dis¬ 
placement  mechanism  which  was  responsible  for  the  overall 
floodability  of  a  reservoir.  The  two  areas  of  primary 
concern  were: 

-  a  definition  of  the  ma croscopic  variables  control¬ 
ling  the  areal  sweep  efficiency  during  and  after 
water  breakthrough  into  the  producing  well  and 

-  the  definition  of  microscopic  variables  control¬ 
ling  the  displacement  of  oil  by  water. 

Both  these  effects  are  influenced  by  the 
oil  viscosity  and  relative  permeability  characteristics 
of  the  porous  medium.  The  former  problem,  dealing  with 
the  definition  of  macroscopic  variables,  was  more  emphasized 
in  the  case  of  pattern  floods  where  the  proportion  of 
reservoir  volume  swept  to  breakthrough  was  a  variable. 
Significant  contributions  along  this  line  evolved  in  the 
prediction  methods  of  Muskat  (7),  Hurst  (8),  and  a  little 
later  Slobod-Caudle  (9),  Dyes  -  Caudle  -  Erickson  (10), 
and  Caudle-Witte  (11) . 

The  fundamental  theory  developed  by  Buckley  - 
Leverett  (12)  for  linear  flood  behavior  has  been  helpful 
in  analyzing  physical  problems  associated  with  the  mic¬ 
roscopic  displacement  mechanism.  This  theory,  however, 
was  not  put  into  a  convenient  form  for  pattern  flooding 
until  incoroporated  in  the  prediction  method  of  Craig, 
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Geffen  and  Morse  (13) . 

Since  then,  the  method  proposed  by  Higgins  and 
Leighton  (14)  has  also  used  this  approach  with  a  practical 
result.  This  method  includes  a  microscopic  formulation 
of  the  Buckley-Leverett  model  and  is  able  to  predict  the 
behavior  of  a  flood  developed  with  almost  any  type  of 
injection/production  well  arrangement. 

Each  published  prediction  method  has  been  more 
or  less  slanted  towards  a  definition  of  one  or  more  pri¬ 
mary  variables  known  to  affect  the  water  flood  performance. 
Pertinent  variables  which,  however,  must  be  considered  in 
any  comprehensive  prediction  method  are: 

-  Permeability  distribution 

-  Injectivity  or  Injection  rate 

-  Areal  Sweep  Efficiency 

-  Mobility  Ratio 

-  Displacement  Mechanism 

-  Non  -uniformities 

PERMEABILITY  DISTRIBUTION 

Early  prediction  methods  were  concerned  with 
permeability  changes  due  to  reservoir  stratification.  The 
availability  of  computers  today  has  virtually  eliminated 
this  problem  as  predictions  can  now  be  made  rapidly  on 
single  layers  and  superimposed  to  get  the  combined  strati¬ 
fied  reservoir  performance.  Higgins  and  Leighton  (15) 
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extended  this  approach  by  subdividing  a  single  homogeneous 
bed  into  individual  flow  channels  or  ducts.  Performance 
calculations  are  made  for  each  individual  channel,  summed 
to  get  the  layer  performance  and  superimposed  with  the 
performance  from  other  layers  to  get  the  overall  reservoir 
performance. 

INJECTIVITY  OR  INJECTION  RATE 

Time  can  be  introduced  into  the  performance 
prediction  only  through  a  knowledge  of  the  expected  injec¬ 
tion  rate  as  production  continues.  Since  roughly  80%  of 
the  total  pressure  drop  in  a  five  spot  pattern  occurs 
in  about  4%  of  the  volume  near  the  injection  well  and  4% 
near  the  producer,  injection  rates  calculated  from  dif¬ 
ferential  pressure  are  sensitive  to  viscous  effects  in 
the  oil  and  water  banks.  An  excellent  analysis  of  this 
problem  has  been  presented  by  Prats  et.al.  (16).  Perhaps 
the  largest  single  source  of  error  in  predicting  reservoir 
performance  is  the  need  for  complete  accountability  of 
injected  water.  Fractures,  casing  and  packer  leaks,  im¬ 
properly  abandoned  wells,  aquifer  injection,  directional 
permeability,  and  regional  pressure  gradients  pose  possible 
reasons  for  differences  between  the  actual  and  predicted 
injectivity  and  production  performance. 


6 


AREAL  SWEEP  EFFICIENCY 

Considerable  ingenuity  was  required  to  solve 
the  areal  sweep  problem.  Prediction  methods  involving 
multi-layered  reservoirs  have  generally  resorted  to  ver¬ 
tical  conformance  or  experience  factors  to  account  for 
volumetric  sweepout  by  the  displacing  fluid.  Single 
layer  prediction  methods  have  relied  more  on  model 
studies  for  determining  areal  sweep  efficiency  and  its 
dependence  on  mobility  ratio.  Earlier  prediction  methods 
were  more  concerned  with  areal  sweep  at  breakthrough  than 
the  variation  in  sweep  out  as  a  function  of  water  cut 
following  breakthrough.  It  was  left  to  Dyes  et.al.  (17) 
to  incorporate  this  factor  into  the  water  flood  prediction 
scheme . 

MOBILITY  RATIO 

The  ratio  of  water  mobility  to  that  of  oil  is 
a  primary  variable  in  predicting  water  flood  performance. 
Basically  it  provides  a  measure  of  the  relative  ease  with 
which  the  injected  water  can  displace  oil.  Prediction 
methods  use  mobility  ratio  as  a  correlative  variable 
determined  from  laboratory  studies.  A  related  variable  - 
conductance,  defined  as  the  ratio  of  injection  rate  q  to 
pressure  drop  AP^-accounts  for  the  overall  fluid  conduc¬ 
tivity  of  the  reservoir  and  has  been  utilised  by  Caudle 
and  Witte  ( 18 ) . 
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DISPLACEMENT  MECHANISM 

The  final  primary  variable  to  influence  water- 
flood  performance  is  the  actual  displacement  mechanism 
itself  which  operates  to  force  the  oil  out  of  the  reser¬ 
voir.  Studies  of  this  variable  in  the  laboratory  indicate 
that  for  values  of  injection  rate  lower  than  the  critical*, 
water  flood  is  rate  sensitive  and,  therefore,  the  amount 
of  oil  recovered  varies  with  the  injection  rate.  However, 
for  injection  rates  higher  than  the  critical;  the  perfor¬ 
mance  of  a  water  flood  is  no  longer  influenced  by  the  rate. 

Most  all  prediction  methods  which  account  for  the  displacement 
mechanism  in  an  acceptable  and  comprehensive  manner  use 
some  form  of  Buckley-Leverett  (20) frontal  advance  theory. 

NON-UNIFORMITIES 

It  has  been  said  that  the  reservoir  non-unifor¬ 
mity  effects  are  the  predominent  factors  in  controlling 
the  overall  water  flood  performance.  Therefore,  thoroughly 
evaluated  reservoir  non-uniformities  make  up  the  most 
important  data  in  any  prediction  method.  Two  basic  methods 
for  characterizing  permeability  variations  in  a  reservoir 
are  available. 

*  critical  injection  rate  is  the  rate  obtained  from  the 
scaling  coefficient  proposed  by  Rapoport,  Carpenter  and 
Leas  ( 19 )  . 
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-  the  frequency  distribution  method  ,  and 

-  the  positional  approach. 

The  frequency  distribution  method  involves  ar¬ 
ranging  the  core  analysis  permeabilities  in  order  from 
largest  to  the  smallest  value.  These  are  ordered  regard¬ 
less  of  location  of  the  core  samples  in  the  well.  Then 
permeabilities  are  divided  into  ten  or  more  groups  and 
the  permeabilities  within  any  group  are  averaged  to  get 
the  average  layer  permeability.  This  approach  is  descri¬ 
bed  by  Stiles  (21) . 

The  positional  approach  on  the  other  hand,  was 
first  described  by  Miller  and  Lents  (22)  and  later  by  Elkins 
(23).  This  requires  core  analyses  from  at  least  two  wells. 
The  permeabilities  of  the  core  samples  in  all  wells  one 
foot  below  the  top  of  the  pay  are  averaged  and  similarly 
two  feet  below  the  top,  etc.  As  a  result  of  this  techni¬ 
que  the  high  and  low  permeabilities  are  averaged  together, 
yielding  a  much  less  extreme  permeability  profile  than 
that  obtained  from  the  frequency  distribution  method. 

Summarising  the  goals  of  various  prediction  me¬ 
thods  that  have  so  far  been  proposed,  it  can  be  stated 
that  the  best  prediction  method  is  one  which  includes: 

-  fluid  flow  effects 

-  pattern  effects 

-  heterogeneity  effects. 
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Included  in  the  fluid  flow  effects  are  satura¬ 
tion  gradient  which  can  occur  behind  the  flood  front,  con¬ 
sideration  of  the  effects  of  an  initial  gas  saturation 
and  varying  fluid  conductivity. 

Well  pattern  effects  involve  consideration  of 
areal  sweep  efficiency  at  breakthrough  as  a  function  of 
mobility  ratio  and  the  subsequent  increase  in  areal  sweep 
with  continued  water  injection. 

Heterogeneity  effects  include  consideration  of 
the  non-uniform  nature  of  oil  reservoirs,  permeability 
stratification  etc.  Aside  from  these  three  factors,  in 
the  case  of  communicating  layers,  the  best  prediction 
method  should  also  consider  the  effects  of  cross  flow  due 
to  viscous,  capillary  and/or  gravitational  forces.  Above 
all  it  should  agree  with  the  actual  field  performance. 
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THE  HIGGINS -LEIGHTON  METHOD 

Recently,  Higgins  and  Leighton  (24)  published 
a  numerical  technique  for  tie  calculation  of  fluid  displace¬ 
ment  for  any  irregularly  bounded  porous  medium.  The  me¬ 
thod  employs  the  Buckley-Leverett  (25)  and  Welge  (26) 
concepts  and  is  based  on  the  following  assumptions: 

1.  The  porous  medium  considered  is  homogeneous  and 
has  uniform  permeability  and  porosity. 

2.  A  constant  known  pressure  drop  is  maintained  be¬ 
tween  the  injection  and  the  production  well. 

3.  Gravitational  and  capillarity  effects  are  negligible. 

4.  The  mobility  ratio  of  the  system  is  unity. 

5.  There  is  no  cross  flow  between  the  various  chan¬ 
nels  . 

6.  Gas  saturation  is  identical  in  both  oil  and  water 
banks . 

7.  Fluids  are  incompressible. 

The  approach  used  is  best  illustrated  by  refer¬ 
ence  to  Figure  I,  which  shows  the  isopotential  and  stream¬ 
lines  for  the  five  spot  flooding  network.  Here,  a  potenio- 
metric  model  has  provided  detailed  information  on  the  po¬ 
sitions  of  the  streamlines  and  isopotentials.  It  may  be 
noticed  that  streamlines  (1)  and  (2)  enclose  a  flow  chan¬ 
nel,  and  that  such  a  channel  could  be  divided  into  a  num¬ 
ber  of  cells  of  equal  volume  between  the  injection  and 
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and  producing  wells.  Cells  along  a  particular  channel 
would  then  be  of  varying  length  and  differing  effective 
cross  -  sectional  area.  The  resistance  to  flow  because 
of  the  shape  and  size  of  the  cells  is  represented  by  G, 
the  shape  factor  which  is  a  measure  of  the  geometric  re¬ 
sistivity.  In  equation  form,  the  shape  factor  is  defined 
as  : 

G  =  K  (!) 

Higgins,  Boley,  and  Leighton  (27)  have  calcu¬ 
lated  shape  factors  and  channel  volumes  for  the  direct 
line  drive, five  spot,  seven  spot,  and  staggered  line  drive 
well  patterns. 

In  the  determination  of  the  shape  factors  for 
the  five  spot  pattern  advantage  was  taken  of  its  inherent 
symmetry  and  only  one-eighth  of  a  five  spot  was  used.  Each 
channel  was  divided  into  40  cells  of  equal  volume.  This 
number  of  cells  was  chosen  on  the  basis  of  results  of  an¬ 
other  investigation  which  indicated  that  40  cells  were 
more  than  a  sufficient  number  to  insure  good  convergence 
properties  during  the  calculation  of  a  water  flooding 
performance  (28). 

Where  potentiometric  data  are  available  on  an 
irregularly  shaped  well  pattern,  a  careful  construction 
of  the  flowlines  and  a  division  into  cells  will  result  in 
specific  values  for  the  shape  factor,  G.  In  the  case  of 
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irregularly  shaped  cells,  the  shape  factor  is  defined  as: 


G 


+  L 


+  L, 


(2) 


where  L  =  length,  the  subscripts  1,  r,  t  and  b  refer 
to  the  left  side,  right  side,  top  and  bottom  of  the  cell 
and  the  channel  thickness  is  unity. 

If  the  cells  are  extremely  irregular,  Henley 

(29)  has  proposed  a  technique  where  the  shape  factor  is 
determined  by  using  inscribed  circles.  Higgins,  et.al. 

(30)  have  developed  a  computer  calculation  technique 
whereby  volumes  of  the  channels  and  the  shape  factors  may 
be  determined  from  data  taken  from  a  potentiometric  mo¬ 
del  study  for  any  well  spacing  pattern. 

Kufus  and  Lynch  (31)  have  shown  that  the  aver¬ 
age  permeability  of  a  cell  times  the  shape  factor  deter¬ 
mines  the  resistance  to  flow.  Then,  from  Darcy's  law 
for  a  linear  horizontal  system,  the  producing  rate  at  the 
beginning  of  the  flood  will  be: 


where 

j 


ap 


‘roIW 


G. 

l 


(3) 


index  of  increments  since  beginning  of  the 
flood 


i 


=  cell  number  index 
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N  -  total  number  of  cells  in  each  channel 
^roIW  =  permeability  to  oil  at  irreducible  water 
saturation . 

After  n  number  of  cells  have  been  invaded,  the  producing 
rate  prior  to  displacing  phase  (water)  breakthrough  is: 

K  AP 
a 

q  .  =  - - - 

n  N 


where  K  =  mean  permeability  in  a  particular  cell. 

In  order  to  solve  equation  (4),  it  is  necessary 
to  determine  mean  values  for  the  oil  and  water  relative 
permeabilities  in  each  cell  when  j  number  of  cells  have 
been  invaded.  This  is  accomplished  by  calculating  the 
f  versus  relationship  from  the  following  equation: 


f  = 
w 


1 


,  ,  K  y 

1  +  o  w 


(5) 


K  y 
w  o 


The  mean  permeability  of  the  first  cell  to  water  at 
the  instant  that  water  breakthrough  occurs  at  the  end  of 
the  cell  is  determined  by  dividing  the  slope  of  the  f 

w 

versus  curve  at  breakthrough,  f'^,  by  the  area  under 
a  curve  of  the  resistance  to  water  flow,  1/K  ,  versus 

JL  W 

the  slope,  f'.  This  area  is  that  up  to  the  value  of  the 


A  derivation  of  equation  (4)  is  presented  in  Appendix  A, 
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slope,  f',  corresponding  to  the  total  volume  behind  the 
flood  front.  Figure  II  illustrates  these  relationships. 

When  the  water  front  has  progressed  to  the  end 
of  the  second  cell,  the  mean  permeability  to  water  in  the 
first  cell  is  one-half  of  the  f'^r  value,  divided  by  the 
area  under  the  resistance  curve,  up  to  the  abscissa  va¬ 
lue  of  one-half  f'jDr*  The  mean  permeability  to  water 

in  the  second  cell  then  becomes  one  half  f',  ,  divided 

br 

by  the  remaining  area  under  the  resistance  curve.  The 
following  equation  presents  the  procedure  outlined  in 
an  analytical  form: 

K 


br 


rw .  = 

1 


. -  N  -i 

All  areas  under  1/K 


(6) 


vs .  f '  curve  in  n 
cell 


th 


rw 


A  similar  equation  can  be  readily  written  for 
the  mean  permeability  to  oil.  When  breakthrough  occurs 
(the  subordinate  phase)  at  the  end  of  one  of  the  channels, 
equation  (4)  must  be  modified  to  include  a  water  producing 
rate : 


(q„  + 


qw)j 


(V) 


*  A  derivation  of  equation  (6)  is  presented  in  Appendix  A. 


WATER  SATURATION,  RELATIVE  PERMEABILITY 

TO  WATER,  AND  RESISTANCE  VERSUS 

THE  SLOPE  TERM,  f'  (AFTER  HIGGINS-LE1GHTON) 
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The  mean  permeabilities  to  water  and  oil  in  each 
cell  are  now  modified  due  to  the  decrease  of  the  slope  of 
the  f  versus  S  curve,  as  water  saturations  increase 

W  W 

beyond  that  at  water  breakthrough.  In  equation  form,  the 
mean  permeability  to  water  in  any  given  cell  becomes: 


K 


br  -  m 


(8) 


rw . 
1 


All  areas  under  1/K 
vs .  f '  curve  in 
cell. 


rw 


where  f'  is  the  first  derivative  of  the  f  ver- 

br  -  m  w 

sus  relationship  from  the  fractional  flow  equation  at 
the  m"^  increment  after  water  breakthrough.  The  equation 
for  mean  oil  permeability  takes  a  similar  form.  After 
the  rrr11  increment  of  production  has  occurred  past  water 
breakthrough,  the  water  -  oil  ratio  is: 

(9) 


(WOR) . 

j  =  m+N 


K  y 
rw  o 


y  K 
w  ro 


The  oil  and  water  rates  can  now  be  calculated  by  solving 

equations  (7)  and  (10)  simultaneously. 

qw 

_  m  +  N  ,  v 

qoj  (subordinate)  = - - —  (10) 

(W0R>j=  m  +  N 

Since  the  mean  permeability  values  are  used  to 
calculate  the  oil  and  water  producing  rates  in  the  subor¬ 
dinate  phase,  the  instantaneous  producing  rates  at  the 
well  may  be  approximated  as  follows: 
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(qo>j  -  1/2 


Vj  +  ^o’j  +  1 


(ID 


(qw}j  -  1/2 


(qw' j  +  ^qw^  j  +  1 


(12) 


Here,  the  index  is  equal  to  m  +  N,  or  the  m  increment 
of  production  past  water  breakthrough,  where  N  is  the  to¬ 
tal  number  of  cells  in  each  channel. 

Since  each  channel  was  divided  into  N  cells  of 
equal  volume,  the  average  volume  of  oil  produced  from  each 
cell  during  the  primary  phase  may  be  written  as: 


V 


ol 


V  (s 
p  w 


Br 


S  wIW^ 


(13) 


N 


where 


V 


w 


Br 


wIW 


=  total  pore  volume  of  the  channel 

=  average  water  saturation  behind  the  front 
at  time  of  breakthrough. 

=  connate  or  irreducible  water  saturation. 


During  the  subordinate  phase,  the  Welge  equation  applies. 


The  following  form  of  the  equation  is  convenient: 

1  -  f 


S  S  ,  w 

w  =  w  +  m 

m  m  - 


(14) 


f  ' 


where  S  =  average  water  saturation  in  a  channel,  and 
w 

the  index  m  refers  to  the  number  of  increments  after  ini¬ 
tial  breakthrough  of  water  at  the  producing  well. 
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As  in  a  linear  calculation  of  the  frontal  dis¬ 
placement  mechanism,  the  oil  produced  during  each  separate 
step  of  the  subordinate  phase  may  be  calculated  as  follows: 


m  +  N 


(S 


w 


m  + 


w 


m 


(15) 


The  time  for  each  step  can  then  be  readily  determined  from 
the  following  equation: 


t  . 
3 


(16) 


Utilizing  equations  (1-16)  performance  is  cal¬ 
culated  for  each  channel  individually  until  a  preset  wa¬ 
ter-oil  ratio  is  exceeded.  The  performance  of  these 
channels  is  then  combined,  using  the  time  for  channel  4 
as  a  reference,  to  get  the  total  production  history.  The 
performances  of  channels  3,  2,  and  1  from  their  ending 
times  to  the  ending  time  of  the  4th  channel  are  obtained 
by  linear  extrapolation. 

It  may  be  noted  that  in  this  calculation  tech¬ 
nique,  an  areal  coverage  factor  has  not  been  used.  Such 
a  procedure  is  not  required,  since  breakthrough  at  the 
end  of  each  particular  cell  in  each  channel  defines  the 
position  of  the  front. 
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THE  MODIFIED  SCHEME 

The  theory  and  equations  used  in  this  scheme 
are  basically  the  same  as  those  used  by  Higgins  and  Leighton. 
It  may  be  recalled  that  the  combined  performance  in  the 
Higgins-Leighton  method  is  obtained  by  combining  the  per¬ 
formances  of  various  channels  on  the  same  time  basis. 

This  introduces  an  error  in  the  prediction  since  straight' 
line  interpolation  and  extrapolation  is  used  to  calculate 
the  production  history  at  intermediate  time  intervals. In 
order  to  circumvent  this  problem,  all  the  four  channels 
are  flooded  simultaneously  in  the  modified  approach. 

This  is  accomplished  by  using  channel  one  as  a  reference 
and  invading  complete  cells,  one  at  each  step,  in  this 
channel.  The  time  required  to  invade  a  complete  cell  in 
channel  one  is  then  used  to  locate  the  position  of  the 
front  in  the  remaining  channels. 

Although  the  Higgins-Leighton  method  assumes 
that  there  is  no  cross  flow  between  various  channels,  no 
attempt  was  made  by  them  to  verify  this  assumption.  In 
attempting  to  flood  the  four  channels  simultaneously  in 
the  suggested  modification,  an  iterative  scheme  was  set 
up  to  calculate  the  position  of  the  flood  front  in  chan¬ 
nels  4,  3,  and  2  relative  to  channel  1  by  using  no  cross 
flow  between  the  channels  as  a  convergence  criterion. 
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It  may  be  recalled  that  Higgins-Leighton  method 
is  applicable  to  only  one  case  where  a  constant  pressure 
drop  is  maintained  between  the  injection  and  the  produc¬ 
tion  well.  In  order  to  circumvent  this  limitation,  the 
modified  approach  considers  the  following  cases: 

-  constant  pressure  drop  between  the  injection  and 
the  production  well 

-  constant  injection  rate 

-  variable  injection  rate  and  pressure  drop. 

CONSTANT  PRESSURE  DROP  BETWEEN  THE  INJECTION  AND  THE 

PRODUCTION  WELL 


Primary  Phase 

In  order  to  prevent  cross  flow  between  the  va¬ 
rious  channels,  the  pressure  drop  between  the  injection 
and  the  production  well  at  each  step  in  each  channel  must 
be  identical.  However,  since  the  resistance  to  flow  changes 
as  front  advances  towards  the  production  well,  the  water 
rates  in  each  channel  must  change  to  account  for  the  cor¬ 
responding  change  in  resistance.  Therefore,  the  iterative 
scheme  should  be  such  as  to  account  for  the  variations 
in  water  rates  with  time.  Consequently,  the  following 
steps  were  used  in  implementing  the  iterative  scheme. 

1.  Before  initiating  the  flood,  flow  resistances 
in  each  channel  are  computed  by  using  the  fol¬ 
lowing  equation: 
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N 


i  =  *1-4 


y  G. 
o  k 


K 


(17) 


roIW 


I<  =  1 


where  i  denotes  the  channel  number 
2.  These  flow  resistances  are  then  used  to  calculate 
initial  water  rates  in  the  various  channels,  us¬ 
ing  the  following  realtionship : 


AP  x  K 


R. 

l 


(IB) 


i  =  1  -*  4 


3.  This  information  is  then  used  to  calculate  the 
time  step  required  to  invade  one  complete  cell 
in  channel  1.  This  has  been  accomplished  by 
using  the  following  scheme: 


Volume  of  oil  produced  when  = 
one  complete  cell  is  invaded 


„  (S  -  S  TT7) 
VP  WBR  WlW 

NCELLS 


(19) 


where 


NCELLS 

V 


w 


BR 


wIW 


=  total  number  of  cells  in  a  channel, 

=  pore  volume  of  a  channel, 

=  average  water  saturation  behind  the 
front  at  the  time  of  breakthrough,  and 

=  connate  or  irreducible  water  saturation 


Assuming  that  the  flow  rate  calculated  in  step 
2  remains  constant  for  this  time  period,  the  vol¬ 
ume  of  water  injected  would  then  be  equal  to 
qi  x  ATi,  where  ATi  is  the  time  step. 
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In  the  primary  phase  since  the  volume  of  oil 
produced  is  equal  to  the  volume  of  water  injected, 


ATi  = 


,7  (S  S  TT7) 
V  •  w  -  wIW 
P  BR 

qi  X  NCELLS 


(20) 


4.  Using  this  time  step,  the  flood  front  positions 
in  channels  2,  3,  and  4  are  determined  by  using 
the  following  equation:  {21) 


No.  of  cells  invaded 
(in  channels  2,3,4) 


Avg .  w . r .  X  ATi  X  NCELLS 


c  ) 
~wIW; 


5.  Once  the  flood  has  been  initiated  and  the  first 
cell  in  channel  1  has  been  invaded,  oil  and  wa¬ 
ter  both  flow  behind  the  front,  whereas  only 
oil  flows  ahead  of  the  front.  Consequently,  the 
total  resistance  to  flow  in  each  channel  is  the 
sum  of  the  resistance  due  to  two-phase  flow 
and  that  due  to  single-phase  flow.  This  is  re¬ 
presented  by  the  following  equation: 


where  n 


number  of  cells  invaded. 
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6.  Using  the  flow  rates  calculated  in  sten  2  as  an 
initial  approximation,  the  pressure  drops  in  each 
channel  are  calculated  according  to  equation  (18). 
These  pressure  drops  are  then  compared  against 
the  one  read  into  the  computer  initially.  If  the 
two  values  differ  more  than  that  allowable  by 

the  error  criterion,  the  flow  rates  in  each  chan¬ 
nel  are  modified  using  the  following  equation: 

AP 

a  =  n  v  original  (23) 

i  4  i  .  AP  .  i  j 

new  previous  calculated 

i  =  1  4 

7 .  Steps  3  to  6  (inclusive)  are  repeated  until  the 
error  criterion  is  satisfied. 

8.  Steps  1  to  7  (inclusive)  are  repeated  until  the 
40  cells  in  channel  1  are  invaded. 


Subordinate  Phase 


After  breakthrough  it  is  believed  that  the  average 


water  saturation  behind  the  flood  front  gradually  increases 

d  f  w 

and  that  the  derivative  gradually  decreases.  In  or- 

w 

der  to  simulate  these  conditions  in  the  modified  version 


dfw 
ds 


was  decreased  in  steps  according  to  the  following 


w 


equation : 


NEND  X  40 _ ________ 

(40  +  additional  steps) 


NSEND 


(24) 
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where 

NEND  =  Total  number  of  divisions  in  f'grid  used 
to  calculate  mean  permeabilities  in  the 
various  cells,  and 

NSEND  =  Total  number  of  divisions  remaining  in 
f'grid  at  any  step. 

This  information  is  then  used  to  compute  the 
water-oil  ratio  according  to  the  procedure  outlined  in 
the  Higgins-Leighton  method. 

The  nature  of  the  channels  is  such  that  the 
resistance  to  flow  is  greatest  in  channel  4  and  decreases 
in  order  of  magnitude  from  channel  4  to  channel  1.  As 
a  result  when  channel  1  has  broken  through  channels  2, 

3,  and  4  are  only  partly  invaded. 

Since  in  the  modified  approach  all  the  four 
channels  are  flooded  simultaneously  and  different  schemes 
are  used  for  treating  primary  and  subordinate  production, 
it  is  necessary  to  monitor  the  position  of  the  flood 
front  in  each  channel.  After  channel  1  has  broken  through, 
the  position  of  the  flood  front  in  the  remaining  channels 
is  observed  in  order  to  determine  when  to  use  the  subor¬ 
dinate  phase  calculations  in  other  channels. 

Finally,  the  flood  is  terminated  when  the  com¬ 
bined  instantaneous  water-oil  ratio  from  all*  the  four 
channels  exceeds  a  preselected  value. 


26 


CONSTANT  INJECTION  RATE 


Primary  Phase 


In  order  to  maintain  a  constant  total  injection 

rate,  and  also  to  satisfy  the  no  cross  flow  criterion, 

the  following  conditions  must  be  fulfilled: 

the  sum  of  the  injection  rates  of  the  various 
channels  at  each  step  must  equal  the  total  injec¬ 
tion  rate  read  initially  into  the  computer;  and, 

to  satisfy  the  no  cross  flow  criterion,  pressure 
drop  between  the  injeciton  and  the  production  well 
at  each  time  step,  in  each  channel,  must  be  identi 
cal.  In  view  of  equation  (18),  this  can  be  stated 
mathematically  in  terms  of  the  expression: 


qlRl  _  q2R2  _  q3R3  _  q4R4  _  qR 


K 


K 


K 


K 


K. 


or 


qlRl  =  q2R2 


=  q3R3  =  q4R4  =  qR 


(25) 


(26) 


In  order  that  these  requirements  be  satisfied, 
the  iterative  scheme  employed  consists  of  the  following 
steps : 

1.  First  the  inital  resistance  to  flow  in  each  chan¬ 
nel  is  calculated  using  equation  (17) . 

2.  The  initial  water  rates  in  all  the  channels  are 
then  determined  using  equation  (18) . 

3.  As  soon  as  the  first  cell  in  channel  1  is  invaded 
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the  time  step  and  the  position  of  the  flood  front 
in  the  remaining  channels  is  computed  using  equa¬ 
tions  (20)  and  (21)  . 

4.  Next,  resistances  in  all  channels  are  determined 
using  the  position  of  the  flood  front  in  each 
channel  as  the  boundary  between  the  two-phase 
and  the  single -phase  region.  Equation  (22)  is 
used  for  this  purpose. 

5.  The  average  of  the  product  of  flow  rate  and  resis¬ 
tance  in  all  the  four  channels  is  then  calculated 
by  means  of  the  following  equation: 


PROD 


(27) 


6.  Following  this,  a  check  is  made  to  determine 

whether  or  not  the  product  of  flow  rate  and  resis¬ 
tance  in  each  channel  agrees  with  PROD  within 
the  limits  of  a  chosen  error  criterion.  The  fol¬ 
lowing  equation  is  used  in  making  this  comparison: 


CHECK 


(q.R.  -  PROD) 

l 

PROD 


(28) 


7.  If  the  magnitude  of  CHECK  exceeds  the  error  cri¬ 
terion,  the  flow  rates  in  each  channel  are  modified 
using  the  following  equation: 


q 

i 


new 
.  th 


=  q 


r.  X  R. 

Ll  .  l 

previous 


PROD 


i  previous 


PROD 


channel 


(29) 
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8.  The  injection  rates  in  all  the  channels  are  then 

added  together  accordina  to  equation: 

4 

SUM  =  y  qi  (30) 

i  =  1 

The  quantity  SUM  is  then  compared  to  the  total 
injection  rate,  Q  ,  which  was  read  into  the 
computer  at  the  start  of  the  program.  This  com¬ 
parison  is  made  through  the  use  of  the  equation: 


RERROR  = 


SUM  -  Qr 

Qt 


(31) 


9.  If  the  magnitude  of  RERROR,  is  greater  than 
some  preselected  value,  the  flow  rates  in  each 
channel  are  modified  by  using  the  following 


equation  and. 


=  q 


x 


'new 


'previous 


qt 

SUM 


(32) 


steps  3  to  8  (inclusive)  are  repeated. 

10.  Once  the  conditions  CHECK  and  RERROR  are  satis¬ 
fied,  the  next  cell  in  channel  1  is  invaded. 

11.  Steps  3  to  10  (inclusive)  are  then  repeated  un¬ 
til  40  cells  in  channel  1  have  been  invaded. 

12.  Once  channel  1  has  broken  through  the  position 
of  the  flood  front  in  the  remaining  channels  is 
monitored  as  in  the  constant  pressure  drop  case. 
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Finally,  the  flood  is  terminated  when  the  com¬ 
bined  instantaneous  water-oil  ratio  from  all  the  four 
channels  exceeds  a  preselected  value. 

VARYING  INJECTION  RATE  AND  PRESSURE  DROP 

The  iterative  scheme  in  this  case  is  identical 
to  the  one  employed  in  the  constant  injection  rate  case. 
However,  since  the  total  injection  rate,  Q^,  may  change 
with  time,  an  interpolation  scheme  is  required  to  deter¬ 
mine  this  injection  rate  for  cumulative  times  calculated 
at  each  step  from  time  versus  injection  rate  history. 

This  value  of  injection  rate  determined  at  each  step  is 
then  used  in  place  of  QT  in  the  iterative  scheme  described 
for  the  constant  injection  rate  case. 

A  detailed  listing  of  the  computer  program 
used  in  the  modified  version  is  presented  in  Appendix  B. 

RECENT  WORK  OF  HIGGINS 

Subsequent  to  the  development  of  the  modifica¬ 
tions  just  described  Higgins  presented  (32)  a  method  which 
extends  the  Higgins-Leighton  method  to  cover  the  constant 
rate  case.  His  approach  seems  to  differ  considerably 
from  the  one  suggested  in  the  present  investigation.  Although 
most  of  the  working  details  are  not  obvious  from  the  pre¬ 
print  available  at  the  time  of  writing  this  report,  the 
following  points  can  still  be  emphasized  to  indicate  the 


differences  involved. 


■ 
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1.  The  production  history  for  the  constant  pressure 
drop  case  is  a  pre-requisite  for  the  calculation 
of  the  production  history  for  constant  rate  case. 

2.  The  instantaneous  oil  rate  is  obtained  by  divi¬ 
ding  the  constant  injection  rate  by  the  instant¬ 
aneous  injection  rate  under  constant  pressure 
and  multiplying  the  quotient  by  the  instantaneous 
oil  rate  computed  in  the  constant  pressure  drop 
case . 

3.  The  time  steps  for  the  constant  injection  rate 
case  are  obtained  by  multiplying  the  fractional 
pore  volume  of  water  injected  by  the  number  of 
barrels  of  pore  volume  in  the  five  spot  pattern, 
and  then  dividing  the  product  by  the  constant 
injection  rate.  It  may  be  noted  that  this  is 
on  the  basis  of  the  total  pattern  and  not  any 
individual  channel  which  is  the  case  in  the 
modified  scheme. 

4 .  The  method  does  not  use  an  iterative  scheme  and 
no  attempt  has  been  made  to  check  whether  the 
instantaneous  rate  used  in  a  channel  at  a  parti¬ 
cular  time  step  complies  with  the  no  cross  flow 
assumption . 

Due  to  limitations  of  time  a  detailed  examination 
of  Higgins  modification  (33)  has  not  been  conducted.  Con¬ 
sequently,  his  method  and  the  one  reported  herein  have  not 


been  compared. 


. 
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RESULTS  AND  DISCUSSION 

The  modified  scheme  was  applied  to  a  five  spot 
pattern  studied  experimentally  by  Douglas,  Peaceman,  Rach- 
ford  (34)  and  theoretically  by  Higgins  and  Leighton  (35). 

The  results  obtained  are  tabulated  in  Appendix  C  and  com¬ 
pared  in  Figures  III  and  IV. 

Figure  III  presents  a  variation  in  pore  volumes 
of  oil  produced  as  a  function  of  pore  volumes  of  water 
injected.  It  may  be  observed  that  until  water  injected 
exceeds  2.0  pore  volumes,  the  results  obtained  from  the 
modified  approach  are  closer  to  the  laboratory  data  as 
compared  to  the  Higgins-Leighton  results.  It  may  also  be 
noted  that  the  modified  scheme  yields  identical  results  in 
both  cases  where  injection  rate  or  pressure  drop  is  maintained 
constant.  This  can  be  explained  on  the  basis  that  since 
maintaining  the  injection  rate  or  pressure  drop  constant 
affects  only  the  time  rate  of  oil  production,  the  results 
must  be  compared  against  the  time  parameter. 

Figure  IV  presents  a  variation  in  instantaneous 
water-oil  ratio  as  a  function  of  time.  It  can  be  clearly 
seen  that  the  modified  approach  yields  different  results 
for  both  cases  where  injection  rate  or  pressure  drop  is 
maintained  constant.  Since,  instantaneous  water-oil  ratio 
versus  time  histroy  was  not  reported  by  Douglas  et.al. 

(36) ,  the  results  of  the  modified  approach  could 
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not  be  compared  against  the  laboratory  data.  However,  a 
comparison  with  Higgins-Leighton  results  is  presented  and 
reveals  that  there  are  discontinuities  in  the  results  ob¬ 
tained  from  the  modified  approach.  These  discontinuities 
represent  breakthrough  of  each  channel.  Since  the  results 
obtained  from  the  Higgins-Leighton  method  were  averaged 
before  plotting,  these  discontinuities  are,  therefore, 
absent  in  the  data  calculated  from  their  method. 

The  modified  scheme  requires  relative  permea¬ 
bility  versus  saturation  curves  as  its  basic  data.  Within 
the  limitations  of  time  provided  for  conducting  the  pre¬ 
sent  investigation  representative  relative  permeability 
versus  saturation  data  for  the  laboratory  model  used  in 
the  present  investigation  could  not  be  obtained  and  thus 
the  attempts  at  checking  compatibility  of  laboratory  re¬ 
sults  proved  in  vain. 
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CONCLUSIONS 

The  conclusions  drawn  as  a  result  of  this  study 

are : 

1.  The  Higgins-Leighton  method  has  been  modified 
and  extended  to  predict  the  performance  of  wa¬ 
ter  floods  for  the  following  cases: 

-  constant  injection  rate 

-  constant  pressure  drop  between  the  injection 
and  the  production  well 

-  variable  injection  rate  and  pressure  drop. 

2.  The  modified  approach  was  tested  on  the  labora¬ 
tory  data  reported  by  Douglas,  et.al.  (37). 

The  agreement  obtained  between  the  experimental 
and  calculated  results  was  well  within  the  de¬ 
sired  degree  of  accuracy. 

3.  The  use  of  an  iteration  scheme  in  the  proposed 
method  ensures  that  no  cross  flow  occurs  bet¬ 
ween  the  various  channels.  Thus,  although  the 
shape  factors  used  are  for  unit  mobility  ratio, 
possibilities  of  obtaining  better  correspondence 
between  laboratory  results  and  calculated  pre¬ 
diction  for  adverse  mobility  ratios  look  favour¬ 


able  . 
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APPENDIX  A 


Derivation  of  the  equation  used  to  calculate 
oil  rate  in  the  primary  phase  after  first 
cell  has  been  invaded. 

Derivation  of  the  equation  to  determine  the 
mean  permeability  in  a  cell. 
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EQUATION  TO  COMPUTE  OIL  RATE  IN  THE  PRIMARY  PHASE  AFTER 

FIRST  CELL  INVASION 


By  Darcy’s  Law  for  horizontal  flow 

o  =  ~Ka  KrAp 
L  y  L/A 


If  one  defines  the  shape  factor 
G  =  L/A 


(1) 


than 


Q 


-Ka  Kr  AP 
y  G 


(2) 


AP 
■■  ■  ■  1  ■ 

pG 

Ka  Kr 


By  analogy  to  Ohm's  Law  one  may  write  this  equation 


as  : 


Q  = 


-AP 

R 


(3) 


where  R  =  Tr- — and  is  called  the  resistance  to 

Ka  Kr 

flow.  Thus  for  a  simultaneous  flow  of  oil  and  water,  the 


resistance  to  flow  for  each  phase  would  be  given  by: 

R 


G 


oil  phase  Ka  Kr 


o 


R 


water  phase  = 


y  G 
w 


Ka  Kr 


w 


(4) 


(5) 
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After  the  invasion  of  the  first  cell,  there 
is  two-phase  flow  behind  the  flood  front  in  the  first 
cell  consisting  of  oil  and  water,  and  single-phase  flow 
ahead  of  the  front  consisting  only  of  oil.  Therefore, 
the  total  resistance  to  fluid  flow  will  comprise 


R  =  R  +  R 

total  two-phase  single-phase 


where 


R 


two-phase 


Ka 


Kr 


w 


Kr 


+ 


R 


w  o 

G 


,  and 


single-phase  Ka  Kr 


oIW 


Thus , 


R 


i  =  1 
Y 


N 


Gi 


total 


Ka 


1  =  0 


Kr  Kr 


w 


w 


+ 


1  = 


(6) 


G.1 


Ka  Kr 


oIW 


where, 

Kr  =  relative  permeability  to  oil  at  irreducible 
oIW  ' 


water  saturation. 


Substituting  this  value  of  ©gu^tion  (3) 

yields • 


Q 


°D 

(primary) 


- 
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con ' t .  . 


EQUATION  TO  DETERMINE  THE  MEAN  PERMEABILITY  IN  A  CELL 


The  resistance  to  flow  of  fluids  in  a  linear, 
horizontal  system  is  described  by  the  following  equation: 


where 

Realizing 


q 


R  = 
that 


-AP 

R 

G 

Ka  Kr 
this 


of  the  flood  front, 
value  R  such  that, 


(1) 

and  is  termed  as  the  resistance, 
resistance,  R,  varies  with  the  position 
it  is  convenient  to  define  an  average 


R 


L 

/  RdX 
o _ 

L 

f  dX 
o 


(2) 


but 

R 


y  G 
Ka  Kr 


L 

/ 

o 


y  G 

Ka  Kr 


dX 


L 

f 

o 


dX 


However,  for  the  problem  at  hand  p,  G,  and  Ka 


are  independent  of  X. 
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R  = 


y  G 

K 


a 


L 

/ 

o 

L 

/ 

o 


dX 

Kr 


dX 


(3) 


However,  one  may  also  write, 


R  = 


y  G 


Ra  Kr 


but,  y ^ ,  G,  and  Ka  remain  constant 


R  = 


y  G 
Ka  Kr 


(4) 


so  that, 
y  G 


Kr  = 


(5) 


Ka  R 


Substituting  for  R  from  equation  (3)  one  obtains 


Kr  = 


yG  Ka 

Ka  yG 


Kr  = 


& 


dX 


dX 
f  Kr 
o 


L 

/  dX 
o 


dX 

/  Kr 
o 


(6) 


Using  Buckley-Leverett ' s  material  balance  equation 


X  = 


4.  df 

qt  w 


0A  ds 


w 


=  gt  f. 

0A  W 


(7) 
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Differentiating  both  sides 


dX  =  31  df' 

0A  w 


(8) 


(Assuming  qt ,  0,  and  A  are  constant) 


Substituting  equation  (8)  in  equation  (6)  one  obtains 

L 

/df  'w 


c-S'  §  -r 


rv~  df  ' 

J  Kr  w 


Simplifying 


Kr  = 


L 

f  df  ’ 
o 


w 


L  1 

r  F  df ' 

/  Kr  w 

o 


(9) 


Applying  boundary  conditions, 


.  .  Kr 


At  x  =  o  ; 

f  ' 

w 

=  o 

x  =  L  ; 

f  ' 

=  f 

w 

f ' 

w 

br 

w. 


br 


mean 


w. 


br 


^  Kr 


1  df' 


w. 


br 


(10) 
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Thus,  for  oil  and  water  phases  this  equation  will  take 
the  following  form: 


K 


w 


ro 


br 


mean 


w,  i 
br  _1_ 

/  Kr 


(ID 


df 


w 


br 


f ' 


K. 


w 


rw 


br 


mean 


w 


(12) 
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o 


Kr  df ' 
w  w. 
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Basic  data  required  to  use  the  method 
Nomenclature  for  computer  program 
Listing  of  the  computer  program 
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BASIC  DATA  REQUIRED  TO  USE  THE  METHOD 

The  following  data  are  required  to  predict  the 
performance  of  any  pattern  using  the  modified  scheme: 

1.  Relative  permeability  versus  saturation  curves. 

2.  Pore  volumes  and  shape  factors  for  each  channel. 

3.  Initial  water  saturation  (connate  or  irreducible) 
and  irreducible  oil  saturation. 

4.  Maximum  water  -  oil  ratio  at  which  cut-off  is 
desired . 

5.  Pressure  drop  between  the  injection  and  the 
production  well  for  a  constant  pressure  drop  case 

6.  Injection  rate  for  a  constant  injection  rate  case 

7.  Time  versus  injection  rate  history  where  both 
injection  rate  and  pressure  drop  vary. 

8.  Absolute  permeability  of  the  medium. 

9.  Oil  and  water  viscosities. 


NOMENCLATURE  FOR  COMPUTER  PROGRAM 


ALPA  (  ) 

-  Left  part  of  a  unit  area,  Af'  times  recipro¬ 
cal  of  permeability. 

ALPAPT 

=  Left  part  on  the  base  of  a  unit,  Af'. 

AVGSAT 

=  Average  saturation  behind  the  flood  front 
after  breakthrough 

AVPERM 

=  Name  of  subroutine  to  determine  the  average 
permeability  to  water. 

AVPER2 

=  Name  of  subroutine  to  determine  the  average 
permeability  to  water  for  fractional  cells. 

D1  (  ) 

First  difference.  used  in  Newton  and  Stir- 

D2  (  ) 

=  Second  difference.  ling  equations  for  inter- 

polating  for  oil  permea- 

D3  (  ) 

=  Third  difference.  bility  values  in  lookup 

table . 

Dll  (  ) 

— 

D22  (  ) 

=  Same  as  above  but  for  water  permeability 

D33  (  ) 


DEL 

=  Width  of  an  interval  into  which  the  saturation 
range  between  irreducible  water  and  oil  is 
divided;  used  for  determining  s  at  break¬ 
through  . 

DEL2 

=  Width  of  an  interval  into  which  the  satura¬ 
tion  range  between  f '  at  breakthrough  and 
irreducible  oil  is  divided. 

DEL3 

=  Width  of  a  division  into  which  range  of  f' 
between  f '  at  breakthrough  and  f '  at  the 
irreducible  oil  saturation  is  divided. 

DELN 

=  Number  of  divisions  into  which  the  water 
saturation  range  between  the  interstitial 
water  and  irreducible  oil  saturation  is  to 
be  divided  in  order  to  attain  the  f  function 
and  the  derivative. 

DJ 

=  Decimal  expression  of  J. 
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DN 

=  Decimal  expression  of  N. 

DNCELL 

=  Decimal  expression  for  the  integer  NCELLS. 

DNDIVS 

=  Decimal  expression  of  NDIVS. 

DS 

=  Number  of  divisions  between  S  at  breakthrough 

and  S  at  irreducible  oil. 
w 

DSMN1 

=  Decimal  expression  of  NDSMNl . 

ELTIME 

=  Total  elapsed  time. 

ENDPO 

=  Last  permeability  to  oil  in  the  permeability 
lookup  table;  irreducible  oil. 

ENDPW 

=  Last  permeability  to  water  in  the  permeability 
lookup  table;  irreducible  oil 

ENDSW 

=  Water  saturation  at  irreducible  oil,  same 
value  as  SWIO. 

F  (  ) 

=  Fractional  flow  of  water  equals  1/(1  +  k  y  / 
i  \  r;  ^  '  ro  w' 

k  y  )  =  f 

rw  o  w . 

FA  (  ) 

=  Average  of  two  f^  values. 

FCELL  (  ) 

=  Retains  the  value  of  FRCELL  (  ) . 

FACTOR 

=  Proportionality  factor  to  interpolate  for 
saturation  at  breakthrough. 

FB 

=  f  at  breakthrough  adjusted  by  interpolation, 
w 

FLORTE  (  )  = 

=  Subroutine  used  to  interpolate  for  QT  in  a 
case  where  both  pressure  drop  and  Injection 
rate  vary. 

FLOCHK  (  )  = 

=  Subroutine  used  to  check  cross  flow. 

FP  (  ) 

=  Derivative  of  f  . 

w 

FPB 

=  f '  at  breakthrough. 

FPP  (  ) 

=  f '  ;  it  is  the  same  as  FP  but  the  computer 
program  requires  a  double  indexing  during 
part  of  the  calculation. 

FRCELL  (  )  = 

=  Total  number  of  cells  invaded  including 
fractional  cells. 
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G  (  ) 

=  Shape  factor,  L/A,  which  is  a  function  of 
the  geometry  of  the  system. 

GEMTRY 

=  Subroutine  for  geometry  -  shape  factors  and 
volumes  -  of  the  system. 

GGG 

=  Single  shape  factor  used  for  inear  case  only. 

H 

=  Difference  between  succeeding  X  values;  used 
in  STIRLG  subroutine. 

I 

=  Index  for  channel  number. 

IGEOM 

=  Test  if  linear  or  nonlinear  geometrical  shape 
factors  are  to  be  read  in  or  previous  shape 
factors  to  be  reread  and  used.  If  zero,  val¬ 
ues  from  previous  case  are  used,  if  1,  non¬ 
linear  values  are  read,  if  2,  linear  value 
is  used. 

II 

=  Talley  for  channel  number. 

IPAVPM 

=  Print  out  in  the  AVPERM  mean 
permeabilities . 

If  zero. 

IPSETC 

=  Print  out  Sw ' s ,  krw ' s ,  krQ ' s ,  do  not  print; 

f  's,  and  f^'s,  at  selected  positive, 

w  .  ,  ,  print, 

interval . 

1ST 

=  Number  of  channels. 

J 

=  Index. 

JJ 

=  Index. 

JPNTM1 

=  JPOINT  -  1. 

JPNTM2 

=  JPOINT  -  2. 

JPNTM3 

=  JPOINT  -  3. 

JPOINT 

=  Ending  index  of  permeability  data  points 
taken  from  empirical  curve. 

K 

=  Index. 

KABS 

=  Absolute  permeability. 

KCASE 

=  Index  to  indicate  if  another  case  is  to  follow. 
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KOUNT 

=  Number  of  iterations  required  for  convergence. 

L 

=  Index. 

LEND 

=  Retains  L  for  later  indexing. 

LFPND 

=  GO  TO  index,  first  sets  index  of  saturation 
distribution  to  end  and  then  reduces  the 
end  by  NSTEPS.  If  1,  set  to  end;  if  2, 
reduce  by  NSTEPS. 

M 

=  Index. 

N 

=  Index. 

NCELLS 

=  Number  of  cells  into  which  a  channel  is 
divided . 

NDELL 

=  Integer  part  of  FRCELL  (  ) . 

NDIVS 

=  Number  of  divisions  (cell  length)  into  which 
piston  is  to  be  divided. 

NDL 

=  Used  to  read  in  only  once  the  points  for  the 
permeability  lookup  curve. 

NDS 

=  Integer  for  decimal  DS;  number  of  divisions 
of  saturation  range  between  irreducible 
oil  and  breakthrough. 

NDSM1 

=  NDS  -  1. 

NDSPl 

=  NDS  +  1. 

NDVSMl 

=  NDIVS  -■  1. 

NEND 

=  Ending  index  for  number  of  units  used  in  AVPERM 
subroutine . 

NENDM1 

=  NEND  -  1,  used  for  summing  PAW's. 

NIN 

=  Integer  equivalent  to  decimal  DELN,  or  num¬ 
ber  of  saturation  increments  between,  SWI'T*7,. 
and  SWIO  for  determining  breakthrough. 

NINM1 

=  NIN  -  1. 

NINPl 

=  NIN  +  1. 

NN 

=  Index. 
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NNQ 

NSEND 

NSPCE 

NSWBR 

NT 

NTMl 
NTTP 
NTTP1 
OILRAT  ( 
OILRCV  ( 
PA  (  ) 

PABR  (  ) 

PABRL  (  ) 

PABRO  (  ) 

PABRW  (  ) 

PALE 
PAM  (  ) 


-  Index  to  determine  how  many  channels  have 
been  invaded. 

=  Number  of  divisions  in  f'grid  for  complete 
cells  invaded. 

=  Number  of  indices  between  printing  of  satu¬ 
ration  distribution  values. 

=  Index  routing  for  loose,  tight,  or  normal 
piston.  If  1,  use  Welge's  equation;  if  2, 
use  trapezoidal  equation. 

=  Index  for  finding  ends  and  beginning  of  the 
units  in  the  subroutine,  AVPERM. 

=  NT  -  1. 

=  NT  retained  before  advancincr  again  on  NT. 

=  NTTP  +  1. 

=  Oil  rate  from  Ith  channel  at  Nth  time. 

=  Oil  recovered  from  a  channel. 

=  Average  of  two  reciprocals  of  permeability; 
used  in  AVPERM. 

=  Mean  permeability  determined  from  mean  reci¬ 
procal  permeability  in  AVPERM  subroutine. 

=  Reciprocal  of  average  permeability  used  in 
AVPERM. 

=  Mean  permeability  to  oil  in  a  cell,  krQ 

mean . 

=  Mean  permeability  to  water  in  a  cell,  k 

mean . 


=  Pore  volume  of  a  channel,  V  . 

P 

=  Average  reciprocal  permeability  in  one  of 
the  units  into  which  reciprocal  permeability 
of  the  piston  is  divided. 

=  Average  reciprocal  permeability  to  oil  of 
a  unit. 


PAO 
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PAW  (  ) 

=  Average  reciprocal  permeability  to  water 
of  a  unit. 

PBWN  (  ) 

=  Reciprocal  of  permeability  at  cell  boundaries; 
used  for  determining  mean  permeability  at 
high  water  -  oil  ratios. 

PERMA 

=  Subroutine  for  calculating  permeabilities. 

PO  (  ) 

=  Permeability  to  oil,  k 

POIW 

=  Oil  permeability  at  irreducible  water  saturation 
k 

r°TT7 

IW. 

PRDCAL 

=  Calculated  pressure  drop  between  the  injection 
and  production  well. 

PRDRP  (  ) 

=  Total  pressure  drop  at  each  step. 

PRSS 

=  Pressure  drop  across  entire  channel  times 
the  absolute  permeability;  from  input  to 
output  well,  APxKa. 

PVINJ 

=  Pore  volumes  injected  in  Ith  channel  at  Nth 
time . 

PW  (  ) 

=  Permeability  to  water,  k  . 

3TW 

PWC 

=  Permeability  to  water  at  connate  water  satur¬ 
ation,  k 

rw  . 
c 

Q  (  ) 

=  Rate  of  flow  in  a  channel. 

QT  (  ) 

=  Total  injection  rate. 

RCELL 

=  Fractional  cell  invaded  in  Ith  channel  at 
any  step. 

RERROR 

=  Relative  error  between  the  injection  rate 
and  sum  of  flow  rates  in  all  the  channels. 

R  (  ) 

=  Resistance  to  flow  in  a  channel  at  Ith  step. 

RP  (  ) 

=  Reciprocal  of  permeability;  used  in  AVPERM 
subroutine . 

RPA  (  ) 

=  Right  part  of  a  unit. 

B-8 


RPAPT 

=  Right  proportional  part  of  a  unit  along  the 
f'  abscissa. 

RPO  (  ) 

=  Reciprocal  of  permeability  to  oil. 

RPW  (  ) 

=  Reciprocal  of  permeability  to  water. 

RSPWM 

=  Retain  value  of  the  past  sum  of  divisions. 

SAT 

=  Difference  between  breakthrough  and  connate; 
water  saturation. 

SAV  (  ) 

=  Average  water  saturation,  S  . 

SGMNT  (  ) 

=  Segment,  length  of  intercept  on  the  S 
abscissa  made  bv  tangent  f',  used 
to  determine  S 

w,  . 
br 

SGMPRL 

=  Value  of  SGMNT  for  previous  division,  used 
for  interpolating. 

SMUNIT  (  )  : 

=  Cumulative  sum  of  widths  of  divisions  =  !  uts . 

SPWNDM 

=  Same  as  WDIV,  used  to  retain  value  of  WDIV 
when  WDIV  is  changed. 

SR 

=  Sum  of  resistances,!  R. 

STIRLG 

=  Stirling  subroutine  name. 

SUMDIV 

=  Cumulative  sum  of  widths  of  cell,  or  distance 
from  input  well. 

SUMG 

=  Summing  terms  for  G'S  of  invaded  cells  in 
a  channel,!  G. 

SW  (  ) 

=  Water  saturation,  S  . 

w 

SWB 

=  Water  saturation,  S 

w, 

br 

SWBR 

=  Average  water  saturation  at  breakthrough, 

S 

wbr ' 

SWC 

=  Connate  water  saturation,  S  . 

w 

c 
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SWIO 

SWIW 

SWW  (  ) 

T  (  ) 
TIME  (  ) 

U 

UNITS 

UO 

UPA 

UPB 

UW 

UWUO 

VINAL 
VO  (  ) 

VOO 

WAT RAT  ( 

WDIV 
WOR  (  ) 


=  Water  saturation  at  irreducible  oil 
saturation,  S 

w 

10. 

=  Irreducible  water  saturation,  S 

w 

IW 

=  Water  saturation  for  second  indexing. 

=  Time,  t. 

=  Instantaneous  time  for  1st  channel  at  each 
step . 

=  Term  used  in  Stirling  equation. 

=  Length  of  division  of  f'  used  in  calculating 
average  areas  for  mean  permeability. 

=  Viscosity  of  oil,  Mo. 

=  Term  used  in  Newton  forward  equation,  in 
Stirling  (STIRLG)  subroutine. 

=  Term  used  in  Newton  backward  equation,  in 
Stirling  subroutine. 

=  Viscosity  of  water,  y  . 

w 

=  Viscosity  ratio  of  water  to  oil,y  /y  . 

wo 


=  Initial  volume  of  oil. 

=  Total  volume  of  oil  produced  during  an  incre¬ 
ment  . 

=  Total  volume  of  oil  produced  at  each  step 
in  channel  1. 

)  =  Water  rate  from  Ith  channel  at  Nth  time, 

qw  .  /0 

3-1/2. 

=  Width  of  a  division  in  the  AVPERM; 


=  Water  -  oil  ratio 


WORMAX 


Maximum  water 


oil  ratio. 
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WORT 10  ( 
X  (  ) 

XX 

YY 

YO  (  ) 

YW  (  ) 


=  Water  -  oil  ratio  for  Ith  channel  at  Nth  time 
=  Abscissa  in  permeability  lookup  table. 

=  Interpolated  X,  not  on  the  coincidence  values 
=  Interpolated  Y,  in  permeability  lookup. 

=  Ordinate  Y  for  oil  in  permeability  lookup. 

=  Ordinate  Y  for  water  in  permeability  lookup. 
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MAIN 


MAIN 


C  THIS  PROGRAM  CALCULATES  TWO  PHASE  FLOW  IN  ANY 

C  IRREGULARLY  BOUNDED  POROUS  MEDIUM  FOR  THE  FOLLOWING 

C  CASES 

C  -  CONSTANT  PRESSURE  DROP  BETWEEN  THE  INJECTION  AND 

C  THE  PRODUCTI  ON  WELL 

C  -  CONSTANT  INJECTION  RATE 

C  -  VARYING  INJECTION  RATE  AND  PRESSURE  DROP 

C 


COMMON  SW ( 5 10 ) *  SWW (510)  >PW ( 510 )  *PO ( 510  )  *RPW ( 510 )  ♦ 

1 RPO ( 5 10 ) »  PAW ( 510 )  ♦  PAO ( 510 )  *PABRW (510)  »PABRO( 510 )  ♦ 

IF ( 510  )  »  FP ( 5 1 0  )  *  F  P  P  (  5  1 0  )  *G ( 4  »4C )  »Q(4*500)*TIME(50Q)  » 

1 R ( 4  *  500 )  * FRCELL (4)>VCELL(4),FCELL(4)*QT(500)* 

1PRDRP (500 ) 

COMMON  SAT*NCELLS*UO*UW*POlW*UNlTS»NEND* IPAVPM* 
1NDL»K0UNT  »KCASE ♦ COMTIMtMNQ 
DIMENSION  PALE ( 4 ) »WOR (4*500)  ,SR(500)  *OILRCV ( 500 )  » 

1 WORT  10 ( 500  )  *PVINJ ( 5  00 )  »VO( 5  00 )  *VOO( 500 )  » WATRAT ( 500 ) * 
101 LRAT ( 500 ) *T ( 200 ) 

REAL  KABS 

C  INITIATION  AND  DETERM  I NA I  ON  OF  BREAKTHROUGH 

C 

NN  =  0 
SGMNT  =  0.0 
PR  SS  =  100.0 
KARS  =  11.1 
MNQ  =  0 
NDL  =  0 
WR  I  TE ( 6  *2009) 

2009  FORMAT (‘O1*'  PERFORMANCE  IN  PRIMARY  RECOVERY  ') 
READ ( 5 ♦ 1001 ) POIW  *UO»UW* 1ST  »KCASE 

1001  FORMAT ( 3E12. 5 ♦ 16  *  16 ) 

READ ( 5  *  1002 ) SW I  0  * SW I W » S WC >PWC ♦ ( PALE ( I  )  » I  =  1  *  I  ST  ) 

1002  FORMAT (4E12. 5 »4F6.0 ) 

READ (5»10C3)DELN»DS»NQ BEG »NSTLPS*NSPCE»NCELLS*IPAVPM 
10  0  3  FORMAT ( 2E 12 • 5 *5  I  5 ) 

DO  3  I  =  1  *  1ST 

READ ( 5  *  1004 )  ( G ( I  *N )  »N  =  1*NCELLS) 

1004  FORMAT (6E12. 5 ) 

3  CONTINUE 
UWUO  =  uw/uo 

DEL  =  (SWIO  -  SWIWJ/DELN 
NIN  =  DELN 
NINP1  =  NIN  +  1 
NINM1  =  NIN  -  1 
SW(1)  =  swiw 


;  i 


Y,  I  .  -ciA-v  OvVT  831  JUOJO  -■  aiHT  3 

I  J  J  -  JHT  •-  ■!.  I  d'  3  ■  Y  J  Jl  ,  1 1 

I  I  F  3  L-  1  M  v.><j  Uc.c  -  Tl  ATcSKO  - 

n  DU'.  i  3wj 

■'I  TD3u.  I  Ti  AT6 /  CO  -  3 

IT  A  I  I)L/I  C/IYflAV  - 


•  (  I  XC  )  „'HM«  (  id)  .«!  1  )  t  (  Old)  (  Id),-  ';»/  "CO 

*(olc)  <  (  -  X  d  i  . ,  «  (  '.  X  .  .  ic;»(  Id  )  i  AC,  (  0X5  )  09H I 

*  i  «  (  ■  *  •  ■  •  .  •'* .  *  (  .  *  (  «  (  J 

•(  -d  )  TO.  (  A  j  Jj  3  3,  (  A  )  Jj  ;  ;v.  (  A  )  3*(.  Jd*A)Hl 

( coe) q  'aqqi 

•  H  V  ‘  \  H  I  *  Q  IA  :  »2  I,  u  ♦ i  )•  »  U*  )l  IUJJ13.  f  TA2  ‘Ovy  O 

♦  i  l  )  ■  »  843  ♦  T  :UC  N,  jn/i 

•  tv  d  )  J 1 0 1  t  . j  .  ♦  (  .  *  •  j  )  w  •  t 

*  t  vU  A  .V *  (  1 : c  )U  V  *  l  .  0  c  )  v  ♦  (  <  d  )  U  I  V c  ♦  l  C  .  d  )  C  I  T P. 0  a  I 

(  s) '• t ooe ) TAqjiox 

J  '  3  1 

I  ‘  I  3T3C  \  .',01  T  A I  I  I 

C  = 

-.0  =  TWDZ 
c.oo x  =  zzpp 

I  •  X  X  =  2cfA> 

0  =  J(M 
( f;:  O.  d)3TI'  •: 

•  *».'•)  T  AMfl03  90 OS 
•  ■ 3 > »  T  I  *  a'  *  '>u*  :I0  •  (  XOCX  •  c  )  \3$ 

(NM*  . si 3£ )  TAwfloq  i  ;oi 

19  1  :  ‘  o  AU  )  t .  q  o . •  8  «  .1.8*  I  a  t  s o  cx  •  a  J  0 A 3 : 

(  .  JA*  C  .  SI  =♦* )  TA'  ■  q  v  £ 
1  *  -  •  *  ‘  •  *  -  *  *  ) 

(  c  I  ci  d  .  SI  3^  )  TA"  !  £0 01 
T 8 1  *  I  =  I  f  CO 
(  6  J  J  3D.'  <  I  =  1  ♦  U!  *  I  )  0  )  (  A0  0 1 «  8  )  C  A  3  ‘ 

ta.SI3d)T4V*i03  A0 ox 
3U:  IT  'CD  C 

•  l\  j  =  CU.CJ 

•  i  J 3 <  \  (  .'» I  :  -  I  ,vc  )  =  ;  3' 

J3<  -  I  ' 

i  +  ,  I  .  *  19/ I  ■' 

I  ~  , .  I , .  =  X  v  I  ■ 

I  8  =  (  X  i  .• 


n  n  n 


APPENDIX 


B 


12 


MAIN  ...  (CONT'D) 


SW(NINPl)  =  SWIO 

DO  20  N  =  2  » N  I N 

SW  (  N )  =  DEL  +  SW(N-l) 

CALL  PERMA ( N  tSW »  PO  *PW »NDL ) 

20  F(N)  =  1.0/(PO(N)/PW(N)*UWUO  +  1.0) 

F ( 1  )  =  0.0 

F(NINPl)  =  1.0 

DO  40  N  =  liNINMl 

K  =  N  I N  P  1  +  1  -  N 

SGMRPL  =  SGMNT 

FP ( K )  =  (F(K)  -  F ( K-l ) ) /DEL 

FA  =  (  F  (  K  )  +  F ( K-l  )  ) /2 • 0 

SGMNT  =  FA/FP(K) 

IF ( SGMNT- ( SW ( K) - ( SW I W+DEL /2 • 0 ) ) ) 45 » 45 >40 
40  CONTINUE 

WR I TE (6  *4002  )  K 

4002  FORMAT  (  'O’  »  '  LOOSE  PISTON  K  =  '*16) 

GO  TO  990 

43  I F ( K  -  NINP 1)50*48* 48 
48  WRITE (6 *4802 ) K 

4802  FORMAT (  *  0  *  *  *  TIGHT  PISTON  K  =  '*16) 

GO  TO  990 

50  FACTOR  =  ( SW ( K ) -SW I W— DEL / 2 • 0-SGMNT ) / ( SW ( K ) -SW I W— 
1 DE  L/2 • O-SGMNT+SGMRP  L- ( SW ( K+l ) -DEL/2 . 0-SW I W )  ) 

FPB  =  (FP(K+1)-FP(K) )*FACTOR+FP(K) 

SWB  =  (SW(K+1 )-SW(K) >#FACTOR+SW(K)-DEL/2.0 

N  =  505 

SW ( N )  =  SWB 

CALL  PERMA ( N  »  SW  »  PQ  *  PW  *  NDL ) 

FR  =  1.0/(PO(N)/PW(N)*UWUO  +  1.0) 

SWBR  =  (1.0  -  FB  )  / FPB  +  SWB 

WRITE (6  *5005 ) FPB  * FB * SWB ♦ SWBR » UWUO * K 

5005  FORMAT (' 1 '»* SLOPE  =  *  *  E 1 2  •  5  *  1  F  AT  TANGENCY  =  '» 
1E12.5 »' SATURATION  AT  TANGENCY  =  »,F10.5/'  EST 
1SATURAT ION  AT  BREAKTHROUGH  =  • » F10 • 5  * • UW/UO  =  '* 

1E14.5 » 'NORMAL  PISTON  K  =  '*14//) 

DETERMINATION  OF  DISTRIBUTION  OF  SATURATION 
OF  OIL  AND  WATER  RELATIVE  PERMEABILITIES 


DEL2  =  (SWIO  -  SWB)/DS 
NDS  =  DS 
NDSP1  =  NDS  +  1 
NDSM1  =  NDS  -  1 
SW(1)  =  SWIO 
SW(NDSPl)  =  SWB 
DO  110  N  =  1 » NDS 

110  SW ( N+ 1 )  =  SW ( N )  -  DEL2 
DO  120  N  =  1  * N D S P 1 


AND 


si  a  -  x  i  g. 

(  <_  *  '  \  j  )  ■  •  •  /  I A 


I  =  (  i  c  I  )  2 

II.  *  S  =  '/  C  OC 
(  I -.  >  )  ,v2  +  J3Q  =  ( .  )  W2 
i  J  Ji  *  3 .  C  -  ■*  d./)  AV.H3C  J  J  A  3 

( v.x  +  ouwu#(  )  q\u  )oq)\o.i  =  (/)q  cs 

0.0  =  (1)3 
o.i  =  (iqniiM3 
i  .  I/.tl  =  ,  c A  OQ 
vi  -  i  +  iq/i/  =  > 

T/.KC2  =  033/02 
J3  J  \  (  (  i->  )  -  0)3)  =  (  v  )  33 

O.S\ ( ( I-  ) -  +  ( X ) 3 )  =  A 3 
( > ) 33 \ A3  =  T//02 
( (  ,.s\jbc+  i  a)-m-A’c)- T;  vaam 

3UHIT  <00  Of* 

^  (  S 00+^*  d)  3TISW 

( cs i •  •  =  x  gt a  i  q  3ac  j  soMTA-qc^  soda 

099  OT  00 

A.8^*ce ( iq/K»  -  >)  =  i  a  a 
X  (  S08  A  •  8  )  3  T  I .  9A 

( d i ♦ '  -  uTar  t  it  •  *  * c 1 ) t a/303  see  a 

C99  OT  GO 

“  I  “  )  \  t  -  •  .  \ . J  -  .  -  )  W£ )  • 

((•i  -  .s\j30-(  i+  * )  ..a  )-jq>iv,,  (.f-T  :  Oc-c#s\j3cx 

( > ) 3H  +  * 0TDA3* ( ( > ) 33- ( 1+ » ) 33 )  =  333 
0  •  'A  J3G-  (  *  )  ..a  +  >  .  T0A3*  (  O  )  ,'8-  (  1  +  3  )  wa  )  =  3  A  a 

cue  =  /, 

hw a  =  (/)wa 

(JO  f  v  q«U  .  V  if  ,  )  A?*  -I3r  J  JAO 
(O.I  +  )U\'IA  (  /  )  <\'C\  (  /I  )  r,q  )  \0.  I  -  H3 

«va  +  flqq \ (  3  -  o.i)  =  3  awe 
f  u  o  *  a  *  l  1 . .  a  <  q  3 *  ^cT(aoca«d)3Ti  ; 

*  1  1 #e.si3« •  *  •  •  i  •  i  qo •  ccod 

i  cii  =  Y  0  ;•  '■)  r) ,  /  T  TA  01  TAHUTA2  •  •  2  •  S 131 

♦  '  *  OlArtU  •  1.  ’  -  H  01  .X 

(  \  \  A  I  «  1  =  >  TCIO  JA'FIC.  ' t  2  •  A 1 3 1 

^  '  '  •  I  TAV'JT/  a  .  GIT  J.\J  I  S3 T a  I : j  30  /Cl  TA/'lV)f!?T3Q 

Cil  i  I  J  I  :/>jV:3jq  JVITAJ3  1  X  jTA'a  C,  A  JIC 


d.\(-  a  -  c i wa )  =  s j 3G 

dG  =  20/ 

i  +  £0.  =  i^ao/i 

i  -  a  o  /  =  i-ao: 

ciwa  =  (i). a 
ewa  =  (iqac/jwa 

dC.’.I  =  a  on 

-  -  <1  o  -  •  { !♦  ,  a  o  1 1 

iHci  i./'  fl  =  /  OS  I  DG 
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MAIN  ...  (CONT’D) 

CALL  PERMA ( N * SW * PO > PW * NDL ) 

120  F ( N )  =  1 • 0/ ( PO ( N ) /PW ( N ) *UWUO  +  1.0) 

FP (  1  )  =  0.0 
FP(NDSPl)  =  FP3 
DO  130  N  =  2  » NDS 

130  FP { N )  =  (F(N-l)  -  F(N+1 ) ) /2.0/DEL2 
DEL3  =  FPB/DS 
FPP(l)  =  0.0 
FPP(NDSPl)  =  FPB 
DO  140  N  =  liNDSMl 
140  FPP ( N+l )  =  DEL3  +  FPP(N) 

SWW(l)  =  SW<1) 

SWW(NDSPl)  =  SW(NDSPl) 

DO  150  N  =  2  *  NDS 
DO  146  K  =  2 » NDSP 1 
I F ( FPP { N )  -  FP ( K >  )  148 » 144 , 146 
144  SWW(N)  =  SWU) 

GO  TO  150 
146  CONTINUE 

148  SWW(N)  =  SW ( K-l ) + ( SW ( K ) -SW ( K-l ) ) / ( FP ( K ) “FP ( K-l ) ) 

1* ( FPP ( N  )  -FP ( K-l ) ) 

150  CONTINUE 

DO  180  N  =  1 » NDSP 1 
FP(N)  =  FPP(N) 

S W  (  N  )  =  SWW(.N) 

CALL  PERMA ( N  *  SW  »  PO  *  PW  *  NDL ) 

180  F ( N )  =  1*0/ ( P  0 ( N ) /PW(N)*UWUO  +  1.0) 

190  WRITE (6 »1902 ) 

1902  FORMAT (  '  0  *  » 13X , ' SW '  ,  17X ♦ ' PW '  ,  18X  ,  ' PO '  ,  18X  ♦  ' F •  ,  16X  ,  ' FP ' 
1//  ) 

WR I TE (6  ♦  1904 )  ( SW ( N) »PW( N )  *PO( N )  *F ( N ) »FP ( N)  »N  =  lnNDSPl 
1 »NSPCE ) 

1904  FORMAT ( 5F20. 7 ) 

COMPUTATION  OF  WATER  OIL  RATIOS*  OIL  VOLUMES* 

AND  OTHER  RE  LATED  INTERMEDIATE  DATA 


200  DO  210  N  =  1  *  NDSP 1 
RPW(N)  =  1  • 0 / P  W ( N ) 

I  F { N*  EQ • 1 ) GO  TO  210 
210  RPO(N)  =  1  •  0 /PO ( N ) 

RPO(l)  =  RPO ( 2 ) 

UNITS  =  FP (2 )  -  FP ( 1 ) 

DO  220  N  =  1 # NDS 

PAO(N)  =  (RPO(N)  +  RPO(N+l))/2. 0*UN ITS 
2  20  PAW(N)  =  (RPW(N)  +  RPW ( N+  1  )  ) / 2 • 0*UN I TS 
GO  TO  (326  *327  *325  )  *KCASE 
326  QT ( 1 )  =  0.10000 
NDL  = 


1 


(  •  '  i  *  •  )  A"'  -  -1-  J  J  A3 
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MA I N  ...  (CONT'D) 


GO  TO  323 

327  OT ( 1 )  =  0.0 
NDL  =  0 

GO  TO  323 
325  OT ( 1 )  =  0.0 
NDL  =  1 

323  DO  328  I  =  1,IST 
DO  328  J  =  1,500 
R  (  I  »  J  )  =  0.0 
WOR ( I »  J  )  =  0.0 

328  CONTINUE 

DO  329  I  =  1,500 
SR (  I  )  =  0.0 
01 LRCV ( I  )  =  0.0 
WORT  I  0 ( I )  =  0.0 
PVINJ (  I  )  =  0.0 
01 LRAT ( I )  =  0.0 
VO ( I )  =  0.0 

WATRAT (  I  )  =  0.0 
OT ( I )  =  QT ( 1 ) 

329  CONTINUE 

DO  221  I  =  1 , I  ST 
SUMG  =  0.0 
DO  222  J  =  1 »N CELLS 
SUMG  =  SUMG  +  G ( I , J ) 

222  CONTINUE 

R ( I » 1 )  =  SUMG*UO/PO I W 
221  CONTINUE 

DO  223  I  =  1  *  I S T 

223  SR (  1  )  =  SR ( 1 )  +  1.0/R (1,1) 

IF ( KCASE.EQ. 1 )GO  TO  60 

GO  TO  68 

60  PRDRP(l)  =  QT ( 1 ) / (KABS*SR ( 1 ) ) 

GO  TO  65 

68  PRDRP ( 1 )  =  PRSS/KABS 

QT ( 1 )  =  PRDRPll ) *KABS*SR( 1 > 

65  VINAL  =  0.0 

DO  224  I  =  1  ,  I S T 
Q (  I , 1 )  =  (KA8S*PRDRP(1)  )/R(  1,1) 

VCELL(I)  =  PALE ( I ) /NCELLS 

224  VINAL  =  PALE(I)  +  VINAL 
SAT  =  SWBR  -  SWC 

WR I TE ( 6  * 2 25  ) 

22  5  FORMAT (  '  1  •  , 2 X , ' T I  ME '  , 3X *  'OIL-RATE'  »3X,  'WATER-RATE'  , 

1 3  X  » 'WATER-OIL-RATIO'  *3X> 'OIL-RECOVERY'  ,3X» ' PORE- VOL 
1-INJ'  ♦ 

13X* 'PRESS-DROP' *3X* ' ITRS  REQD ', 3X ,' CELLS  INVADED  IN') 
WR  I  TE ( 6  *226) 

2  26  FORMAT (  'O'* 104X  »  ' CHNL-1 '  ,  IX  » ' CHNL-2 '  » IX  »  '  CHNL-3 '  ♦ IX 
1  * ' CHNL-4 '  ) 


!  '  T  \  , ) 


/>  I  A 


:  s  £  OT  00 
.  0  =  (I) TO  YSC 


£S>E  OT  00 
0.0  =  ( I ) TO  eS£ 
I  =  JCir; 

re i . i  =  i  fse  oa  css 
i  =  c  esf  oa 

•  ■  =  (  L  f  I  )  c 

J.C  =  (Ltl)ROW 

3U/.IT  '00  8  S  £ 

ooe*!  =  i  esc  oa 

C.C  =  (1)98 
0.0  =  ( I ) VOPJ 10 
C.C  =  (DOITHOW 
C.C  =  (  I  )UMV9 
.0  =  ( I ) TA9J 10 
C.C  =  ( I ) OV 

o.O  =  (I)TAflTAW 

( I ) TO  =  ( I ) TO 

3UI  IT  00  9S£ 
T8I*X  =  I  KS  oa 
0.0  =  0MU8 

■zj .i?Ok « I  =  l  sss  oa 

(l.d-  +  ovjz  =  rO'.'ua 

3U,  IT  00  SSS 

.  IC9 \OU*cOVU8  =  (  I  ♦  I  )  9. 

30,  ITHOO  ISS 
T8I*X  =  I  ESS  oa 

(1*1)  \0.I  +  (1)93  =  (1)98  ess 

C8  OT  0  ■>(  I  .03  •  38 A 0^  MI 

8 d  OT  00 

Hi)  -  ’  =  ( i )  qac »q  oa 

da  OT  00 
8"A>\  8^  =  (I)QCoqg 

(  X  )  V8*8^A»*  (  f)  cpQ$>q  =  (  f  )  TO 

0.0  =  JA /.IV  8d 
T8I  •  X  =  I  ass  00 
•  1 1 ) a\ ( < x )qHa«9*88A;  >  =  (x# j ) 
zj  ; 3 o n \  ( i )  jja'-  =  ( i )  jjbov 

JA  IV  +  (  I  )  3JA9  =  J A1/  IV  ASS 


■  1 3TA9-X3TAW 


*  \  «  ’  Y.  JVC 


OV  2  -  99W8  =  ^  A3 
(  <;SS#  d)  3T  I  9‘‘! 

*  '  ;  "  :  *  *  X  £  « 1 3M I T  * • X  S  t  •  X  *  )  TAvqo^l  c 

■  JVOO  Jfc-JIQ'  «X£t  *01  '  *  x 

«  X  £  •  »  I  '  •  X  £  « *  • .  •  -  : 

(  dSS*  d)  9T  I  P  •/ 

*  t  XI « ’ X-J  HD  *  *  XAQX  «  IQ*  !  T  . 

( • A-J/HO ’ • I 


n  n 
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MAIN  ...  (CONT'D) 

CALCULATION  OF  PRODUCTION  HISTORY  IN  THE  PRIMARY  PHASE 


CQIYTIM  =  0.0 
M  =  0 
ICOND  =  1 
330  M  =  M+l 

MMM  =  ICOND  -  1 
JP  =  M 
NDIVS  =  M 
NEND  =  NDSP1 

CALL  AVPERMt NDIVS * NEND* RPW* PAW »PABRW»FP» UNITS* IPAVPM) 
CALL  AVPERM(NDIVS*NEND*RPO»PAO*PABRO»FP*UNITS» IPAVPM) 
TIME(M)  =  (VCELL(1)*SAT)/Q(1»M) 

MP1  =  M  +  1 
SUMG  =  0.0 

DO  350  K  =  MP1 *NCELLS 
350  SUMG  =  SUMG  +  G (  1  » < ) 

DO  360  jj  =  1 *ND I  VS 

360  R(ltMPl)  =  R ( 1  * MP 1 )  +  G( 1 *jj) / (PABRW(jj) / U W  + 

1  PABRO ( JJ ) /UO ) 

R(1»MP1)  =  R(liMPl)  +  SUMG*UO/POIW 
CALL  FLOCHK(M*JP  » MMM  * ICOND) 

VOO(M)  =  SAT*VCELL ( 1 ) *M 
OILRAT(M)  =  0.0 
SR (MP1  )  =  0.0 
VO ( M )  =  0.0 
DO  370  J  =  1 » I  ST 

OILRAT(M)  =  OILRAT(M)  +  (Q(JtMPl)  +  Q(J»M))/2.0 
SR  (  MP  1  )  =  SR  (  MP  1  )  +  2  •  0/  (  R  (  J  *MP1  )  +  R  (  j  t  M. )  ) 

I  F ( J.EQ.  I  ST ) GO  TO  370 

VO(M)  =  VO ( M )  +  (SAT*VCELL( J+l )*FRCELL( J+l ) ) 

370  CONTINUE 

VO ( M )  =  VO ( M )  +  VOO(M) 

IF (M-2 ) 377  »378  »379 

377  COMTIM  =  COMTIM  +  { T I  ME ( M )  ) / 2 . 0 

OILRCV(M)  =  V  0 ( M ) /2.0/VINAL 

GO  TO  371 

378  COMTIM  =  COMTIM  +  (TIME(M-l)  +  TIME(M) )/2.0 
OILRCV(M)  =  OILRCV(M-l)  +  ( VO ( M ) / 2 . 0 ) / V  I NA L 
GO  TO  371 

379  COMTIM  =  COMTIM  +  (TIME(M-l)  +  TIME(M))/2.0 
OILRCV(M)  =  OILRCV(M-l)  +  (  ( VO ( M )  -  VO ( M-2 )  ) / 2 . 0  ) 

1/VINAL 

371  PVINJ(M)  =  OILRCV(M) 

GO  TO  ( 376 *342 » 340 ) »KCASE 
342  CALL  FLORTE (COMT IMi MPl *QT »NDL ) 

GO  TO  376 

340  OT(MPl)  =  OILRAT(M) 


;  —  A  I  ;  A  -  A 

I  *  T  /((.  0  )  •  •  •  I  A 

2,^q  v  •  ■  .I-.S  -|,  i  v.,j  lM  I  TD’jr  J>iq  3  /CITAJ’OJAO 


O.C  =  MITMCO 
0  =  a 
I  =  0/001 
I+v  =  m  off 
I  -  /CO  I  =  MW 

V  =  qc 

v  =  2  V 1 0 '/ 

IS2C/  =  0/3/ 

(  AS  I #  2 T  •  •  »  3Z  *  2  V I  )  (  1 0 

(  S  V  A  i  I  »  TliUfSSf  A  H  A"-i  *  0  Ac  *  *  0  ' Ti  »  2  V  I  0/  )  MS  3  S  V  A  JJAO 

(■•»!)  \  (  TA2*<  I )  JJ30V)  =  (’  )3MIT 

I  +  y  =  xqw 

0«U  =  0VU2 

2j.j -jo  *ism  =  >  car  oo 

(  *  I  )  -  +  0MU2  =  "MU3  oe^ 

avi  .1  =  a  cae  oo 

+  -1 ,  \  (  L'„»  I  )  2  +  (IS  *  I  )  F  =  (IS/»I)S  Od  £ 

( 0U\ ( ll )  (  CHAS  I 

>  ioq\c  •  )' 02  +  ( i q  •i)‘-  =  nqvtijfl 

(  :oi <  <  )>:hooj3  jjao 

•  • ( I > JJ30V*TA2  =  (  ) nov 

;  .C  =  (M)TAflJIO 
o.j  =  (iqv)9£ 

C.O  =  (V)0 V 


T  2 1  *  I  =  L  CYf  00 

^  •  S  x  (  I'OU  v)  +  (IS  *c)l  )  +  ( ’  )  T  A  1 J  I  f .  =  (  C  )  T  AS  J  1 0 
(  (  «  L  )  >  +  (IS  «  L ) S ) \0 • S  +  (1 SV  }  H2  =  (  I S  -  ' )  S2 

0  \  £  OT  CO ( T  2  I .03.2 ) 31 

(  (  i  - u  )  2 _  .  -  0 '■  3  (  I  +  o  )  .  J -3 0 V * T A 2  :  *♦*  (  )  J V  =  ( ,# )  cv 


2'J  IT  oo  OY £ 

(  ■' )  oov  +  r  )cv  =  (f')ov 

v  Y £  f  8  Y £  «  Y  V  £  (  S-M)  31 

0  •  S \ (  (  ) 3 •  I T )  +  ITcro  =  IT MOO  Y Y £ 

J  A .  I  V \ J • S \ (  ■  ■ ) C V  =  (M)VOSJIO 

I Y  £  OT  00 

•  v'!(  O  I"  f  (i-  )  3  IT)  +  IT  00  =  IT^OO  8  Y  £ 

J  .  !  O  U  •  S  \  l  ’  )  O  +  (  I  -  )  .O  110  =  (  v  )  VO^1  1  I  0 

XYfc  OT  00 

C  •  A  {  (  '  )  3'-l  T  +  (I  -  )  IT)  +  IT  CO  0  =  1 1  Tm  00  9Y£ 

(  j. c\ ( ( :  -  )  c  /  -  (  )'jvn  +  (i-oyof  iic  =  (^)vosjio 

JA/IVM 

(  ) VO  J 10  =  { M u/ I  VS  I  YE 

32^0>.  {  OAc  *  S*i£»  > Y t  )  OT  00 
( J  0 C «  TO  *  X  Sc • M I T  MO  0 ) 3  T  SO J  3  JAO  Sa£ 

dY£  OT  00 

(  C  )  T  A J  I  C  =  (iqv)TO  OAE 


n  n 
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MAIN  ...  (CONT'D) 

376  PRDRP(MPl)  =  QT ( MP1 ) *1 .0/ ( SR { MP1 ) *KABS ) 

T ( M )  =  COMTIM 

PRINT  OUT  OF  RESULTS  IN  THE  PRIMARY  PHASE 


WR I TE (6  >400 ) COMT IM*OI LRAT ( M ) ♦ WAT RAT ( M )  » WORT  I  0 ( M }  ♦ 

lOILRCV(M) »PV IN J ( M ) >PRDRP(MP1)  * KOUNT  * ND I  VS  * ( FRCELL ( I )  » 
II  =  2  » 4  ) 

400  FORMAT (  « 0 '  » 1 X  * E 10 . 3  ♦  1 X  * F7 . 5 » 2 X » F7 . 5  * 6X » E 10 . 3 . 3X » 
1E12*5#3X»E12.5»3X#E11*4*5X*I4*3X*I6»1X#F6«3*1X# 

1F6 • 3  » IX  »F6 . 3  ) 

WRITE (7  »700 ) COMT IM#OILRCV (M) *PVINj(M)  ♦ WORT  10 (M) 

1 »PRDRP (MP1 ) 

7  00  FORMAT  (E12.5»4X»E12.5»4X*E12.5»4X#E12.5>4X»E11.4) 

NN  =  NN  +  1 
I F £ NN.EQ. 28 ) GO  TO  30 
GO  TO  33 
30  WRITE (6 *223) 

WRITE (6  *226) 

WRITE (6  » 2 2 7 ) 

227  FORMAT ( 1 0 ' >55X» ' PRIMARY  -  RECOVERY') 

NN  =  0 

35  IF (M.EQ.40 )GO  TO  401 
GO  TO  330 

CALCULATION  OF  PRODUCTION  HISTORY  IN  THE  SUBORDINATE 
PHASE 


401  I F ( WORT  1 0 ( M ) #  GE • 1 2 • 0 ) GO  TO  431 
M  =  M  +  1 
MM  =  M 
JP  =  M 

MMM  =  ICOND  -  1 
GO  TO  461 
431  ICOND  =  0 
M  =  M  +  1 
NSEND  =  K TP  -  2 
MMM  =  ICOND  -  1 

JK  =  I F  I  X  <  ( FLOAT ( NEND ) *FL OAT ( NCELLS )  )/ FLOAT ( NSEND )  ) 

JP  =  jK  -  M 

MM  =  M* ICOND  +  (M+jP  -  1)*IABS(MMM) 

461  NSEND  =  I  FIX (( FLOAT ( NEND ) *FLOAT ( NCELLS )) /FLOAT ( MM )  ) 

KTP  =  NSEND 

IF(N5END.LE.2)GO  TO  888 
ND  I  VS  =  NCELLS 

CALL  AVPERM (NDI VS# NSEND  »RPW*PAW  »PABRW*FP  *UNI TS  *  I PAVPM ) 
CALL  AVPERM (NDI VS • NSEND * RPO » P AO  * PABRO # FP » UNI TS» I PAVPM) 
TIME(M)  =  (VCELL ( 1 )*SAT*FLOAT ( ICOND) ) / Q ( 1 »M) +(VCELL ( 1 ) 


-  ; 1  ,  •  IA 


(  '  :>)  ...  .  i  a 

{  .  .  \  - .  i  x  q  yf  =  '  ■  • 

ivo:  -  r  jt 

V  A  I  5  I  T-  it  ■■■  “J,v  Tl  T  I'.-  9 


«(  )^IT.  •(  )  T  A ,  T  A  t<  ),  !  J  1 0  *  v'  I  T  j:  (  0A««)3"Ir 

*  (  I  )  J  J3D  -< )  *8  „■  I  -o.  *  T.  J0>*  (  13’  )  (  •  )  o,  I  VS»  (  •  )  VDP  J  I  01 

( A.S  =  II 

«>  *  .  .  I  3«  •  -  1 «  •  . '  - i*  X  t  £ •  01  t  XX«  *  C  *  )  T 

•  XX*  . .  *  I*  I*X  •  .  » X  «A«XX3*  *8.SX3«X€«€» 

(£.83* YX . r .831 

(  )0IlHv  ■#,{  )t:ilvq*(  )VD/J1  ♦  I  . ;  3  (  0  C  \  *  T  )  3  T  I  '  ■! 

( I 3v }  qciQcc .  x 

(  a  •  1 1  *  .  *  •  i  •  *  .  . .  .  «  ■  •  : .  I  *  •  .  X 

1  +■  l>  .  * 

£  OT  00  (  8S* 2./  '  )  31 
ce  oT  oo 
( ciSSt  d)  BTIF'.V  0£ 
(as^«d)iTinw 
(VSS*0)3TIflV; 

(  'Y;  d  v  03-  -  Yr  \  .  J  '  ,  X  £  d  ,  »  C  1  )  T  A'  $01  X  S 

0  =  V  ‘ 

I0A  OT  00(  0A.CH  Ml  8E 

T  A/i  1 0  ■'  J  .HI  .1  Y,'.,-Tai  i  .OlTDJCO/-  3C  >  . ,  I  T  A  JUO  J  O 

38AH3 


.  £A  •:  I  01  •  SX  .  if  •(•):.  I  T  )  3  I  1 0  A 

X  +  K  =  U 

M  =  ^  L 
I  -  0/001  =  M'lM 
IdA  0T  00 
0  =  OZODI  X  £ A 
X  +  v  =  V 
S  -  3T>.  =  C  38/ 

i  -  a/ oo  i  =  m  "■ 

;(Cj,  ,6  )  ■  AOJ  3  \  (  (  8J_. .O*  )  TAQJ3#  (  a/  j.,  )  T  AGJ3  )  )  XI  3  I  =  >w 

- 

(  )  £  H  A  i  *  (  X  -  )  +  Q 1  0  0  I  *  f  =  ' 

(  (  vv  )  TAUJ  3\  (  (  cJJBO:  )  T  A  :J3+  (0.3  )TAQJ3))XI3I  =  0/381 i  X0A 

CJZ38k  =  3  T  > 

OT  O  (  ; .  -i  J  .  10.  )  31 

8  j  j  o.  =  avion 

l  HVmM  I  «  r  I,  f  A  H  »  H»  3  p  *  *2Via/)r  33  iVA  JJA0 

*  *  •  -  ■  *  ' : *  i  •  •  ... 

(  I  )  J  J  3  3  V  )  +  (  *X)  \((  - I  )  T  0  J  3  *  T  A8*  <  X ) J  J3 2V )  =  (vj^j; 


n  n  n  o  o 
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MAIN  ...  (CONT'D) 


1*SAT*  JP  * 

1  FLOAT ( I ABS (MMM )  )  > / Q ( 1  * M ) 

MP1  =  M  +  1 

DO  402  JJ  =  1  *  ND  I  VS 

402  R ( 1  *  M  P 1  )  =  R { 1  ♦  MP  1  )  +  G( 1 »JJ) /(PABRW(jJ) / U W  + 

1  PABRO ( Jj ) /UO ) 

AVGSAT  =  SW(NSEND)  +  ((1.0  -  F ( NSEND )) /FP ( NSEND ) ) 
WOR(liM)  =  PW (NSEND )/UWUO/PO( NSEND) 

VOO(M)  =  (AVGSAT  -  SWBR ) *P ALE ( 1 )  +  VOO(NCELLS) 

CALL  FLOCHK(M*jP»  MMM  * ICOND) 

NNQ  =  1 
VO ( M )  =  0.0 
WATRAT(M)  =  0.0 
OILRAT(M)  =  0.0 
DO  409  KL  =  2*4 

I F ( (FRCELL(KL)  -  NCEL LS ) . GE . 0 . 0 ) GO  TO  411 
GO  TO  409 

411  NNQ  =  KL 
FRCELL(KL)  =  40.0 

409  CONTINUE 

IF (NNQ.EQ. 1 )GO  TO  417 
DO  412  N  =  2* NNQ 

NSEND  =  I  F  I  X (  (FLOAT(NEND)*FLOAT(N CELLS)  ) /FCELL ( N )  ) 
AVGSAT  =  SW( NSEND)  +  (1.0  -  F ( NSEND ))/ FP ( NSEND ) 

WOR ( N  »M )  =  PW( NSEND )/UWUO/PO( NSEND) 

VO  (M)  =  VO ( M )  +  (AVGSAT  -  SWBR ) *PALE ( N )  +  SAT *PALE ( N ) 

412  CONTINUE 

I F ( NNQ • EQ • 4 ) GO  TO  450 

417  NNQP1  =  NNQ  +  1 

DO  418  K  =  NNQP1 » 1ST 

OILRAT(M)  =  OILRAT(M)  +  ((Q(K»MP1)  +  Q(K»M))/2.Q) 

VO(M)  =  V  C ( M )  +  (SAT*VCELL(K)*FCELL(K)  ) 

418  CONTINUE 

450  VO ( M )  =  VO ( M )  +  VOO(M) 

OILRCV(M)  =  OILRCV(M-I)  +  (  <  VO ( M )  -  VO ( M-2 )  )  /  2  •  0  ) 
1/VINAL 

DO  421  KK  =  1 » NNQ 
QAVG  =  ( Q ( KK  *  MP 1 )  +  Q(KK*M))/2.0 
AVGWOR  =  (  WOR ( KK  » M )  +  WOR ( KK  *M-1 ) ) / 2 • 0 
OILRAT(M)  =  OILRAT(M)  +  (  QAVG/(  1.0  +  AVGWOR)) 
WATRAT(M)  =  WATRAT(M)  +  (  QAVG*AVGWOR  )/(  1.0  +  AVGWOR 
1  ) 

421  CONTINUE 

COMTIM  =  COMTIM  +  (TIME(M)  +  T I  ME ( M-l )  ) /2 . 0 
GO  TO  ( 422 »440 »445 ) »KCASE 
440  CALL  FLORTE ( COMT IN >MP1 *QT *NDL ) 

GO  TO  422 

445  GT(MPl)  =  OILRAT(M)  +  WAT  RAT ( M ) 

42  2  PVINJ(M)  =  PVI.NJ(M-l)  +  (  (  ( T I  ME ( M ) *WAT  RAT ( M )  + 

IT  I  ME  (  M-l  )  *WATRAT  (  M-l  )  )  / 2  •  0  )  +  {(VO(.M)  -  V 0  (  M - 2  )  ) 


■.I  y; 


Vi 


T  .3) 


I  A 


*TAe*i 


(  *  I  )  \ (  (  (  •  ) S  1  ) TA0J9I 

.  *  =  I qv 

ev:  o/t  i  =  oc  si  a  oq 

+  ,U\(lJ  as)  \  (lo*  I  )D  +  (!■!  ♦  I  )  ■  =  (ic  **I)fl  SQA 

(  JU\ { lUCJ  *A9  I 

(  (  -  J c,  >  -i  \  (  {  .  :u  .  )  3  —  •  1  )  )  •  ■+•  (  .  ■  )\-c  -  T ASDVA 

)  .  9  \  j  0  \  (  C. .  0 .  )  . L ;  -  (  f  X  )  ^  0  W 

i  .J  J  ;  -  -  .  .  :  -  '  '  /A)  =  (  '  , 

(  A  ON  t9u«*)  >O0_H  J  J  AO 

I  =  QlM 
C.G  =  l")OV 
0.0  =  (  •  ) T  A  V  T  A  W 
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APPENDIX  -  B 
MAIN  ...  (CONT'D) 

1/2.0)  )  /VINAL 
DO  423  K;N  =  1*1  ST 

SR(MPl)  =  SR(MPl)  +  2.0/ (R(KN*MP1 )  +  R ( K  N  *  M )  ) 
423  CONTINUE 

PRDRP(MPl)  =  QT ( MP1 ) #1 • 0/ ( SR ( MP 1 ) *KABS ) 

WORT  I  0 ( M )  =  WATRAT (M) /OILRAT (M) 

PRINT  OUT  OF  RESULTS  IN  THE  SUBORDINATE  PHASE 


49  5  WRITE (6 *400 ) COMTIM* OILRAT (M)  »  W  A  T  R  A  T ( M )  *WORTIO(M)  , 
IQILRCV(M)  *  PV I N J ( M )  »PRDRP(MP1)  * KOUNT  * ND I  VS  » 

1  ( FRCELL (  I )  » I  =  2*4) 

WRITE (7*700 )COMT I M *01 LRCV(M) *PVINj(M) » WORT  10 (M) 

1 *PRDRP (MP1  ) 

ML  =  0 

TIM)  =  COMTIM 
NN  =  NN  +  1 
IFINN.EQ.28 ) GO  TO  31 
GO  TO  37 
31  WRITE(6*225) 

WRITE (6  *226  ) 

WRITE (6  *335  ) 

335  FORMAT {'  0  '  *53X ♦ ' SECONDARY  -  RECOVERY') 

NN  =  0 

37  I F ( PVINJ ( M ) • GE. 3. 24 ) GO  TO  991 
GO  TO  401 

991  GO  TO  (990*994*994)  »KCASE 
994  KCASE  =  1 
GO  TO  326 

888  WRITE (6 *889) 

889  FORMAT ('O'*  'MAX.  WATER  OIL  RATIO  UNREACHABLE' 
990  STOP 

END 
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SUBROUTINE  PERMA 


SUBROUT INE  PERMA ( N  » SW , PO  * PW ♦ NDL ) 

this  subroutine  interpolates  the  relative 

PERMEABILITIES  ACCORDING  TO  THE  SATURATION  GRID 


DIMENSION  S W { 5 1 0  )  »  P  0 ( 5 1 0 )  *  P  W  (  5  1 0  )  >  Y  0  (  1 1 0  )  » Y  W ( 11 0 )  » 
1 X ( 110)  *  D 1 ( 110 )  » D2 ( 1 10 )  »  D3 ( 1 10  )  *  D1 1 { 1 10 )  , 

1D22 ( 110  )  » D3 3 ( 110 ) 

IF(NDL) 100 » 100  »  200 
100  NDL  =  1 

READ { 5  » 1 00 1 ) JPO I  NT  »ENDPO»  ENDPW*ENDSW 

1001  FORMAT ( I  3  *  3E 12  •  5 ) 

READ ( 5  »  1002 )( YO ( j )  ,  j  =  lOtJPOINT) 

READ ( 5  *  1002 )  ( YW ( J )  »  J  =  lOtJPOINT) 

READ ( 5 » 1002 ) (  X(j)»  j  =  10»JPOINT) 

1002  FORMAT (6E12. 5 ) 

JPNTM1  =  JPOINT  -  1 

JPNTM2  =  UPOINT  -  2 

JPNTM3  =  UPOINT  -  3 

DO  110  J  =  10  *  JPNTM1 

110  D 1  (  j  )  =  YO(  j+1  )  -  YOU) 

DO  120  J  =  10 » JPNTM2 
120  D2 ( J )  =  Dl( J+1)  -  D 1 ( J ) 

DO  130  J  =  10 » JPNTM3 
130  D3 ( j )  =  D2 ( J+1 )  -  D2 ( J ) 

DO  160  J  =  10,JPNTM1 
160  Dll(j)  =  Y W ( J  + 1 )  -  YW ( J ) 

DO  170  J  =  10  f  JPNTM2 
170  D2  2 ( j  )  =  DlltJ+1)  -  Dll(J) 

DO  180  j  =  10 » JPNTM3 
180  D33  (  J  )  =  D22  (  J+ 1  )  -  D22U) 

200  IF ( ENDSW-SW (N ) 1300,300,210 
210  XX  =  SW ( N ) 

CALL  STIRLG( JPOINT, YO ♦ XX, X»D1 »D2*D3»YY) 

PO ( N )  =  YY 

CALL  STIR LG (JPOINT  *YW»XX»X,D11 #D2  2 ,D3  3 , YY ) 

PW ( N )  =  YY 
GO  TO  500 
3  00  S W ( N )  =  ENDSW 

PO ( N )  =  ENDPO 
PW(N)  =  ENDPW 
500  RETURN 
END 
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SUBROUTINE  ST  I R  LG 

SUBROUTINE  ST  I RLG ( JPO I  NT » Y  » XX  * X  * D 1 »D2 » D3  * Y Y ) 

C  THIS  SUBROUTINE  DETERMINES  WHICH  OF  THE  FOLLOWING 

C  EQUATIONS  SHOULD  BE  USED  TO  INTERPOLATE  THE  RELATIVE 

C  PERMEABILITIES  FROM  THE  RELATI  VE  PERMEABILITY  VERSUS 

C  SATURATION  DATA 

C  -  NEWTONS  FORWARD  DIFFERENCE  EQUATION 

C  -  STIRLING  EQUATION 

C  -  NEWTONS  BACKWARD  DIFFERENCE  EQUATION 

C 

DIMENSION  Y(11Q)  »  X ( 1 1 0 )  *  D  1  ( 1 1 0 )  #02(110)  *  D3 ( 1 1 0 ) 

H  =  X  (  1 1  )  -  X  (  10  ) 

DO  70  J  =  11 *jPOINT 
I  F ( X ( J )  -  XX)70*30*90 
70  CONTINUE 
GO  TO  90 
80  YY  =  Y ( J ) 

GO  TO  500 
90  MM  =  J-l 

I F ( MM  -  12)100*200*300 
100  U  =  (XX  -  X (MM) ) /H 

UPA  =  (U  -  1.0)*U/2.0 

YY  =  Y ( MM )  +  U*D1(MM)+UPA*D2(MM)+UPA*(U-2.0)*D3(MM) 
1/3.0 

GO  TO  500 
200  U  =  ( XX-X ( MM ) ) /H 

yy  =  y ( MM ) + ( D 1 ( MM-1 )+Dl (MM) ) *U/2 . 0+U**2*D2 (MM-1 ) /2.0+U 
1 * ( U##2- 1 • 0  ) 

1 ( D 3 ( MM- 2 )+D3(MM-l  )  )/12.0 
GO  TO  500 

300  IF ( MM+3-J POINT ) 200  *200 *400 
400  U  =  ( XX-X ( MM+1 ) ) /H 
UPB  =  U* ( U+l . 0 ) /2 .0 

YY  =  Y (MM+1 ) +  U*D 1 ( MM ) +UPB*D2 ( MM-1 ) +UPB* ( U+2 . 0 ) *D3 ( MM 
1-2 ) /3.0 
500  RETURN 

END 
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SUBROUTINE  AVPERM 

SUBROUT INE  AVPERM ( ND I  VS  * N END  » R P  >PA  ,  PABR » FP  »UNI TS 
1  »  I  PAVPM ) 

THIS  SUBROUTINE  IS  USED  TO  DETERMINE  THE  MEAN 
PERMEABILITY  TO  OIL  AND  WATER  WHEN  THE  NUMBER  OF  CELLS 
INVADED  IS  A  WHOLE  NUMBER 


DIMENSION  RP ( 510  )  *  PA ( 510 )  »PABRL ( 510 )  *FP ( 5 10 )  » 

IP AM (510)  *  RP  A ( 5 1 0 )  *  ALP A (510)  t PABR ( 5 10 )  ♦ RSWD I  V ( 5 10 )  , 

1 PB WN ( 510 )  »  SMUN I T (510) 

1  DNDIVS  =  ND I  VS 

WDIV  =  FP (NEND) /DNDIVS 

IF ( 2.0*UNITS  -  WDIV) 300*300*200 

CALCULATION  OF  MEAN  PERMEABILITIES  WHEN  NUMBER 
OF  DIVISION-  S  IS  LOW  -  HIGH  WATER  OIL  RATIO 


200  NDVSM1  =  ND I  VS  -  1 
RSWD I  V ( 1 )  =  WDIV 
SMUN I T ( 1 )  =  UNITS 
DO  205  N  =  2  » NEND 

205  SMUNIT(N)  =  SMUNIT(N-l)  +  UNITS 
DO  210  N  =  2  » ND  I  VS 
2  10  RSWD I  V ( N )  =  RSWDIV(N-l)  +  WDIV 
DO  260  K  =  1  * NDVSMl 
DO  240  NN  =  1  * N E N D 
N  =  NN 

IF ( SMUN I T ( N )  -  RSWDIV(K) ) 240,240,250 
240  CONTINUE 

250  IF ( N-l ) 251 ,251 ,252 

251  PBWN ( K  +  l  )  =  RP ( 1 )  -  (RP( 1 )-RP (2 )  ) *RSWD I  V ( K ) /UNI TS 
GO  TO  260 

252  PBWN ( K  +  l )  =  RP(N)-(RP(N)-RP(N+1 )  ) * ( RSWD I  V ( K ) -SMUN I T ( N 
1-1  )  ) /UNITS 

260  CONTINUE 

PBWN(l)  =  RP ( 1 ) 

PBWN ( ND I VS  +  1 )  =  RP ( NEND  +  1) 

DO  270  K  =  1 » ND I  VS 

270  PABRL(K)  =  (PBWN(K)  +  PBWN ( K+l )) /2 . 0 
GO  TO  600 

CALCULATION  OF  MEAN  PERMEABILITY  WHEN  NUMBER  OF 
DIVISIONS  -  IS  HIGH 


300  IF ( ND IVS  -  1)310,310*340 
310  J  =  1 

NENDM1  =  NEND  -  1 
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SUBROUTINE  AVPERM. • •  (CONT'D) 


PAM(J)  =  0.0 
DO  320  K  =  1  * N E N D M 1 
3  20  PAM(J)  =  PAM(J)  +  PA ( K  > 

PABRL(J)  =  P AM ( j ) /WD I V 
GO  TO  600 
340  J  =  1 

SPWNDM  =  WD I V 
DO  350  K  =  1»NEND 
IF ( F  P ( K  )  -  WD I  V ) 350  *360  » 360 
350  CONTINUE 
360  NT  =  K  -  1 

RP APT  =  WD I V  -  FP(NT) 

RP A ( j  )  =  ( RP ( NT ) -( RPAPT /UNI TS* (RP (NT ) -RP ( NT+ 1 )  )  ) / 2 . 0 ) 
1*RPAPT 
PAM(J)  =  0.0 
NT Ml  =  NT  -  1 
I F ( NTM1  -  1)380,370,370 
370  DO  375  K  =  1*NTM1 
375  PAM(j)  =  PAM(J)  +  PA ( K ) 

380  PABRL(J)  =  (PAM(j)  +  RPA(J))/WDIV 
400  DO  500  J  =  2  » ND I  VS 
PAM(J)  =  0.0 
DJ  =  J 

RSPWM  =  SPWNDM 

SPWNDM  =  WDIV*DJ 

ALPAPT  =  FP ( NT  +  1)  -  RSPWM 

ALPA(J)  =  (RP(NT+1)+{ALPAPT/UNITS*(RP(NT)“RP(NT+1) ) ) 
1/2.0) *ALPAPT 
NT  TP  =  NT 

IF ( j-NDI VS )430  *42  0  >420 
420  RPA(J)  =  0.0 

NT Ml  =  NEND  -  1 
NT  TP  1  =  NTTP  +  1 
GO  TO  450 

430  DO  440  N  =  NTTP ♦ NEND 

I  F ( FP ( N  )  -SPWNDM ) 440  >445  *  445 
440  CONTINUE 
445  NT  =  N  -  1 

RPAPT  =  SPWNDM  -  FP(NT) 

NT  TP  1  =  NTTP  + 1 
NTM1  =  NT  -  1 

RPA(j)  =  (RP(NT)-(RPAPT/UNITS)*(RP(NT)-RP(NT+1) )/2.0) 
1* RPAPT 

450  I F ( NT-NTTP  -  1)460*490*470 
460  IF ( NEND-NT+1 ) 490 *490*470 
470  DO  480  K  =  NTTP1  »NTM1 
480  PAM(J)  =  PAM(J)  +  P  A ( K ) 

490  PABRL(J)  =  ( P AM ( J ) +ALP A ( J ) +RP A ( J )  ) / WD I  V 

500  CONTINUE 
6  00  DO  610  N  =  1  *  N D I V S 
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SUBROUTINE  avperm... 

610  PABR(N)  =  1.0/PABRUN) 

I F ( 1  “  I PAVPM ) 620 » 620 » 630 
620  WRITE (6 *6200 JNDIVS* (PABR(N)  »N 
6200  FORMAT!//'  MEAN  PERMEABILITIES 
1  INVADED  '  »  '  =  '  ,  I 
1/ <  5E15.5  )  ) 

630  RETURN 
END 
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SUBROUTINE  FLQRTE 

SUBROUTINE  FLORTE ( E LT I  ME  * MP 1  * GT  * NDL ) 

THIS  SUBROUTINE  IS  USED  TO  INTERPOLATE  FOR  INJECTION 
RATE  FROM  TI-  ME  VERSUS  INJECTION  RATE  HISTORY  AND  IS 
USED  IN  A  CASE  WHEN  BOTH  INJECTION  RATE  AND  PRESSURE 
DROP  VARY 


DIMENSION  T  I  ME ( 5  00 )  »QT (300)  *  X ( 1 10 ) *  Y ( 1 10 ) *011110) 
1  *D2 (110 ) *  D3 ( 1 10 ) 

IF ( NDL ) 100  * 100  *190 
100  NDL  =  1 

READ ( 5*1001 ) JPO I  NT  » ENDT I M  * ENDQT 

1001  FORMAT (  I3.2E12.5  ) 

READ ( 5 »  1002  )  (  X(J)»  J  =  10*jPOINT) 

READ! 5  *1002  )  (  Y(j)»  J  =  10»jPOINT) 

1002  FORMAT ( 6E12. 5  ) 

JPNTM1  =  JPOINT  -  1 
JPNTM2  =  JPOINT  -  2 
JPNTM3  =  JPOINT  -  3 
DO  110  J  =  10 » JPNTM1 

110  D 1 ( J)  =  Y( J+l)  -  Y( J) 

DO  120  J  =  10 » JPNTM2 
120  D2 ( J )  =  D1 ( J+l )  “  D 1 ( J ) 

DO  130  J  =  10  *  JPNTM3 
130  D3 ( J )  =  D2 ( J+l )  “  D2 ( J ) 

190  IF ( E L T I M E  -  X ( 1 0  )  )  5  ♦  5  ♦  2 00 
5  QT ( MP 1  )  =  Y ( 10 ) 

GO  TO  500 

200  I F ( ENDT I M  -  ELT I  ME ) 300  » 300 ♦ 2 1 0 
210  XX  =  ELTIME 

CALL  ST IRLG ( JPOINT »Y  » XX  *  X » D 1  * D2 » D3 ♦ Y Y ) 

QT(MPl)  =  YY 
GO  TO  500 
300  QT(MPl)  =  ENDQT 
500  RETURN 
END 
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SUBROUTINE  FLOCHK 
SUBROUTINE  FLOCHK ( M  » JP  *MMM  * ICOND) 

this  subroutine  is  used  to  determine  flow  rates  in 

VARIOUS  CHANNE-  LS  IN  A  MANNER  SO  THAT  CROSS-FLOW 
BETWEEN  THE  CHANNELS  IS  MINIMIS-  ED 


COMMON  SW ( 5 10 ) »SWW ( 510 )  • P W ( 5 1 0 )  *PO ( 510 )  *RPW ( 510 )  » 

1RP0( 510  )  ,PAW ( 510  )  ♦ P AO ( 5 1 0 )  ♦  PABRW ( 5 10 )  * PABRO ( 510 )  t 
IF (510)  »  FP ( 5 1 0 ) *  FPP ( 5 1 0 )  »G(4»4U)  *  Q  (  4  »  5  0  0  )  »  T I M  E ( 5  0  0 )  * 

1 R  (  4  »  5  0  0  )  *FRCELL(4)  *  V  C  E  L  L ( 4  )  »  F  C  E  L  L  (  4  )  »QT  (500)  » 

1PRDRP ( 500  ) 

COMMON  SAT*NCELLS*UO*UW»POIW»UNITS*NEND» IPAVPM» 

1NDL  *KOUNT * KCASE  *  COMT I M * MNQ 
REAL  KABS 
KABS  =  11.1 
KOUNT  =  1 
I  I  I  =  0 
NS END  =  NEND 
MP 1  =  M  +  1 
DO  5  I  =  1 » 4 
5  Q  (  I  » M  P 1  )  =  Q  (  I  *M  ) 

IF ( M-2  )  10  *21 *21 
10  DO  15  I  =  1*4 
15  FCELL (  I  )  =  0.0 
21  NO  =  2 
20  DO  25  L  =  NQ  * 4 

FRCELL(L)  =  FCELL ( L )  +  (  (  ( Q ( L *M ) +Q ( L »MP1 )  ) /2 .0 ) *T I  ME ( M )  ) 
1 / ( VCELL ( L  )  * 

1SAT) 

NSEND  =  NEND 

IF ( FRCELL ( L) .GE.40.0)GO  TO  300 
GO  TO  301 

300  NSEND  =  I  FIX (( FLOAT ( NEND ) *FLOAT ( NCELLS )) /FRCELL ( L )  ) 
NDIVS  =  40 

GO  TO  19 

301  NDELL  =  I F I X ( FRCELL( L )  ) 

RCELL  =  FRCELL(L)  -  FLOAT (NDELL) 

DEV  =  FRCELL ( L ) /501 . 

IF (RCELL. LE. DEV) RCELL  =  0.0 
751  IF (FPCELL(L) *LT.1.0)GO  TO  17 
IF (RCELL. NE. 0.0 )GO  TO  30 
NDIVS  =  NDELL 
NDLP1  =  NDELL  +  1 
GO  TO  19 
17  NDIVS  =  1 

19  CALL  AVPERM ( ND I  VS  » NSEND » RPW *  PAW  » PABRW * FP  » UN  I TS ♦ I PAVPM ) 
CALL  AVPERM ( ND I  VS  * NSEND » RPO » P AO  » PABRO »  FP  » UN I TS » I PAVPM ) 
IF(FRCELL(L).LT.1.0)G0  TO  71 
IF(FRCELL(L).GE.40.0)GO  TO  303 
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SUBROUTINE  FLOCHK.  •  •  •  (CONT'D) 

R ( L.MP1 )  =0.0 
SUMG  =  0.0 

DO  100  K  =  NDLP 1  * NCELLS 
100  SUMG  =  SUMG  +  G ( L  *  K ) 

GO  TO  302 
303  SUMG  =  0.0 

R ( L ♦ M  P 1 )  =  0.0 

302  DO  110  Jj  =  1  *  ND  I  VS 

110  R ( L  »MP1 )  =  R ( L ♦ MP 1 )  +  G(L*JJ) / (PABRW(JJ) /UW  + 

1  PABROl JJ ) /UO) 

R ( L  » M P 1  )  =  R ( L  »  MP 1  )  +  SUMG  *  UO/POIW 
GO  TO  25 

30  CALL  AVPER2(NDELL*NEND*L*FRCELL*RPW*PAW*PABRW*FP*UNITS 
1 » I PAVPM ) 

CALL  AVPER2 ( NDELL * NEND * L  > FRCELL * RPO ♦ PAO » PABRO ♦ FP  *  UN  I TS 
1  *  I PAVPM ) 

71  K2  =  NDELL  +  2 
K 1  =  NDELL  +  1 
400  SUMG  =  0.0 

I F ( K2 .GT.4Q.0 )GO  TO  205 
DO  200  I  =  K2  ♦  NC E LLS 
200  SUMG  =  SUMG  +  G ( L  ♦  I  ) 

205  SUMG  =  SUMG  +  G ( L  »  < 1 ) * ( 1.0  -  RCELL ) 

410  R ( L  *  MP1  )  =0.0 
DO  210  J  =  1  *  K 1 
IF ( J. EQ.K1 )GO  TO  215 

R ( L  »  MP 1  )  =  R ( L  *  M  P 1 )  +  G ( L » J ) / ( PABRW ( J ) /UW  +  PABRO(J) 

1/UO) 

GO  TO  210 

215  R ( L  *  MP1  )  =  R(L.MP1 )+G(L.J )*RCELL/( PABRW (J) /UW  + 

1  P ABRO ( J ) /UO ) 

210  CONTINUE 

R ( L  » M P 1 )  =  R ( L  » MP 1 )  +  SUMG*UO/PO I W 

25  CONTINUE 

GO  TO  (2*2*3)  > KC ASE 

2  PROD  =  0.0 

DO  27  K  =  1*4 

QAVG  =  (  Q ( K  *M )  +  Q ( K  »  MP 1 )  ) / 2 • 0 

RAVG  =  (  R(K*MP1)  +  R(K*M))/2«0 

PROD  =  PROD  +  QAVG*RAVG 

I  F (  ( KOUNT.EQ. 10 )  . AND. ( K. EQ* 1 )  ) GO  TO  700 
27  CONTINUE 

PROD  =  PROD/4 • 0 
AVPROD  =  PROD 
GO  TO  725 

700  PROD  =  (PROD  +  AVPRODJ/2.0 
KOUNT  =  0 
III  =  III  +  1 
725  DO  32  N  =  1*4 

QAVG  =  (  Q(N»M)  +  Q(N*MP1))/  2.0 
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SUBROUTINE  FLOCHK...  (CONT'D) 

RAVG  =  (  R ( N  *  M )  +  R ( N  *  M P 1  )  )  / 2  •  0 
CHECK  =  ( QAVG*R A VG  -  PROD) /PROD 
IF ( I  I  I .GE. 1 ) GO  TO  28 
GO  TO  29 

28  WRITE (6  ,18 )CHECK 

18  FORMAT ('  0 '  »4X » '  CHECK  =  ’,E13.6) 

IPAVPM  =  1 

WR I TE ( 6 ,45  )  ( FRCELL ( L )  ,  L  =  2,4) 

45  FORMAT ('  0  '  >3 ( 2X » E13.6 )  ) 

WRITE(6»777)PROD» (Q(K»MP1 ) *K  =  1 » 4 ) , ( R ( j ♦ MP 1 ) , j  =  1,4) 
777  FORMAT ('0',1X,E13.6,4(1X»E13.6),4(1X,E13.6)) 

29  IF (ARS(CHECK) •GE.1«E-04)GO  TO  22 

32  CONTINUE 
GO  TO  280 

22  DO  33  I  =  1,4 

QAVG  =  (  Q ( I  , M )  +  Q ( I • MP 1  )  )  / 2 • 0 

RAVG  =  ( R ( I , M )  +  R ( I »MP  1  )  ) /2 • 0 

QAVG  =  QAVG  -  ((  QAVG*RAVG  -  PROD)/  RAVG  ) 

Q (  I , M P 1 )  =  2 • 0*QAVG  -  0(1, M) 

33  CONTINUE 

705  TIME(M)  =  (VCELL(l)*SAT*FLOAT(ICOND) )/ ( (Q(1»M)+ 

1Q(  1  ,MP1  )  ) /2.0  )  +  ( VCELL ( 1 ) *SAT*jP*FLOAT (  I  ARS ( MMM 
1 )))/(( Q ( 1 , MP 1 )  +  Q(1»M)  ) / 2  •  0 ) 

IF ( KCASE.EQ. 1 )GO  TO  711 
EL TIME  =  COMTIM  +  TIME(M) 

IF (MNQ.GE.l )GO  TO  710 
IF ( ELTIME.GT.O. 13893E  04)GO  TO  703 
GO  TO  710 
703  MNQ  =  MNQ  +  1 
NDL  =  0 

710  CALL  FLORTE ( ELTIME*MP1 »QT »NDL ) 

711  SUM  =  0.0 

DO  43  II  =  1,4 

43  SUM  =  SUM  +  Q ( I  I  »  MP 1 ) 

IF (KCASE.EQ. 1 )GO  TO  712 

RERROR  =  ((SUM  -  QT ( MP1 ) ) /QT ( MP 1 ) ) 

IF ( ABS ( RERROR ) . LE . 1 . E-04 ) GO  TO  707 

712  DO  44  KK  =  1,4 

Q ( KK , MP 1 )  =  (Q( KK,MP1 )*QT (MP1 ) ) /SUM 

44  CONTINUE 

I F ( KCASE • EQ. 1 ) GO  TO  707 
GO  TO  705 

707  KOUNT  =  KOUNT  +  1 

IF ( KOUNT.EQ.80 ) GO  TO  310 
GO  TO  20 

280  KOUNT  =  KOUNT  +  111*10 
GO  TO  130 
3  DO  120  K  =  1,4 

QAVG  =  ( Q ( K , M )  +  Q ( K , M  P 1 )  ) / 2 • 0 

RAVG  =  ( R ( K , M  )  +  R ( K , MP 1 )  ) /  2.0 
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SUBROUTINE  FLOCHK...  (CONT'D) 

PRDCAL  =  (QAVG*RAVG ) /KABS 

PRCHK  =  (PRDCAL  -  PRDRP ( 1 )  ) /PRDRP ( 1 ) 

IF (ABS(PRCHK) .GE . 1 . E-05 ) GO  TO  122 
120  CONTINUE 
GO  TO  130 

122  TIME(M)  =  ( VCELL ( 1 )*SAT*FLOAT ( ICOND) ) 

1/ (  (Q( 1 »M)+Q ( 1 »MP1 )  )  /2.0  )  + ( VCELL ( 1 ) *SAT*jP* 
1  FLOAT  (  I  ABS  (  MMM  )  )  )  /  (  (  0  (  1  » MP 1  )  +  Q  (  1  » i-  )  )  / 2  •  0  ) 
QAVG  =  (QAVG*PRDRP ( 1 )) /PRDCAL 
Q ( K  ♦  MP 1 )  =  ( 2  •  0*QA VG  -  Q ( K  * M )  ) 

KOUNT  =  KOUNT  +  1 
IF(KOUNT.EQ.80)GO  TO  310 
GO  TO  20 

130  DO  290  L  =  1»4 

290  FCELL(L)  =  FRCELL(L) 

310  RETURN 

END 
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SUBROUTINE  AVPER2 

SUBROUT  INF  AVPER 2 ( NDELL  * NEND >  L  > FRCELL  > RP  »  PA  » PABR  * FP 
1 » UN  I TS  > I PAVPM ) 

THIS  SUBROUTINE  IS  USED  TO  DETERMINE  THE  MEAN 
PERMEABILITY  TO  OIL  AND  WATER  IN  A  FRACTIONAL 

CELL 


DIMENSION  FRCELL ( 4  )  »RP ( 5 10  )  *PA ( 5 10 )  »PABRL ( 510 )  » 

1FP ( 510  )  »PAM( 510 )  » RP  A  ( 510 )  » ALP  A ( 5 1 0 )  » PABR ( 5 1 0 ) ♦ 

1RSWDI V ( 51 )  » PBWN (51)  » SMUN  I  T { 510 )  * 

1RSWDIV ( 51 )  » PBWN ( 5 1 )  >SMUNIT(510) 

I  END  =  IFIX (FLOAT (NEND)*(  ( FRCELL ( L ) -FLOAT ( NDELL )  ) 

1 /FRCELL ( L ) ) ) 

ND END  =  NEND  -  I  END 

WDIV  =  FP (NDEND) /FLOAT (NDELL) 

NPDEL 1  =  NDELL  +  1 

IF (2*0# UNITS  -  WDIV) 300 »300 »200 

CALCULATION  OF  MEAN  PERMEABILITIES  WHEN  NUMBER 
OF  DIVISION- 

S  IS  LOW  -  HIGH  WATER  OIL  RATIO 

200  NDLLMl  =  NDELL  -  1 
RSWDIV(l)  =  WDIV 
SMUN I T ( 1  )  =  UNITS 
DO  205  N  =  2  # NDEND 
205  SMUNIT(N)  =  SMUNIT(N-l)  +  UNITS 
DO  210  N  =  2  » NDE L L 
210  RSWDIV(N)  =  RSWDIV(N-l)  +  WDIV 
DO  260  K  =  1  ♦  NDLLMl 
DO  240  NN  =  1  *NDEND 
N  =  NN 

I F ( SMUN I T ( N )  -  RSWDIV(K) )240»240»25Q 
240  CONTINUE 

250  I  F ( N-l  ) 25  1  *251  ♦ 252 

251  PBWN ( K  +  l )  =  RP ( 1  )  -  ( RP ( 1 ) -RP (2 )  ) *RSWD I  V ( K ) /UNI TS 
GO  TO  260 

25  2  PBWN (K  +  l)  =  RP ( N ) - ( RP ( N ) -RP ( N  +  l )  ) * ( RSWD I  V ( K ) -SMUN I T ( N 
1-1 )  ) /UN  I TS 
260  CONTINUE 

PBWN (  1 )  =  RP (  1  ) 

PBWN ( NDELL  +  1 )  =  RP ( NDEND+ 1  ) 

DO  270  K  =  1 »  NDE L L 

270  PABRL(K)  =  (PBWN(K)  +  PBWN ( K+l )) /2 . 0 
275  ALOK  =  0.0 

NENDM1  =  NEND  -  1 
DO  280  K  =  NDEND *NENDM1 
280  ALOK  =  ALOK  +  PA ( K ) 

PABRL ( NPDEL1  )  =  ALOK/ ( FP ( NEND )  -  FP ( NDEND )  ) 
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SUBROUTINE  AVPER2...  (CONT'D) 


GO  TO  600 

CALCULATION  OF  MEAN  PERMEABILITY  WHEN  NUMBER  OF 
DIVISIONS  - 
IS  HIGH 

300  IF ( NDELL  -  1)310.310*340 
310  J  =  1 

NDNDM1  =  NDEND  -  1 
PAM(J)  =  0.0 
DO  320  K  =  1 » NDNDM1 
320  PAM(J)  =  PAM(J)  +  PA ( K ) 

PABRL(J)  =  PAM ( j  )  /WD I  V 
GO  TO  275 
340  J  =  1 

SPWNDM  =  WD I V 
DO  350  K  =  1, NDEND 
I F ( FP ( K  )  -  WDIV) 350.360.360 
350  CONTINUE 
360  NT  =  K  -  1 

RPAPT  =  WDIV  -  FP(NT) 

RPA(J)  =  (RP(.NT)-(RPAPT/UNITS*(RP(NT)-RP(NT+1)  )  )  /  2 . 0  ) 
1* RPAPT 
PAM(J)  =  0.0 
NTM1  =  NT  -  1 
IF ( N  T  M 1  -  1)380.370,370 
370  DO  375  K  =  l.NTMl 
375  PAM(J)  =  PAM(j)  +  PA ( K ) 

380  PABRL(J)  =  (PAM(J)  +  RPA(J))/WDIV 
400  DO  500  J  =  2, NDELL 
P AM ( J )  =  0.0 

DJ  =  J 

RSPWM  =  SPWNDM 

SPWNDM  =  W  D I V  *  D  J 

ALPAPT  =  FP ( NT  +  1)  -  RSPWM 

ALPA(J)  =  (RP(NT+1)+(ALPAPT/UNITS*(RP(NT)-RP(NT+1) ) ) 
1/2.0) *ALP APT 
IF ( J-NDELL ) 43 0,42  0 ,420 
420  RPA(J)  =  0.0 

NTM1  =  NDEND  -  1 
NT  TP  1  =  NTTP  +  1 
GO  TO  450 

430  DO  440  N  =  NTTP, NDEND 

I F ( FP ( N  )  -  SPWNDM ) 440 .44  5 , 445 
440  CONTINUE 
445  NT  =  N  -  1 

RPAPT  =  SPWNDM  -  FP ( NT ) 

NTTP 1  =  NTTP+1 
NT Ml  =  NT  -  1 

RPA(j)  =  (RP(NT)-(RPAPT/UNITS)*(RP(NT)-RP(NT+1) ) /2.Q) 
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SUBROUTINE  AVPER2...  (CONT'D) 


1*RPAPT 

450  I  F ( NT-NTTP  -  1  ) 460  * 490  > 47 0 
460  IF (NDEND-NT+1 )490»490»47Q 
470  DO  480  K  =  NTTP1  ♦  N T M 1 
480  PAM(J)  =  PAM(j)  +  PA ( K ) 

490  PABRL(U)  =  (PAM(U)+ALPA( J)+RPA(J) J/WDIV 
5 CO  CONTINUE 
GO  TO  275 

600  DO  610  N  =  1  » NPDEL 1 
610  PABR(N)  =  1.0/PABRLIN) 

IFtl.O  -  I P  A  V  P  M ) 620  »620»630 
620  WRITE (6  >6200 ) N  D  E  L  L  * ( P  A  B  R ( N) »N  =  ltNDELL) 

6200  FORMAT (// 'MEAN  PERMEABILITIES'  *  5  X  » 1 C  OM  P .  CELLS 
1INV.  =  '  » I6*//(5E15.5)  ) 

WRITE (6  »6  300) PABR ( NPDEL  1  ) 

6300  FORMAT!//  '  MEAN  PERMEABILITY  IN  FRACTIONAL  CELL  =  ' 

1  *  E 1 5 . 5  ) 

WRITE (6  >6400 ) L 

6400  FORMAT (' 0 '  »4X  ♦  '  CHANNEL  NUMBER  =  '>14) 

630  RETURN 
END 
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SUBROUTINE  CAPLOT 

SUBROUTINE  CAPLOT ( X 3  * Y3  * X4  » Y4  * X6 , Y6  » X7 , Y7  » KCPD  » KCFR ) 
THIS  SUBROUTINE  IS  USED  TO  PLOT  THE  RESULTS 


DIMENSION  X(32)  *  Y (32)  »  X 1 ( 32  )  *  Y 1  (32)  »X2(32)  »  Y2 ( 32  )  » 
1X3(200)  *  Y  3  (200)  »X4(200)  »Y4(20G)  »X5 ( 32  )  *  Y  5 ( 3  2 )  » 
1X6(200)  *Y6(200)  »X7(200)  »Y7(200)  .WORK ( 102  4 ) 

CALL  PLOTS (WORK ( 1 ) *4096 ) 

READ ( 5  » 10 )  (XI  (I)  »  I  =  1,30) 

READ ( 5 ,10 )  ( Y1  (  I  )  ,  I  =  1,30) 

RE AD ( 5 • 10 )  ( X2  (  I  )  *  I  =  1,30) 

RE AD ( 5  » 10 )  ( Y2 ( I )  *  I  =  1,30) 

READ ( 5  , 10  )  (X5 ( I )  » I  =  1*30) 

READ ( 5,10 )  (Y5 (  I  )  »  I  =  1,30) 

10  FORMAT ( 6E12. 5 ) 


DO  20 

J  =  1,30 

X  5  (  J  ) 

=  X5 ( j ) /10.0 

20  CONTINUE 

CALL 

PLOT (  0.0,  0.0, 

3  ) 

CALL 

PLOT (  0.0,12.0, 

2  ) 

CALL 

PLOT ( 18.0,12.0* 

2  ) 

CALL 

PLOT ( 18.0,  0.0, 

2  ) 

CALL 

PLOT (0.0  ,0.0 ,2  ) 

CALL 

PLOT(l. 0,1. 0,-3) 

CALL 

SCALE (XI ,16.0, 

30 , 1 

,20.0) 

CALL 

SC ALE ( Y 1 , 1 0 . 0  , 

30 , 1 

,20.0 ) 

CALL 

SCALE ( X2  *  16 • 0  * 

30 , 1 

,20.0) 

CALL 

SCALE ( Y2 , 10 • 0  , 

30 , 1 

,20.0 ) 

CALL 

SCALE (X3 , 16 .0  ,  KCPD » 1 

,20.0 ) 

CALL 

SCALE (Y3 , 10.0  , KCPD,  1 

,20.0) 

CALL 

SCALE(X4, 16.0, KCFR, 1 

,20.0 ) 

CALL 

SCALE(Y4, 10.0, KCFR, 1 

,20.0 ) 

XMIN 

=  0.0 

XDELTA  =  0.25 
YMIN  =  0.0 
YDELTA  =  0.08 
XI  (31)  =  XMIN 
X 1 ( 32 )  =  XDELTA 
Y 1 ( 3 1  )  =  YMIN 
Y 1 ( 3 2  )  =  YDELTA 

CALL  AXISIO. 0,0.0, 'PORE  VOLUMES  WATER  INJECTED', 
1-27, 16. Q*0.0,X1(31)*X1(32)» 20.0) 

CALL  AX  IS ( 0.0 ,0. 0 , ' PORE  VOLUMES  OIL  PRODUCED', 
125, 10.0, 90. 0»Y1(31)*Y1(32)» 20.0) 

X 2  (31  )  =  XMIN 
X 2 ( 3 2  )  =  XDELTA 
X3 ( KCPD+ 1  )  =  XMIN 
X 3 ( KC PD  +  2  )  =  XDELTA 
X4(KCFR+1)  =  XMIN 
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SUBROUTINE  CAPLOT...  (CONT'D) 

X4 ( KCFR+2 )  =  XDELTA 

V2  (  3 1  )  =  YMIN 

Y2(32)  =  YDELTA 

Y3 ( KCPD+1 )  =  YMIN 

Y3 ( KCPD+2 )  =  YDELTA 

Y4 ( KCFR+1 )  =  YMIN 

Y4 (KCFR+2 )  =  YDELTA 

CALL  L  INE ( XI , Y1 » 30 » 1  *  1 , 5 ) 

CALL  LINE(X2,Y2, 30*1, 1,2) 

CALL  L INE < X3 , Y3 »KCPD* 1  , 1 »  11  ) 

CALL  SYMBOL ( 1 1 *0  *3. 0 .0 • 15 ♦ ' LEGEND  * >0.0  »6 ) 

CALL  SYMBOL (10. 0*2*0. 0.10.5*0 *0.-1) 

CALL  SYMBOL (10. 0.1. 7.0. 10. 2.0. 0,-1) 

CALL  SYMBOL ( 10.0  *1.4*0.10.11 *0.0, -1 ) 

CALL  SYMBOL( 10.0 . 1.1. 0,10. 3. 0.0* “1) 

CALL  SYMBOL( 10.5*2.0*0.10* 'DOUGLAS  AND  PLACEMANS  DATA' 
1*0.0*26) 

CALL  SYMBOL ( 10. 5 » 1.7*0. 10 *' HIGGINS  AND  LEIGHTONS 
1  RESULTS  '  *0.0*29  ) 

CALL  SYMBOL ( 10. 5 » 1 .4 *0. 10 » ' CONST. PRESS  DROP- 
1  PROPOSED  METHOD '  *0.0*32) 

CALL  SYMBOL( 10.5 *1.1*0.10* 'CONST. FLOW  RATE- 
1  PROPOSED  METHOD'  *0.0*31  ) 

CALL  SYMB0H4. 0,10. 2*0. 20,  'PLOT  OF  OIL  RECUV 
1ERY  VS  PORE-VOL- INUECTED '  *0.0  ,41 ) 

WR I TE (6  *50  ) 

50  FORMAT (' 0 ', 10X, ' AT  THE  MIDDLE  OF  THE  SUBROUTINE') 
CALL  PLOT (23.0*-1.0,-3) 

CALL  PLOT (  0.0*  0.0  ,  3  ) 

CALL  PLOT (  0.0,12.0*  2  ) 

CALL  PLOT ( 18.0*12.0*  2  ) 

CALL  PLOT ( 18.0*  0,0*  2) 

CALL  PLOT (0.0»0.0*2) 

CALL  PLOT( 1.0*1. 0,-3) 

CALL  SCALE (X5*16.0»  30,1,20.0) 

CALL  SCALE(Y5*10.0*  30,1*20.0) 

CALL  SCALE (X6*16.0»KCFR*1 *20.0) 

CALL  SCALE(Y6*10.0,KCFR,1  ,20.0) 

CALL  SCALE (X7* 16.0 »KCPD » 1 *20.0) 

CALL  SCALE ( Y 7 , 10. 0»KCPD> 1 *20.0 ) 

XMIN  =  0.0 
Y^IN  =  0.0 
XDELTA  =  70.0 
YDELTA  =  5.0 
Y  5 ( 3 1  )  =  YMIN 
Y5 ( 3 2  )  =  YDELTA 
X5  (  3  1  )  =  XMIN 
X  5 ( 3  2 )  =  XDELTA 

CALL  AX IS(0. 0*0.0* 'TIME  IN  SECS  ( *10**X )  '  ,-21 , 
116.0»0.0*X5(31)  *X5(32)  ,20.0) 
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SUBROUTINE  CAPLOT...  (CONT'D) 


CALL  AX IS(0. 0*0.0* 'WATER-OIL-RATIO' *15*10.0*90.0* 
1 Y  5  <31  )  *  Y5 ( 32  )  »  20 • 0 ) 


X6 ( KCFR+1  ) 
X6 ( KCFR  +  2  ) 
X  7 ( KCPD+ 1 ) 
X7 ( KCPD+2 ) 

Y  6 ( KCFR  +  1  ) 

Y  6 ( KCFR  +  2 ) 
Y7 ( KCPD+1 ) 
Y7 ( KCPD+2 ) 
CALL 
CALL 


CALL 

CALL 

CALL 

CALL 

CALL 

CALL 


XMIN 
XDELTA 
XMIN 
XDELTA 
YMIN 
YDELTA 
YMIN 
YDELTA 
LINE(X5»Y5*3Q»1*1»2) 

LINE (X6*Y6»KCFR*1 *1  *  11 ) 

LINE (X7  *  Y7  »  KCPD  *  1 » 1 *3) 

SYMBOL (12.0*3.0*0.15*' LEGEND ' *  0 • 0  *  6 ) 

0*0. 10*2*0. Q*-l) 

6*0. 10* 11*0. 0,-1 ) 

2*0.10 *3  *0.0  *-l  ) 

0*0. 10 » 'HIGGINS  AND  LEIGHTONS 


SYMBOL ( 11.0*2 
SYMBOL ( 11.0*1 
SYMBOL (11.0*1 
SYMBOL ( 11.5*2 
Ire sul t s  *x  =  1 • »o .0 *36 j 
CALL  SYMBOL! 11.5*1.6*0.10* 

1  ED  METHOD  *X  =  2 '  *0.0  *38  ) 

CALL  SYMBOL (11.5*1.2*0.10*' CONST 
1  ED  METHOD  *X  =  1 '  *0.0  *38 ) 

CALL  SYMBOL! 5. 0*10. 5*0. 20* 'PLOT  OF 
1-OIL-RATIO'  *0.0*31  ) 

CALL  PLOT ( 24.0»-l .0 *-3 ) 

CALL  PLOTfO. 0*0. 0*999) 

RETURN 

END 


CONST.  FLOW  RATE'-PROPOS 


PRSS  DROP-PROPOS 


TIME  VS  WATER 


aIC 


Ac  l  - 


;  •  i 


T 

•  •  •  J 


J  S  A  D 


i  r  •  ie 


:  I t  1  C  r  ;•  -  J  i 


-  1  T  A  '  ,  .  t 

(  u  •  f 

:  x 
ATJ  if  • 
I  X 
AT  J  •<  < 
I  Y 
A  T  J  -  Y 
I  ‘  Y 
AT jjc V 
(  S  .  I  *  ..  £  t  d  Y 


.  j  a  i  x  j  j  /  ^ 
if.:  )  <? v  •  (  I  )  3  Y  i 
( mm  1 )  dx 
(  T  +  '■!  d > ;  A  v 
(  I  +  O  v  )  v  X 
(  S  +  GrO>  )  rX 
(  I  +  P-3  D  •  J  3  Y 

(  s  +  s  o  ) 

(i+CPD  )w 
(  S  +  '  D )  v  y 
3 ! '  I  J  J  J  A  D 
J  J  AD 
J  J  A  D 


«  f*X  ) 

(  I  -  *  i  *  1 1  3  ’  «  d  Y  *  c  )  3  I  J 
(  ♦  X  «  I »  -  D ,  t  '  Y  t  VX  )  3 :  I  J 


'  1.  i .  1  t  c:  1  .  • »  .  •• «  .  [  xa  1  JAD 


!  *  ”  t  •  J  * 

(  1  “  ♦  •»  *  «  ^ 

(  I  —  *  (  .  » 

1 

-  —  3  T  /' 

j-.To,  D 1 

i  - '  '  U  /  C 

•  T  i'  D' 

:TA.  Is!  :  I  F  3<  1  Jjq*  « 


. .  »  .  . : i )  jcrr-'Yr  j jad 

* • - » i.if c.xx ) j  Yd  jad 
.  «S.I»  • i I )  MY3  . 

I .  «  . »  .  i ; )  j 

( J  £  t  0 •  «  * i  =  X*  £T  JU23C I 

X#Of d»X«  c«XI ) JO  Y 

rmc.c  *  •  s  =  x »  a.;  ■  t 3"  cm 

>X.J*S.X*  .  1 1  )  JO  -■  '  Yc:  J  J  AD 

(  *  -  C  *  1  I  -  At  c>T3r  031 

S  •  -  t  c?  •  1 1  J.  3  )  JOT;  •  Yc  JJAD 
(  i  *  .  *  ♦ 0 1 T  '  -JI  -I 

-•  • . -»  .  ; 

(  »  .  «  .  )  1 


'3 


APPENDIX 
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Laboratory  data  on  a  five  spot  pattern 
reported  by  Douglas,  Peaceman,  and  Rachford. 

Calculated  results  using  Higgins-Leighton 
method . 

Calculated  results  using  the  modified  scheme. 
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HIGGINS  -  LEIGHTON  RESULTS 
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CONSTANT 

PRESSURE 

DROP 

PROPOSED  METHOD 

OIL  RECOVERY 

PORE  VOL. 

I  N  j  • 

W.O.R.  PRESS.  DROP 

0. 38406E 

00 

0.36499E 

00 

0.36 

0.900E 

01 

0. 3886  7E 

00 

0.39165E 

00 

0.49 

0.900E 

01 

0.39320E 

00 

0. 39844E 

00 

0.50 

0.900E 

01 

0.39770E 

00 

0.40523E 

00 

0.51 

0.9Q0E 

01 

0.40216E 

00 

0.4U99E 

00 

0.52 

0.900E 

0  1 

0.40658E 

00 

0.41873E 

00 

0.52 

0.900E 

01 

0.41097E 

00 

0.42543E 

00 

0.53 

0.900E 

01 

0.41532E 

00 

0.43212E 

00 

0.53 

0.900E 

01 

0.41973E 

00 

0.43889E 

00 

0.53 

0.900E 

01 

0. 42422 E 

00 

0 • 44644E 

00 

1.10 

0.900E 

01 

0 • 428  05  E 

00 

0.45435E 

uo 

1 .38 

U.900E 

01 

0.431 17E 

00 

0 • 46 1 8  3  E 

00 

1.41 

0.900E 

01 

0.43422E 

00 

0.46928E 

00 

1.43 

0.900E 

01 

0.43728E 

00 

0.47675E 

00 

1.46 

0 . 900E 

01 

0.440 30 E 

00 

0.48422E 

00 

1.48 

0.900E 

01 

0 • 443  30  E 

00 

0.49167E 

00 

1.50 

0.900E 

01 

0.44626E 

00 

0.4991  IE 

00 

1.53 

0.9C0E 

01 

0.44919E 

00 

0.50653E 

00 

1.55 

0.900E 

01 

0.45209E 

00 

0.51394E 

00 

1.57 

0.900E 

01 

0.45491E 

00 

0.52130E 

00 

1.60 

0.900E 

01 

0.45774E 

00 

0.52868E 

00 

1.62 

0.900E 

01 

0.46053E 

00 

0.53604E 

00 

1.64 

0.900E 

01 

0 • 463  29  E 

00 

0.54338E 

00 

1.66 

0 . 900E 

01 

0.46604E 

00 

0.55074E 

uO 

1.69 

0 . 900  E 

01 

0.46877E 

00 

Q.55808E 

00 

1.71 

0.900E 

01 

0.47145E 

00 

0. 56540E 

00 

1.73 

0.900E 

01 

0.4741QE 

00 

0.57269E 

00 

1.75 

0.900E 

01 

0.47678E 

00 

0. 58004E 

00 

1.77 

0.900E 

01 

0.47944E 

00 

0.58736E 

00 

1.79 

0.900E 

01 

0.48201E 

00 

0.59462E 

00 

1 .80 

0.900E 

01 

0.48459E 

00 

0.60190E 

00 

1.82 

0.900E 

01 

0.48717E 

00 

0.60920E 

00 

1.84 

0.900E 

01 

0.48976E 

00 

0.61652E 

00 

1.85 

0.900E 

01 

0.49232E 

00 

0.62381E 

00 

1.86 

0.900E 

01 

0.49486E 

00 

0*631 10E 

00 

1.86 

0.900E 

01 

0.4-9747E 

00 

0.63847E 

00 

1.79 

0 . 900E 

01 

0. 50023E 

00 

0 . 6460UE 

00 

1*64 

0.900E 

01 

0. 50270E 

00 

0.65385E 

00 

3.11 

0.900E 

01 

0 •  5  04  5  5  E 

00 

0.66184E 

00 

3.81 

0.900E 

01 

0.50618E 

00 

0.66977E 

00 

3.86 

0.900E 

01 
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CONSTANT  PRESSURE 
OIL  RECOVERY  PORE  VOL. 


0.50777E 

00 

0.67767E 

0. 50933E 

00 

0.68556E 

0.51094E 

00 

0.69352E 

Q.51255E 

00 

0.70149E 

0. 51413E 

00 

0.70945E 

0. 51564E 

00 

0.71733E 

0.51715E 

00 

0.72527E 

0.51872E 

00 

0.73326E 

0.52020E 

00 

0.741 17E 

0.52168E 

00 

0.74910E 

0. 52317E 

CO 

Q.75705E 

0. 52463E 

00 

0.76499 E 

0.52610E 

0  0 

0.77294E 

0.52752E 

00 

0.78C36E 

0.52900E 

00 

0.78886E 

0. 53G48E 

00 

0.79688E 

0. 53187E 

00 

0.80482E 

0. 53322E 

00 

G.81274E 

0.53458E 

0  0 

0.82068E 

0.53594E 

00 

0.82863E 

0.53731E 

00 

0.83661E 

0. 53868E 

00 

0 . 8446 OE 

C. 54002E 

00 

G.85256E 

Q.54136E 

00 

G.86054E 

0. 54265E 

00 

G.86847E 

0.54394E 

00 

0.87643E 

0.54523E 

0  0 

0.88439E 

0. 54654E 

00 

0.89238E 

0. 54785E 

00 

0.90039E 

0. 549G5E 

00 

0.90830E 

0.55031E 

00 

G.91627E 

0. 55163E 

00 

0.92431E 

0. 55285E 

00 

0.93226E 

0*  55407E 

00 

0.94023E 

0. 55530E 

00 

0.94821E 

0.55653E 

00 

0.95620E 

0. 55777E 

00 

0 . 9642  IE 

0.55900E 

00 

0.97224E 

0.3602UE 

00 

0.98022E 

0.56133E 

00 

0.98813E 

DROP  PROPOSED  METHOD 


INj. 

W  •  0  .  R  • 

PRESS.  DROP 

00 

3.91 

0.900E 

01 

00 

3.96 

0.900E 

01 

00 

4.01 

0.900E 

01 

00 

4.06 

0.9C0E 

Cl 

00 

4.11 

0.900E 

01 

00 

4.16 

0.900E 

01 

00 

4.21 

0.900E 

01 

00 

4.26 

0.900E 

01 

00 

4.31 

0.900E 

01 

00 

4.36 

0.900E 

01 

00 

4.41 

0 • 900E 

01 

00 

4.46 

0.900E 

01 

00 

4.51 

0.900E 

01 

00 

4.56 

0.900E 

01 

00 

4.62 

0.900E 

01 

00 

4.67 

0.900E 

0  1 

00 

4.71 

0.900E 

01 

00 

4.76 

0.900E 

01 

00 

4.80 

U.900E 

01 

00 

4.8  5 

0.900E 

01 

00 

4.90 

0.90CE 

01 

00 

4.95 

0.900E 

01 

00 

4.99 

0.900E 

01 

00 

5.04 

0.900E 

01 

00 

5.09 

0.900E 

01 

00 

5.13 

0.900E 

01 

00 

5.18 

0.900E 

01 

00 

5.23 

0.900E 

01 

00 

5.28 

0.900E 

01 

00 

5.32 

0 . 9  0  0  E 

01 

00 

5.36 

0  •  v  0  0  E 

01 

00 

5.41 

0.900E 

Cl 

00 

5.46 

0 .900E 

01 

00 

5.50 

0.900E 

01 

00 

5.55 

0.900E 

Cl 

00 

5 . 59 

0.900E 

01 

oc 

5.64 

0.900E 

01 

00 

5.69 

0 .900E 

01 

00 

5.73 

0.900E 

01 

00 

5.77 

0.900E 

01 
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C  -  6 

CONSTANT  PRESSURE  DROP  PROPOSED  METHOD 


OIL  RECOVERY 

PORE  VOL. 

• 

“> 

2! 

i — i 

W  •  0  •  R  • 

PRESS,  i 

DROP 

0.56252E 

00 

0.99615E 

00 

5.81 

0.900E 

01 

0.56372E 

00 

0. 10042E 

01 

5.85 

0 • 900E 

01 

0*  56493E 

00 

0. 10122E 

01 

5.89 

0.900E 

01 

0 •  5  660  8  E 

00 

0. 10202E 

01 

5.94 

0.900E 

01 

0.56717E 

00 

0. 10281E 

01 

5.97 

0.900E 

01 

0. 56832E 

00 

0. 10361E 

01 

6.02 

0.900E 

0  1 

0.56948E 

00 

0. 10441E 

01 

6.07 

0.900E 

01 

Q.57065E 

00 

0. 103 2  IE 

01 

6.12 

0.900E 

01 

0. 57174E 

00 

0. 10600E 

01 

6  •  1 6 

0.900E 

01 

0. 57279E 

00 

0. 1U680E 

01 

6.19 

0 .900E 

01 

Q.57397E 

00 

0. 10760E 

01 
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0.900E 

01 

0.57515E 

00 

0. 10641E 
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0.900E 
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00 
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0.900E 
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00 
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C  -  9 

CONSTANT  INJECTION  RATE  PROPOSED  METHOD 


OIL  RECOVERY 

PORE  VOL. 

INj. 

W .  0 .  R . 

PRESS. 

DROP 

0 • 38402  E 

00 

0. 38494E 

00 

0.36 

0.818E 

00 

0 . 3  8  8  6  2  E 

00 

0.39160E 

00 

0.49 

0.812E 

00 

Q.39314E 

00 

0.39839E 

00 

0.50 

0  •  807E 

00 

0.39765E 

0  0 

0.40518E 

00 

0.51 

0 • 80 1 E 

00 

0.4Q211E 

00 

C.41194E 

00 

0.52 

0.795E 

00 

0.40653E 

00 

0.41867E 

00 

0.52 

0. 789E 

00 

0.41091E 

00 

0.42538E 

00 

0.53 

0. 783E 

00 

0.41526E 

00 

0.43207E 

00 

C.  54 

0. 776E 

00 

0.41967E 

00 

0.43883E 

00 

0.53 

0. 755E 

00 

0.42417E 

00 

0.44639E 

00 

1.10 

0.706E 

00 

Q.42902E 

00 

0.4543  IE 

00 

1 .38 

0.670E 

00 

0.431 14E 

00 

0.46180E 

00 

1.41 

0.665E 

00 

0.43419E 

0  0 

Q.46924E 

00 

1.43 

0.660E 

00 

0.43725E 

00 

0.47672E 

00 

1.46 

0.656E 

00 

0.44027E 

00 

0.48418E 

00 

1.48 

0.651E 

00 

0.44327E 

0  0 

0.49163E 

00 

1.50 

0.647E 

00 

0.44623E 

00 

Q.499Q7E 

00 

1 .53 

0.642E 

00 

0.44912E 

00 

0. 50646E 

00 

1.55 

o  .  6  3  8  E 

00 

0.45202E 

00 

0.51387E 

00 

1.57 

0.634E 

00 

0.45488E 

00 

0.52126E 

00 

1.60 

0.629E 

00 

0.45771E 

00 

0.52864E 

00 

1.62 

0.625E 

00 

0.46050E 

00 

0.53600E 

0  0 

1.64 

0.621E 

00 

0.46322E 

0  0 

0.54330E 

00 

1.66 

0.617E 

00 

0.46597E 

00 

0. 55066E 

00 

1.69 

0.613E 

00 

C.46873E 

00 

0 . 5  5  6  0  3  E 

00 

1.71 

0.609E 

00 

0.47 142 E 

00 

0.56534E 

00 

1.73 

0.606E 

00 

0.474u6E 

00 

0. 57264E 

00 

1.75 

0.602E 

0  0 

0.47671E 

00 

0.57994E 

00 

1.77 

0.598E 

00 

0.47936E 

00 

0.58727E 

00 

1 . 79 

0.595E 

00 

0.48197E 

00 

0.59457E 

00 

1.80 

0.591E 

0  0 

0.48455E 

00 

J.60134E 

00 

1.82 

LU 

r~- 

u> 

. 

00 

0.48714E 

00 

0.60914E 

00 

1.84 

0 . 5  84  E 

00 

0.48973E 

00 

0 .61645E 

00 

1.85 

0  •  5  8  0  E 

00 

C.49228E 

00 

0.62375E 

00 

1.8  6 

0.576E 

00 

0.49482E 

00 

0.63103E 

00 

1.86 

J . 570E 

00 

0.49743E 

00 

0.63840E 

00 

1.79 

0 . 5  5  8  E 

0  0 

0. 50018E 

00 

0.64593E 

00 

1.65 

0.536E 

00 

0.50267E 

00 

0.65379E 

00 

3.11 

0.518E 

00 

0.50454E 

00 

0.66180E 

00 

3.81 

0.512E 

00 

C.  50617E 

0  0 

0.66972E 

00 

3.86 

0.509E 

00 
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CONSTANT 

INJECTION 

RATE 

PROPOSED  METHOD 

TIME 

OIL  RECOVERY 

PORE  VOL. 

I  Nj . 

W.O.R.  PRESS.  DROP 

0. 21354E 

05 

0. 50776E 

00 

0.67763E 

00 

3.91 

0.507E 

00 

0.21604E 

05 

0.50932E 

0  0 

0.68552E 

00 

3.96 

0.504E 

00 

0.21855E 

05 

0.51093E 

00 

0.69348E 

00 

4.01 

0. 501E 

00 

0.22105E 

05 

0. 5125QE 

00 

0.70141E 

00 

4.06 

0.499E 

00 

0.22356E 

05 

0. 51409E 

00 

0.70937E 

00 

4.11 

0.496E 

00 

0.22606E 

05 

C.  51563E 

00 

0.71730E 

00 

4.16 

0.494E 

00 

0.22856E 

05 

0.51714E 

00 

0.72522E 

00 

4.21 

Q.491E 

00 

0.23107E 

05 

0.51871E 

00 

0.73321E 

00 

4.26 

0.489E 

00 

0.23357E 

0  5 

0.52018E 

00 

0.74112E 

00 

4.31 

0 • 48  6  E 

00 

0.23608E 

05 

0.52167E 

00 

0.74905E 

00 

4.36 

0.484E 

00 

0.23858E 

05 

0. 523  16E 

00 

0.75700E 

00 

4.41 

0.482E 

00 

C.24108E 

05 

0. 52462E 

00 

0.7649 3 E 

00 

4*46 

0.479E 

00 

0. 24358E 

05 

0. 52609E 

00 

0.772S9E 

00 

4.51 

0.477E 

00 

0  •  24609 E 

05 

0 . 5  27  5 1 E 

0  0 

0.78081E 

00 

4.56 

0 .475E 

00 

0. 24859E 

05 

0. 52894E 

00 

0.78876E 

00 

4.61 

0.473E 

00 

0.25110E 

05 

0. 53042E 

00 

0.79678E 

00 

4.67 

0.471E 

00 

0. 25361E 

05 

0. 53181E 

00 

0.80472E 

00 

4.71 

0.469E 

00 

0. 25612E 

05 

0.53316E 

0  0 

0.81263E 

00 

4.75 

0 . 4  6  7  E 

00 

0. 25863E 

05 

0. 5 3457 E 

00 

0 . 8  2  0  6  2  E 

00 

4.80 

0.465E 

00 

0.26114E 

05 

0. 53593F 

00 

Q.62857E 

00 

4.85 

0.463E 

00 

0.26365E 

05 

0. 53730E 

00 

0.83654E 

00 

4.90 

0.461E 

00 

0.2661AE 

05 

0. 53863E 

00 

0  •  8444  9 E 

00 

4.94 

0.459E 

00 

0.26867E 

05 

0. 53996E 

00 

0.P5245E 

00 

4.99 

0.457E 

00 

0.27119E 

05 

0. 54134E 

00 

Q.86047E 

00 

5.04 

C.455E 

00 

0. 27370E 

05 

0. 54263E 

00 

0.P6840E 

00 

5.09 

0 . 4  5  3  E 

00 

0  •  2762  IE 

05 

0. 54392E 

00 

0.87635E 

00 

5 .13 

0.452E 

00 

0.27872E 

05 

0. 54522E 

00 

0.88432E 

00 

5.18 

0.450E 

00 

0. 28124E 

05 

0.54652E 

0  0 

0.89231E 

00 

5.23 

0 . 448  E 

00 

0.28375E 

05 

0. 54783E 

00 

C.9Q031E 

00 

5.28 

0 • 446  E 

00 

0 . 28627E 

05 

0. 54904E 

00 

0.9Q822E 

00 

5 .32 

0 . 445  E 

00 

0. 28879E 

05 

0.55025E 

00 

0 .91615E 

00 

5.36 

0.443E 

00 

0.29130E 

05 

0. 55157E 

00 

0.92419E 

CO 

5.41 

0 • 44 1 E 

CO 

0.29381E 

05 

0. 55284E 

0  0 

0.93219E 

00 

5.46 

0.440E 

0  0 

0.29633E 

05 

0. 55406E 

00 

0.94015E 

00 

5 . 50 

0.438E 

00 

0.29884E 

05 

0.55529E 

00 

0.94813E 

00 

5.55 

0.437E 

00 

0.30136E 

05 

0.55652E 

00 

0 . 9  5  6  1 2  E 

00 

5.59 

0.435E 

0  0 

0 . 30  3  88  E 

05 

0. 55775E 

00 

0.96413E 

00 

5 . 64 

Q.433E 

00 

0.30640E 

05 

0. 55894E 

00 

0.97210E 

00 

5.69 

0.432E 

0  0 

0. 3089  IE 

05 

0.56013E 

00 

0.98009E 

00 

5.73 

0 . 4  3  0  E 

00 

0. 31 143E 

05 

0.56125E 

00 

0.98801E 

CO 

5.77 

0.429E 

00 
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C  -  11 

CONSTANT  INJECTION  RATE  PROPOSED  METHOD 


OIL  RECOVERY 

PORE  VOL. 

I  N  J  • 

W .  0  •  R  . 

PRESS.  1 

DROP 

0. 56245E 

00 

0.99600E 

00 

5.81 

0.428E 

00 

0  •  5  6365 E 

00 

0. 1 0  0  4  0  E 

01 

5.85 

0 .426E 

00 

0. 56486E 

00 

0. 10120E 

01 

5 .89 

C.425E 

00 

0. 56606E 

00 

0. 1Q201E 

01 

5.94 

0.423E 

00 

0. 56715E 

00 

0. 10280E 

01 

5.97 

0  •  4  2  2  E 

00 

0. 56831E 

00 

0. 10350E 

01 

6.02 

0 .420E 

00 

0 • 56947E 

00 

0. 10440E 

01 

6.07 

0 . 4 1 9  E 

00 

0. 57063E 

00 

0. 1Q520E 

01 

6.12 

0 .417E 

00 

Q.57173E 

00 

0. 10599E 

01 

6.16 

0.416E 

00 

0.57275E 

00 

0. 10679E 

01 

6.19 

0  •  4 1  5  E 

00 

0.57395E 

00 

0. 10759E 

01 

6.24 

0.414E 

00 

0. 575  13E 

00 

0. 10839E 

01 

6.28 

0.412E 

00 

0. 57619E 

0  0 

0. 10919E 

01 

6.32 

0.411E 

00 

0.57718E 

0  0 

0. 10998E 

01 

6.35 

0 . 4 1 0  E 

00 

0. 57830E 

00 

0 . 1 1 0  7  8  F 

01 

6.39 

0.409E 

00 

0. 57942E 

00 

0. 11158E 

01 

6.43 

0.407E 

00 

0. 53049E 

00 

0. 11238E 

01 

6.46 

0 • 406  E 

00 

0.58156E 

00 

0. 11318E 

01 

6.50 

0.405E 

00 

0.58257E 

0  0 

C.  11397E 

01 

6.53 

0.404E 

00 

0. 58358E 

0  0 

0. 1I476E 

01 

6.56 

0.403E 

00 

0*  58466E 

00 

0. 11556E 

01 

6.60 

0.401E 

0  0 

0. 58575E 

00 

0. 1163  7E 

01 

6.63 

0.400E 

00 

0. 58677E 

00 

0. 11716E 

01 

6.66 

0. 399E 

00 

0.58779E 

00 

0. 11796E 

01 

6.69 

0.398E 

00 

0 . 5  8  8  8  9  E 

00 

0. 11676E 

01 

6.7  3 

C.397E 

00 

0. 58999E 

00 

0. 11957E 

01 

6.76 

0.395E 

00 

0. 59103E 

00 

0. 12037E 

01 

6.79 

0.394E 

00 

0. 59201E 

00 

0.12117E 

01 

6.81 

0.393E 

00 

0. 59300E 

00 

0. 12196E 

01 

6.83 

0.392E 

00 

0. 59404E 

00 

0. 12277E 

01 

6.85 

0.391E 

0  0 

0.59518E 

00 

0. 12358E 

01 

6*74 

0.388E 

00 

0. 59629E 

00 

0. 12439E 

01 

6.48 

0.385E 

00 

0. 59738E 

00 

0. 1252UE 

01 

6.22 

0.381E 

00 

Q.59852E 

00 

0. 12601E 

01 

5.88 

0.377E 

00 
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